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USER  AND  ANALYST  MANUAL  FOR  A  FORTRAN  COMPUTER  PROGRAM 
SIMULATING  THE  ENGAGEMENT  OF  A  STATIONARY  POINT  TARGET 
BY  A  STATIONARY  DIRECT  FIRE  WEAPON 


1.  INTRODUCTION 


This  report  describes  a  computer  program  developed  under  the 
sponsorship  of  the  Joint  Munitions  Effectiveness  Manual  Surface-to-Surface 
(JMEM-SS)  Methodology  and  Evaluations  Working  Group,  with  a  substantial 
contribution  by  the  Army  Materiel  Systems  Analysis  Activity  (AMSAA).  The 
program  has  been  used  by  JMEM-SS  and  AMSAA  to  make  effectiveness  calculations 
for  armor  and  infantry  weapons  engaging  tanks,  armored  personnel  carriers, 
and  other  targets.  It  utilizes  Monte  Carlo  techniques  to  simulate  the 
engagement  and  has  been  based,  to  a  large  extent,  on  a  similar  deterministic 
program  previously  developed  by  AMSAA.  An  earlier  version  of  the  JMEM-SS 
direct  fire  program  has  been  documented  in  1976  by  the  then  US  Army 
Armament  Command  (Reference  1). 


The  JMEM-SS  direct  fire  program  provides  a  method  for  calculating 
the  effectiveness  of  weapons  that  fire  aimed  rounds  one  at  a  time  against 
point  targets.  Such  weapons  include  tank  main  guns  and  infantry  antitank 
weapons.  The  required  input  includes  performance  parameters  related  to 
the  rate  of  fire,  delivery  accuracy,  and  reliability  of  the  firing  weapon. 
A  detailed  target  description  and  data  indicating  vulnerability  of  the 
target  to  the  type  of  round  fired  are  also  part  of  the  input.  Other 
influencing  factors  specified  for  a  particular  calculation  are  the  range, 
the  orientation  of  the  target,  whether  the  target  is  fully  exposed  or 
only  partly  exposed,  and  the  aimpoint.  To  simulate  a  single  engagement, 
the  computer  determines  the  effects  of  the  first  round,  taking  into 
account  several  kill  criteria  of  interest.  Then,  as  necessary,  second 
and  succeeding  rounds  up  to  a  limit  of  15  rounds  are  evaluated.  Although 
this  limit  is  a  practical  one  for  the  calculations,  it  clearly  exceeds 
the  number  of  rounds  normally  fired  in  real  combat  engagements.  The 
calculations  for  the  first  engagement  result  in  an  actual  possible  outcome 
being  established.  For  example,  the  target  might  be  determined  to  be  hit 
on  the  third  round  38  seconds  after  the  start  of  the  engagement.  Similar 
specific  numerical  values  are  also  assigned  to  the  round  number  and  time 
correspond!' ng  to  the  target  being  killed  for  each  kill  criterion  of 
concern.  The  process  is  repeated  until  10,000  engagements  have  been 
considered,  and  data  for  all  engagements  are  combined  to  give  plausible 
outputs  not  unduly  influenced  by  particular  engagements.  Outputs 
are: 


Michels,  Herman  W.,  An  Analyst/User  Manual  for  the  Direct  Fire  Computer 
Program,  DRSAR/SA/R-18,  May  1976,  US  Army  Armament  Command,  Rock  Island, 
IL,  61201. 
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a.  Hit  probabilities  and  kill  probabilities  corresponding  to 
various  numbers  of  rounds  fired. 

b.  Average  numbers  of  rounds  needed  to  hit  or  to  kill  the 

target. 


c.  Hit  probabilities  and  kill  probabilities  for  various  engagement 
times  up  to  a  practical  limit,  established  as  2  minutes  for  the  calculations. 

The  engagement  simulated  by  the  direct  fire  program  is  restricted  to 
being  one-sided.  This  means  that  the  firing  weapon  is  not  subjected  to 
return  fire  from  the  target  or  some  other  opposing  element. 


This  report  has  been  prepared  primarily  for  users  and  analysts 
with  an  interest  in  documentation  far  exceeding  in  scope  the  very  general 
explanations  in  this  introduction.  The  material  included  is  as  follows: 

a.  Basic  definitions  (section  2). 

b.  General  flow  of  program  logic  (section  3). 

c.  Input  data  (section  4). 

d.  Output  data  (section  5). 

e.  Computer  aspects  (section  6). 

f.  Program  listing  and  detailed  explanations  of  program  statemerts 
(Appendixes  A  and  B). 


Input  data  used  to  exercise  the  JMEM-SS  direct  fire  program  as 
well  as  outputs  of  the  calculations  frequently  require  a  CONFIDENTIAL 
classification.  To  make  this  user  and  analyst  manual  UNCLASSIFIED,  input 
and  output  data  not  bearing  any  identifiable  relation  to  classified 
information  have  been  fabricated  to  illustrate  inputs  and  outputs  associated 
with  the  program. 


The  direct  fire  program  contains  some  options  that 
or  no  concern  to  most  users  or  analysts.  These  options  are 
this  report  but  not  necessarily  discussed  thoroughly.  If  a 
detailed  explanations  concerning  any  program  option  arises, 
provided  by  AMSAA  upon  request. 


are  of  little 
identified  in 
need  for 
these  can  be 
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2.  BASIC  DEFINITIONS 
2.1  Introduction. 

This  section  contains  definitions,  as  well  as  related  assumptions 
and  explanations,  intended  to  clarify  the  meaning  of  terms  used  in  this 
report.  The  following  terms  are  covered: 


K  Kill  (2.22) 

M  Kill  (2.23) 

F  Kill  (2.24) 

M  Or  F  Kill  (2.25) 

Kill  Probability  for  Specified 


Firing  Engagement  (2.2) 
Engagement  Time  (2.3) 
Orientation  (2.4) 

Exposure  (2.5) 

Reliability  (2.6) 

Delivery  Accuracy  (2.7) 

Hit  Probability  (2.8) 

Firing  Occasion  (2.9) 

Fixed  Bias  (2.10) 

Variable  Bias  (2.11 ) 

Random  Error  (2.12) 

Lay  Error  (2.13) 
Round-to-Round  Error  (2.14) 
Sensing  (2.15) 

Sensing  Probability  (2.16) 
Sensing  Error  (2.17) 

Target  Vulnerability  (2.18) 
Kill  Probability  (2.19) 

Kill  Criterion  (2.20) 
K,M,F,M  Or  F  (2.21) 


Number  of  Rounds  Fired  (2.26) 
Average  Number  of  Rounds  Needed 
to  Kil  1  Target  (2.27) 

Expected  Personnel  Casualties  (2.28) 
Rate  of  Fire  (2.29) 

First  Round  Firing  Time  (2.30) 

Flight  Time  (2.31) 

Subsequent  Round  Firing  Time  (2.32) 
Fixed  Time  (2.33) 

Variable  Time  (2.34) 

Median  Time  (2.35) 

Variability  Factor  (2.36) 

Minimum  Time  (2.37) 

Kill  Probability  for  Specified 


Engagement  Time  (2.38) 


2.2  Firing  Engagement. 

The  firing  of  one  or  more  rounds  at  a  target.  Although  combat 
is  generally  two-sided,  engagements  simulated  in  the  direct  fire  program 
are  only  one-sided  in  that  the  firing  weapon  is  not  itself  subjected  to 
any  fire.  It  is  assumed  that  neither  the  target  nor  the  firing  weapon 
changes  its  position  during  the  engagement.  An  engagement  begins  when 
the  decision  is  made  to  fire  at  the  target  with  a  particular  type  of 
round.  The  engagement  ends  when  either  the  target  is  killed  or  15  rounds 
have  been  fired.  This  limit  has  been  established  for  the  calculations  to 
ensure  obtaining  all  data  of  possible  interest.  It  is  realized  this  limit 
exceeds  the  number  of  rounds  available  in  the  ready  rack  of  tanks.  Firing 
engagements  in  combat  could  sometimes  be  deliberately  ended  for  reasons 
not  represented  in  the  computer  simulation  program.  For  example,  an 
armored  force  might  have  a  policy  whereby  each  tank  fires  no  more  than 
some  specified  number  of  rounds,  like  three,  at  any  particular  target  in 
certain  situations. 
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2.3  Engagement  Time. 


» 


Time  elapsing  from  the  beginning  of  the  firing  engagement. 
Engagement  time  increases  as  the  program  calculations  account  for  more 
and  more  specific  events,  for  example,  the  firing  of  the  first  round  and 
the  arrival  of  this  round  at  the  target  range.  Engagement  time  can,  but 
does  not  always,  mean  the  time  needed  for  the  entire  engagement  to  be 
completed. 

2.4  Orientation. 

The  view  that  the  target  presents  to  the  firing  weapon. 

Orientation  is  specified  as  an  angle  with  0  degrees  representing  the 
front  view,  180  degrees  the  rear  view,  and  90  and  270  degrees  the  side 
view  from  two  opposite  directions.  Other  orientations  are  represented 
by  angles  from  0  to  180  degrees  or  from  180  to  360  (same  as  0)  degrees. 

An  observer  facing  the  target  and  located  on  the  circumference  of  a  circle 
with  the  target  at  its  center  should  move  to  the  right  to  see  the  target 
at  orientations  of  0  to  180  degrees,  and  to  the  left  for  orientations 
exceeding  180  degrees. 

2.5  Exposure. 

Whether  the  maximum  possible  target  area  or  only  a  portion  of 
the  maximum  area  is  presented,  for  a  particular  orientation,  to  the  firing 
weapon.  Two  exposures  are  normally  considered  in  applications  of  the 
direct  fire  program: 

a.  Open  Exposure:  The  target  presents  its  maximum  presented 
area  for  a  particular  orientation. 

b.  Defilade  Exposure:  This  condition  has  been  applied  only  to 
tank  targets.  Defilade,  or  hull  defilade,  applies  to  the  condition  where 
a  tank  target  presents  only  its  turret  as  the  target.  The  orientation 
specified  applies  to  the  turret  and  not  the  hull  since  the  hull  is  shielded. 

2.6  Reliability. 

A  measure  of  the  ability  of  a  round  to  function  properly  when  it 
is  fired  and  of  the  fuze  to  function  properly  upon  impact.  Proper 
functioning  basically  refers  to  the  absence  of  catastrophic  malfunctions 
affecting  flight  of  the  projectile  or  fuze  action. 

2.7  Delivery  Accuracy. 

A  measure  of  the  ability  of  rounds  to  hit  a  target.  Delivery 
accuracy  data  applying  to  first  rounds  and  data  for  subsequent  rounds  are 
both  needed  for  the  direct  fire  program  calculations. 
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2.8  Hit  Probability. 

The  chance  that,  under  specific  conditions,  a  target  is  hit. 

2.9  Firing  Occasion. 

The  firing  of  a  small,  continuous  sequence  of  rounds  at  a 
particular  target. 

2.10  Fixed  Bias. 

Errors  of  the  weapon-ammunition-fire  control  system  which  are 
usually  constant  and  predictable  at  any  given  range,  and  are  fixed  for 
all  firing  occasions  by  the  system  design.  Although  the  term  fixed  bias 
may  refer  to  an  individual  error,  it  is  used  in  this  report  only  for  the 
aggregated  fixed  bias  error,  in  each  of  the  horizontal  and  vertical 
directions,  that  accounts  for  all  applicable  individual  error  sources. 

2.11  Variable  Bias. 

Errors  whose  values  remain  very  nearly  constant  during  a  particular 
firing  occasion,  but  which  may  vary  considerably  from  occasion  to  occasion. 
Term  is  used  in  this  report  only  for  the  aggregated  horizontal  or  vertical 
variable  bias  accounting  for  all  contributing  error  components. 

2.12  Random  Error. 

Errors  or  an  error  whose  magnitude  and  direction  vary  randomly 
from  round  to  round  even  during  a  single  firing  occasion.  The  direct 
fire  program  explicitly  involves  the  lay  error  and  the  round-to-round 
error,  in  each  of  the  horizontal  and  vertical  directions,  rather  than 
aggregated  random  errors. 

2.13  Lay  Error. 

Random  error  associated  with  the  fine  lay  made  by  the  gunner 
before  firing. 

2.14  Round-to-Round  Error. 


Random  error  resulting  from  differences  between  individual 

rounds. 


2.15  Sensing. 

Process  whereby  sufficient  information  concerning  the  trajectory 
of  a  missing  round  is  obtained  by  the  gunner,  or  by  the  commander  with 
transmittal  to  the  gunner,  as  a  basis  for  adjustment  of  fire  before  firing 
of  the  next  round.  Sensing  is  applicable  to  unguided  rounds  only. 
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2.16  Sensing  Probability. 


The  chance  that  sensing  occurs  as  a  basis  for  adjustment  of  fire 
on  the  next  round. 

2.17  Sensing  Error. 

The  difference  between  the  weapon  aimpoint  associated  with  a 
missing  round  that  is  sensed  and  the  new  aimpoint  following  adjustment  of 
fire.  Such  differences  account  not  only  for  incorrect  judgments  about 
where  missing  rounds  went  and  possibly  inaccuracies  arising  from  transmittal 
of  information  from  the  commander  to  the  gunner,  but  also  for  any  errors 
made  by  the  gunner  in  relaying  as  part  of  the  fire  adjustment  process. 

The  sensing  error  in  each  of  the  horizontal  and  vertical  directions  is  of 
concern  for  the  direct  fire  program  calculations  involving  unguided  rouncs. 

2.18  Target  Vulnerability. 

A  measure  of  the  target's  susceptibility  to  sustain  a  given 
amount  of  damage  when  hit  by  a  round. 

2.19  Kill  Probability. 

The  chance  that  a  target  is  killed  as  a  consequence  of  damage 
sustained.  The  word  kill  acquires  a  definite  meaning  when  it  is  associated 
with  a  specific  kill  criterion. 

2.20  Kill  (Damage)  Criterion. 

The  function  or  functions  that  a  target  must  lose  to  be  considered 
out  of  action.  A  target  killed  in  one  engagement  could  be  repaired  and 
participate  in  a  subsequent  engagement. 

2.21  K,  M,  F,  M  Or  F. 

Symbols  representing  four  kill  criteria  whose  specific  definitions 
have  been  used,  by  the  United  States  and  several  other  countries,  as  the 
basis  for  vulnerability  analyses  of  armored  vehicles. 

2.22  K  Kill. 


A  target  vehicle  is  subjected  to  a  K  kill  (complete  destruction' 
if  it  sustains  both  an  M  kill  and  an  F  kill  and  is  damaged  to  the  extent 
that  it  is  not  economically  repairable.  A  K  kill  is  more  likely  to  be 
apparent  to  the  crew  of  the  firing  weapon  than  any  other  kill  because  of 
resulting  fires  and/or  detonation  of  ammunition. 
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2.23  M  Kill. 


A  target  is  subjected  to  an  M  kill  (mobility)  if  it  is  incapable 
of  executing  controlled  movement  and  the  damage  is  not  repairable  by  the 
crew  on  the  battlefield.  The  loss  of  this  function  may  be  caused  by 
either  incapacitation  of  the  crew  or  damage  to  the  propulsion  or  control 
equipment.  A  target  that  is  M  killed  could  sustain  other  types  of  damage 
at  the  same  time,  i.e. ,  F  kill  or  K  kill. 

2.24  F  Kill. 

A  target  is  subjected  to  an  F  kill  (firepower)  if  it  is  incapable 
of  delivering  controlled  fire  from  the  main  armament  and  the  damage  is 
not  repairable  by  the  crew  on  the  battlefield.  The  loss  of  this  function 
may  be  caused  by  either  incapacitation  of  the  crew  or  damage  to  the  main 
armament  and  its  associated  equipment.  A  target  that  is  F  killed  could 
sustain  other  types  of  damage  at  the  same  time,  i.e.,  M  kill  or  K  kill. 

2.25  M  Or  F  Kill. 

A  vehicle  is  subjected  to  an  H  or  F  kill  if  it  sustains  either 
an  M  kill  or  an  F  kill,  or  both  an  M  kill  and  F  kill.  A  target  that 
sustains  an  H  or  F  kill  could  also  at  the  same  time  sustain  a  K  kill. 

2.26  Kill  Probability  for  Specified  Number  of  Rounds  Fired. 

Probability  that  the  target  is  killed  when  a  specific  number  of 
rounds,  say  N,  are  fired.  This  probability  is  meaningful  only  in 
conjunction  with  a  kill  criterion.  N  is  considered  to  equal  1  through  15 
for  the  program  calculations.  The  kill  probability  associated  with  N 
rounds  is  an  aggregated  result  that  accounts  for  kills  caused  by  any 
round  up  to  and  including  the  Nth  round.  When  N  equals  2,  the  second 
round  is  fired  only  when  the  first  does  not  kill  the  target.  More 
generally,  any  round  after  the  first  is  fired  only  when  the  target  has 
not  already  been  killed. 

2.27  Average  Number  of  Rounds  Needed  to  Kill  Target. 

The  average  number  of  rounds  needed  to  kill  a  target  is  the 
arithmetic  average  of  the  number  of  rounds  needed  to  achieve  a  specified 
type  of  kill  given  an  unlimited  number  of  individual  firing  engagements. 
Although  engagements  are  not  continued  beyond  15  rounds,  the  program 
calculations  account  for  additional  rounds  needed  in  exceptional  instances 
where  the  target  may  survive  15  rounds.  It  is  not  implied  by  this  that 
such  large  numbers  of  rounds  would  be  involved  in  an  actual  combat 
engagement. 
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2.28  Expected  Personnel  Casualties. 


Refer  to  the  average,  or  expected,  number  of  personnel  casualties 
sustained  by  the  troops  carried  in  an  armored  personnel  carrier  target. 

Crew  members  are  excluded.  Calculations  of  personnel  casualties  have 
been  based  on  the  assault  —  5-minute  criterion,  which  means  a  man 
must  be  rendered  incapable  of  performing  an  assault  role  within  5  minutes 
after  being  wounded  in  order  to  be  considered  a  casualty. 

2.29  Rate  of  Fire. 

A  measure  of  the  ability  of  a  weapon  to  fire  aimed  rounds  at  a 
target.  The  direct  fire  program  requires  input  data  permitting  the 
computer  to  generate  the  distribution  of  times,  not  only  average  times, 
that  crews  would  need  to  fire  the  weapon  in  a  firing  engagement. 
Furthermore,  care  is  necessary  to  ensure  that  the  inputs  correspond  to 
rate  of  fire  under  combat  conditions  rather  than  test  conditions  generally 
tending  to  favor  rapid  fire. 

2.30  First  Round  Firing  Time. 

Time  elapsing  between  the  beginning  of  the  engagement  and  firing 
of  the  first  round.  Operations  involved  in  firing  of  the  first  round 
include,  as  applicable,  the  following: 

a.  Slewing  of  the  weapon. 

b.  Target  recognition  by  the  gunner. 

c.  Ranging,  which  may  be  done  with  equipment,  such  as  a 
rangefinder,  or  by  visual  estimation  of  the  range  to  the  target. 

d.  Laying  of  the  weapon. 

Note  that  some  of  these  operations  can  be  carried  out  at  the  same  time 
that  others  are.  Loading  of  the  first  round  may  be  included  or  excluded. 

In  the  latter  case,  the  appropriate  type  of  round  is  considered  already 
loaded  before  recognition  of  the  target. 

2.31  Flight  Time. 

Time  interval  between  the  time  a  round  begins  to  be  propelled 
forward  in  the  gun  or  launcher  and  the  time  it  reaches  the  range  of  the 
target. 

2.32  Subsequent  Round  Firing  Time. 

Time  that  elapses  between  firing  of  any  round  and  firing  of  the 
next  round  against  the  same  target.  Subsequent  round  firing  time  involves, 
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as  applicable,  the  following: 

a.  Flight  time  of  the  previous  round. 

b.  Round  loading. 

c.  Relaying  of  the  weapon,  which  may  or  may  not  reflect  efforts 
to  improve  the  aimpoint  with  reference  to  the  aimpoint  used  for  the 
previous  round. 

2.33  Fixed  Time. 

Portion  of  total  firing  time,  for  a  first  or  a  subsequent  round, 
that  is  independent  of  environmental  conditions,  crew  skill  and  level  of 
training,  or  other  causes  of  variation  in  firing  time  performance.  This 
time  is  frequently  associated  with  mechanical  operations,  such  as  automatic 
loading,  for  unguided  rounds  or  with  flight  time  of  the  previous  round 
for  missiles. 

2.34  Variable  Time. 

Component  of  total  first  round  or  subsequent  round  firing  time 
that  needs  to  be  represented  by  a  distribution  of  times,  rather  than  by  a 
fixed  time  component.  On  any  particular  firing  occasion,  total  firing 
time  equals  the  sum  of  the  applicable  fixed  time  component  and  a  particular 
time  from  the  distribution  describing  the  variable  time  component. 

Variable  time  is  tied  to  human  operations  and  generally  tends. not  to 
equal  zero,  because  there  exist  so  many  causes  of  variation  that  are  not 
controllable,  at  least  completely. 

2.35  Median  Time. 

Median  time  can  be  defined  as  the  particular  time  which  is 
greater  than  half  the  times  represented  by  the  corresponding  distribution 
and  smaller  than  the  remaining  half.  For  the  special  case  where  a 
1  ogari thmi conormal  distribution  is  fitted  to  data  consisting  of  N  particular 
times,  the  median  time  is  the  Nth  root  of  the  product  of  all  times  in  the 
set  of  data  considered.  Alternatively,  for  this  special  case,  the  median 
time  is  the  anti  1 ogari thm  of  the  arithmetic  mean  of  the  N  logarithms 
corresponding  to  the  times  constituting  the  basic  data.  Detailed  studies 
made  many  years  ago  established  that  logarithmi conormal  distributions 
seemed  to  correspond  more  closely  than  other  known  distributions  to  test 
data  obtained  for  tank  weapons  of  interest  at  the  time.  Distributions  of 
this  type  are  still  considered  useful  to  represent  variable  components  of 
firing  times  for  weapons  of  current  concern. 
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2.36  Variability  Factor. 


Factor  indicating  to  what  extent  the  times  represented  by 
1 ogarithmiconormal  distributions  vary  with  reference  to  the  applicable 
median  times.  A  factor  of  about  0.5  (based  on  calculations  with  natural 
logarithms)  has  been  found  to  apply  consistently  and  is  normally  used, 
except  when  a  specific  reason  for  doing  otherwise  is  identified.  This 
factor  is  the  standard  deviation  of  the  natural  logarithms  of  time  data 
for  a  particular  set  of  conditions  with  reference  to  the  natural  logarithm 
of  the  median  time. 

2.37  Minimum  Time. 


Least  time  considered  possible.  Such  a  least  time,  if  one  other 
than  zero  is  specified,  can  be  made  to  override  any  unrealistically 
shorter  times  that  may  be  implied  by  a  particular  time  distribution  used. 

2.38  Kill  Probability  for  Specified  Engagement  Time. 

Probability  that  the  target  is  killed  when  a  specific  engagement 
time,  say  T  seconds,  has  elapsed  since  the  beginning  of  the  firing 
engagement.  Specification  of  a  kill  criterion  is  also  necessary.  T  is 
considered  to  equal  0  to  120  seconds  for  the  program  calculations.  The 
kill  probability  associated  with  any  specific  time  T  is  an  aggregated 
result  that  accounts  for  kills  that  occur  at  any  time  up  to  T. 
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3.  GENERAL  FLOW  OF  PROGRAM  LOGIC 


3.1  Introduction. 

This  section  provides  a  general  indication  of  the  principal 
steps  involved  in  the  simulation  of  firing  engagements.  Such  an  indication, 
which  may  suffice  for  many  readers,  is  supplemented  by  the  detailed 
explanations  of  the  program  logic  that  are  contained  in  the  appendixes. 

3.2  Principal  Program  Steps. 

The  program  instructions  cause  the  computer  to  perform  a  sequence 
of  steps.  The  flow  can  be  complicated,  principally  when  special  purpose 
options  are  exercised.  The  principal  steps  for  the  basic  program  (exclusive 
of  special  purpose  options)  are  as  follows: 

a.  Read  input  cards  for  the  weapon/round/fire  control 
type,  target  type,  target  exposure  conditions,  and  range  involved. 

b.  Make  all  necessary  adjustments  to  the  card  input  data  just 
read;  for  example,  some  data  may  require  conversion  from  meters  to  inches. 

c.  Skip  any  files  on  the  vulnerability  data  tape  that  precede 
the  kill  probability  and  personnel  casualty  data  needed. 

d.  Read  kill  probability  data  from  the  vul nerabil ity  data  tape 
for  target  orientation  angles  of  0  and  180,  30  and  210,  60  and  240,  and 
90  and  270  degrees.  The  two  angles  in  each  pair  correspond  to  reverse 
directions.  It  is  convenient  to  refer  to  the  first  angle  as  basic  and  to 
the  second  as  the  corresponding  reverse  angle.  Steps  e  through  j  apply, 
in  turn,  to  each  pair  of  basic  and  reverse  directions. 

e.  Simulate  the  first  firing  engagement  for  the  basic  orientation 
angle.  See  3.3  for  details. 

f.  Simulate  the  second  firing  engagement  for  the  basic  orientation 
angle,  as  for  e. 

g.  Likewise,  simulate  additional  engagements  until  a  total  of 
10,000  engagements  for  the  basic  orientation  angle  has  been  reached. 

h.  Repeat  steps  e  through  g  for  the  reverse  orientation  angle. 

i.  Combine  data  from  all  10,000  firing  engagements  for  each 
angle  of  the  matched  pair  of  concern  to  determine  the  hit  probability  and 
kill  probabilities  corresponding  to  1 ,  2,  3,  etc.,  14,  15  rounds,  the 
average,  or  expected,  numbers  of  rounds  needed  to  hit  or  to  kill  the 
target,  and  the  hit  and  kill  probabilities  corresponding  to  engagement 
times  of  0,  2,  4,  etc.,  118,  120  seconds. 

j.  Print  as  output  the  results  obtained  in  step  i. 
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3.3  Program  Steps  for  Single  Engagement  at  Particular 

Angle. 

Simulation  of  a  single  firing  engagement  for  a  particular 
orientation  angle  begins  with  consideration  of  the  first  main  round 
fired.  Steps  involved  are  as  follows: 

a.  Determine  the  number  of  seconds  between  the  beginning  of  the 
engagement  and  the  arrival  of  the  first  round  at  the  target  range. 

b.  Determine  if  the  first  round  flies  reliably.  If  not,  the 
target  cannot  be  hit  or  killed  and  calculations  for  the  next  round  begin 
immediately  without  any  credit  for  sensing  of  the  first  round  as  a  basis 
for  adjustment  of  fire. 

c.  If  the  first  round  has  a  reliable  flight  trajectory,  determine 
the  impact  point  in  the  plane  of  the  target  and  whether  or  not  the  target 
is  hit. 


d.  If  the  first  round  misses  the  target,  determine  whether  the 

round  is  sensed  and  begin  calculations  for  the  next  round.  The  calculations 

for  the  next  round  eventually  reflect  whether  or  not  adjustment  of  fire, 

to  take  advantage  of  information  obtained  by  sensing,  has  been  possible. 

e.  If  the  round  hits,  determine  whether  reliable  functioning  of 
the  fuze  occurs.  If  not,  the  target  cannot  be  killed  and  calculations 
for  the  next  round  begin  without  any  further  consideration  of  possible 
damaging  effects  on  the  target. 

f.  If  damaging  effects  are  to  be  assessed,  determine  whether 

each  of  the  various  types  of  kill,  i.e.,  M  kill,  F  kill,  H  or  F  kill,  and 
K  kill,  is  inflicted  on  the  target.  Continue  with  the  calculations  for 
the  second  round  whenever  there  is  at  least  one  kill  criterion  according 
to  which  the  target  has  not  been  killed  by  the  first  round.  If  the  target 

has  been  killed  according  to  all  criteria,  proceed  to  simulation  of  the 

second  firing  engagement. 

Unless  the  engagement  is  over,  simulate  firing  of  the  second  main  round 
as  follows: 


g.  Determine  the  number  of  seconds  between  firing  of  the  first 
round  and  firing  of  the  second  round.  Add  this  number  to  the  time  of 
arrival  of  the  first  round  at  the  target  range  to  obtain  the  time  between 
the  start  of  the  engagement  and  the  arrival  of  the  second  round  at  the 
target  range. 


h.  Determine  if  the  second  round  flies  reliably.  If  not, 
proceed  as  for  b  above,  i.e.,  begin  calculations  for  the  next  round. 
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i.  If  the  second  round  has  a  reliable  flight  trajectory, 
determine  the  impact  point  in  the  plane  of  the  target  and  whether  or  not 
the  target  is  hit,  as  for  c.  The  second  round  is  assumed  fired  with  the 
same  aimpoint  as  for  the  first  round  if  the  latter  is  not  sensed,  or  with 
a  new  aimpoint  based  on  sensing  information. 

j.  If  the  second  round  misses  the  target,  determine  whether  the 
round  is  sensed  and  begin  calculations  for  the  next  round,  as  for  d. 

k.  If  the  round  hits,  determine  whether  damaging  effects  are 
not  to  be  assessed  because  of  unreliable  functioning  of  the  fuze.  If  the 
fuze  does  not  function  reliably,  begin  calculations  for  the  next  round, 
as  for  e. 


1.  If  damaging  effects  on  the  target  are  to  be  assessed, 
determine  whether  each  of  the  various  types  of  kill  still  being  considered 
is  inflicted  on  the  target.  Continue  with  the  calculations  for  the  third 
round  whenever  there  is  at  least  one  kill  criterion  according  to  which 
the  target  has  not  yet  been  killed,  or  proceed  to  simulation  of  the  second 
firing  engagement. 

If  the  engagement  is  not  yet  over  after  two  rounds,  simulate  firing  of 
the  third  round  and  succeeding  rounds,  as  necessary,  but  not  beyond  the 
fifteenth  round.  Steps  g  through  1  also  apply  to  each  round  following 
the  second  round.  Times  at  which  rounds  arrive  at  the  target  range  are 
ignored  whenever  these  times  are  beyond  2  minutes  from  the  start  of  the 
engagement. 

3.4  Special  Options. 

Listing  of  the  various  special  options  available  is  not 
attempted  here.  Two  examples,  selected  somewhat  arbitrarily,  are  as 
fol lows : 


a.  A  vulnerability  data  tape  or  disc  is  not  absolutely  required. 
It  may  be  possible  to  substitute  a  set  of  input  cards  containing  detailed 
target  shape  information  identical  to  that  on  the  tape  or  disc.  Target 
shape  data  are  sufficient  if  hit  probabilities,  but  not  kill  probabilities, 
need  to  be  computed.  Although  it  would  be  possible  to  use  input  cards, 
rather  than  a  tape  or  disc,  for  detailed  vulnerability  data  also,  this 

is  impractical. 

b.  Calculations  can  be  made  for  simplified  engagements  each 
involving  the  firing  of  only  one  round.  In  particular,  single  round  hit 
probabilities  and  kill  probabilities  can  be  obtained  for  a  moving  target 
or  a  moving  firing  tank.  However,  this  does  not  imply  that  the  program 
can  actually  simulate  in  detail  the  movement  of  the  target  or  firing  tank. 
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4.  INPUT  DATA 


4.1  Introduction. 

The  direct  fire  program  can  make  calculations  for  several 
combinations  of  conditions  considered  in  sequence.  However,  it  is 
sufficient  for  this  report  to  emphasize  a  single  basic  cycle  definable  as 
all  the  calculations  for: 

a.  A  specified  firing  weapon,  round  type,  and  fire  control  system. 

b.  A  particular  target. 

c.  One  target  exposure  condition,  i.e. ,  the  target  is  fully 
exposed  or  in  defilade. 

d.  One  range. 

e.  One  of  the  orientation  angles  0,  30,  60,  or  90  degrees  and 
the  corresponding  angle  for  the  reverse  direction.  Matched  pairs  of 
angles  are  0  and  180,  30  and  210,  60  and  240,  and  90  and  270  degrees. 

The  input  data  required  for  a  single  basic  cycle  consist  of  both  a  set 
of  cards,  and  normally,  a  tape  or  disc  file.  The  cards  contain  control 
quantities  as  well  as  information  related  to  aimpoint,  reliability, 
delivery  accuracy,  and  rate  of  fire.  The  tape  or  disc  provides  target 
vulnerability  data. 


This  report  documents  the  input  data  used  for  a  sample  computer 
run  and  the  corresponding  output  results.  For  this  sample  run,  it  was 
more  important  to  illustrate  input  and  output  quantities  in  detail  than 
to  consider  weapon,  target,  and  other  combat  engagement  conditions  of 
real  concern.  The  inputs  developed  for  the  purpose  just  stated  are 
identified  in  this  section.  They  are  suitable  for  illustrating  input 
requirements  and  formats  in  detail  and,  while  possibly  similar  to  classified 
data  of  actual  interest,  are  not  applicable  to  a  specific  set  of  engagement 
conditions. 

4.2  Card  Inputs. 

Contents  of  the  input  cards  required  to  run  the  direct  fire 
program  are  listed  and  briefly  descrbed  in  Table  4.1.  The  role  of  some 
input  quantities  should  be  evident  from  the  table.  Several  other  quantities, 
however,  can  be  fully  understood  only  in  connection  with  detailed  explanations 
subsequently  provided.  Cards  numbered  1  through  5  contain  controls  needed 
for  each  computer  run  and  applicable  to  each  range  regardless  of  how  many 
ranges  are  grouped  for  the  run.  Integer  controls  are  on  the  first  three 
of  these  cards,  and  real  quantities  on  the  last  two  cards.  Cards  numbered 
R1  through  R8  are  needed  for  each  range.  Except  for  the  integer  controls 
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on  the  R1  card,  all  input  data  on  these  cards  are  real.  The  EXAMPLE 
column  of  Table  4.1  contains  the  data  used  for  the  run  that  yielded  the 
sample  outputs  presented  in  section  5.  Since  IMILS  is  set  to  0,  meters 
are  the  unit  associated  with  cards  R3  through  R6.  Blanks  on  cards  R4 
through  R8  correspond  to  input  quantities  associated  with  special  options 
that  were  not  exercised  in  the  sample  run. 

The  manner  in  which  input  data  are  actually  arranged  on  cards  is 
illustrated  in  Table  4.2.  This  table  contains  the  same  inputs  as  the 
previous  one.  The  first  portion  of  the  table  shows  how  the  basic  control 
cards  (identifiable  by  R15  A,  B,  C,  D,  E  entered  in  the  last  8  columns) 
are  followed  by  the  8  cards  corresponding  to  250  meters  (identifiable  by 
R15  2501,  2,  3,  4,  5,  6,  7,  8  in  these  same  columns).  The  250-meter  cards 
would  be  followed  if  necessary  by  similar  groups  of  8  cards  for  other 
ranges.  The  remainder  of  the  table  repeats  the  card  images,  with  headings 
added  to  facilitate  relating  each  quantity  to  the  explanatory  information 
in  Table  4.1 . 

4.3  Tape  or  Disc  Inputs. 

Vulnerability  data  constituting  needed  tape  or  disc  input 
information  are  generated  by  the  Vulnerability  Division  of  the  Ballistic 
Research  Laboratory  (BRL).  These  data  are  extensive  in  scope  and  consist 
of  kill  probabilities  for  a  large  number  of  specific  impact  points  on 
the  target.  These  impact  points  are  selected  from  the  grid  squares 
obtained  by  dividing  the  target  as  shown  in  Figure  4.1.  These  squares, 
usually  called  cells,  are  either  4  inches  (101.6mm)  or  100mm  on  each 
side.  One  impact  point  is  chosen  from  each  cell  and  data  developed  for 
that  point  are  then  considered  representative  of  the  effects  of  a  hit 
anywhere  on  the  cell.  Table  4.3  illustrates  a  portion  of  the  data.  The 
associated  conditions  are  not  identifiable,  since  fictitious  projectile, 
target,  and  range  headings  appear  in  the  table.  The  coordinates  X  and  Y 
are  those  of  the  center  of  each  cell  and  refer  to  a  suitably  selected 
origin,  e.g.,  at  the  turret  center  line  for  some  tank  targets  or  at  the 
base  of  the  target.  Each  line  of  data  includes  probabilities  of  M  kill, 

F  kill,  M  or  F  kill  (M/F  in  the  table),  and  K  kill,  as  well  as  expected 
casualties  when  applicable,  given  a  hit  on  the  cell  of  interest  for  both 
a  basic  orientation  angle  and  the  corresponding  reverse  angle.  Note 
that  the  reverse  angle  corresponding  here  to  30  degrees  is  150  degrees. 

This  reflects  the  assumption,  established  by  BRL,  that  a  single  set  of 
vulnerability  data  can  reasonably  be  used  for  symmetrically  related 
directions,  like  150  and  210  degrees,  that  represent  attack  of  a  target 
from  the  right  in  one  instance  and  from  the  left  in  the  other. 

The  vulnerability  data  used  for  the  sample  computer  run  differ 
from  those  partly  shown  in  Table  4.3,  but  need  not  be  identified  further. 

The  points  labelled  A,  B,  C,  and  D  in  Figure  4.1  are  the  corners 
of  a  rectangle  just  large  enough  to  enclose  the  target  representation. 

Such  a  rectangle  is  referred  to  in  some  of  the  explanations  in  Appendix  A. 
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CARD  R15  D  PASSN,WCELL,XC*YC>YBASEr YTOPrRELT 
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TABLE  4.2  SAMPLE  INPUT  CARD  DATA  (CONTINUED) 
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Note;  A,  B ,  C,  and  D  are  corners  of  rectangle  enclosing  target 
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flGURE  4.1  MEDIUM  TANK  TARGET  AT  ORIENTATION  OF  60  DEGREES 


TABLE  4.3  SAMPLE  DATA  FROM  VULNERABILITY  TAPE 
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666 


4.4  Requirements  for  Special  Options. 


A  special  option,  already  mentioned  in  3.4,  permits  using  input 
cards  with  target  shape  information  instead  of  the  vulnerability  data 
tape  or  disc  normally  required  for  a  run.  Card  requirements  associated 
with  this  option  are  illustrated  in  Table  4.4.  Numbers  are  grouped  by 
threes.  For  example,  76,  -28,  -8  mean  that  the  target  includes  an  unbroken 
string  of  cells  whose  centers  have  coordinates  (-28,76),  (-24,76),  (-20,76), 
(-16,76),  (-12,76),  and  (-8,76),  while  88,  -24,  -24  simply  identify  an 
isolated  cell  the  coordinates  of  whose  center  are  (-24,88).  The  total 
target  can  thus  be  specified  as  a  combination  of  horizontal  strips  of 
cells  and,  normally,  some  isolated  cells  toward  the  edges.  The  target 
shape  card  option  has  been  restricted  to  only  four  orientations,  namely, 

0,  30,  60,  and  90  degrees,  identified  in  columns  77  and  78  on  the  cards. 

The  sample  run  documented  in  this  report  did  not  involve  target  shape 
data  cards. 

The  direct  fire  program  has  been  designed  to  permit  certain 
calculations  for  a  moving  target  or  a  moving  firing  tank.  Single  round 
hit  probabilities  and  kill  probabilities  are  obtainable  for  sets  of  up  to 
18  conditions.  Input  card  requirements  are  then  as  explained  in  connection 
with  Table  4.1,  except  for  two  differences.  First,  on  the  R1  card  for 
only  the  first  range  involved  in  the  run,  INTPL  does  not  have  the  previously 
explained  meaning  but  should  be  set  to  987.  The  second  difference  is 
that  each  R3  card  is  replaced  by  a  set  of  18  such  cards.  It  is  assumed 
that,  from  other  sources,  biases  and  overall  standard  deviations  including 
contributions  of  lay  and  round-to-round  errors  are  known.  Thus,  these 
overall  standard  deviations  can  be  entered  in  SIGXB  and  SIGYB,  redefined 
accordingly;  redefinition  of  the  biases  XB,  YB  is  not  necessary.  Tape  or 
disc  input  requirements  are  not  affected.  Sections  of  the  computer 
program  that  specifically  apply  to  moving  target  or  moving  firing  tank 
calculations  were  not  exercised  in  the  sample  run  made  for  this  report. 
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TABLE  4. A  SAMPLE  CARDS  WITH  TARGET  SHAPE  DATA 
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TABLE  A, A  SAMPLE  CARDS  WITH  TARGET  SHAPE  DATA  (CONTINUED) 
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5.  OUTPUT  DATA 


5.1  Introduction. 

Output  results  provided  by  the  direct  fire  program  consist 
primarily  of  cumulative  hit  probabilities  and  kill  probabilities  for  up 
to  15  rounds,  average  numbers  of  rounds  expected  to  be  needed  to  hit  or 
kill  targets,  and  cumulative  hit  probabilities  and  kill  probabilities 
versus  engagement  time  for  up  to  2  minutes.  In  addition,  cumulative 
passenger  casualties  expected  for  up  to  15  rounds  are  included  for 
personnel  carrier  targets.  Information  useful  for  verifying  the  input 
data  used  also  constitutes  a  part  of  the  printed  output. 

5.2  Outputs  of  Sample  Computer  Run. 

The  input  quantities  identified  in  Table  4.1  (and  Table  4.2)  anc 
data  from  a  vulnerability  data  tape  were  used  to  make  a  sample  run. 

Table  5.1  presents  selected  results  of  this  run. 

All  data  from  the  input  cards  used  for  a  run  are  rewritten 
exactly  as  they  appear  on  the  cards.  The  same  data  are  generally  repeated 
with  corresponding  program  symbols  and  a  brief  indication  of  what  the 
symbols  mean.  Exceptionally,  quantities  associated  with  the  program 
symbols  may  not  be  strictly  identical  to  the  related  card  quantities. 

For  example,  the  dimension  of  target  cells  may  be  in  millimeters  originally 
and  converted  to  inches. 

Identifying  conditions  associated  with  input  information  provided 
by  the  vulnerability  data  tape  or  disc  are  printed,  but  the  actual  vulner¬ 
ability  data  are  not. 

The  number  of  samples  is  not  calculated  but  has  been  fixed  in 
the  program  itself  at  10,000. 

Angles  1  and  2  designate  target  orientation  directions  of  0  and 
180  degrees  respectively.  Data  already  identified  as  the  primary  output 
results  of  a  run  are  printed  for  each  such  related  pair  of  target 
orientations  (angles  1  and  2,  3  and  4,  etc.). 

Cumulative  hit  probabilities  and  kill  probabilities,  as  well  as 
personnel  casualties  when  applicable,  for  up  to  15  rounds  are  arranged  in 
twelve  columns  of  15  quantities  each.  Reading  horizontally,  one  sees  two 
sets  of  five  columns  followed  by  two  columns.  The  first  five  columns 
correspond  to  the  first  angle  (0  degrees  for  the  first  pair  of  angles) 
and  are,  from  left  to  right,  cumulative  probabilities  of  achieving  an  M 
kill,  an  F  kill,  an  M  or  F  kill,  a  K  kill,  or  a  hit  on  the  target.  The 
next  five  columns  contain  this  same  information,  similarly  arranged,  for 
the  second  angle  (180  degrees  for  the  first  pair).  The  last  two  columns, 
printed  only  when  they  apply,  are  cumulative  passenger  casualties  for  the 
two  angles  (0  and  180  degrees)  respectively. 
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TABLE  5,1  SAMPLE  OUTPUT  DATA 


R159851340 


3  10 

250  2  0  0  15  0 

1 

0  0 

0 

0 

0 

0  0 

0  0  0  0  0  0 

NRSKP 

VULN.  TAPE  FILES  TO  SKIP 

3 

NCASES 

NO.  OF  CASES 

10 

NRFRST 

FIRST  RANGE  FOR  VULN.  DATA 

250 

NTGT 

VULN.  TAPE  FORMAT 

2 

NRDTYP 

VULN.  DATA  DEP.  ON  RANGE 

0 

NRD1 

SINGLE  ROUND  ENGAGEMENT 

0 

NRDS 

MAX.  ND.  OF  RDS. 

15 

NADJST 

SP.  PRDC.  FOR  MISS.  RD. 

0 

NRFBS 

RANGEFINDER  DR  BATTLESIGHT 

1 

NDROP 

STANDARD  DROP  ADJUSTMENT 

0 

NHIT 

ADJ.  HIT 

0 

NPRHIT 

HIT  PRDB.  DUTPUT  ONLY 

0 

NDTRM1 

DETERMINISTIC2  HIT1 

0 

NDTRM2 

DET2  ND  K  KILL  HIT1 

0 

NDTRM3 

DET2  MISS1 

0 

NDTRM4 

DET2  OBSERVED  MISS1 

0 

NDTRM5 

DET2  UNOBSERVED  MISS1 

0 

NHTKLL 

HIT  AND  KILL  PROB.  INPUT 

0 

NTCRD1 

CARDS  FDR  TARGET1 

0 

NTCRD2 

CARDS  FOR  TARGET2 

0 

NTCRD3 

CARDS  FOR  TARGET3 

0 

NTCRD4 

CARDS  FOR  TARGET4 

0 

6.0000 

4.0000  0.0000  -28.0000 

-70.0000  9999.9999 

.9900 

.9280 

HOP,. 

VERT. 

PASSN 

TGT.  PERS.  EXCL.  CREW 

6.0000 

W,W 

CELL  DIM.  IN  INCHES 

4.0000 

4.0000 

XC,  YC 

TGT.  CTR.  CDORD.  INS. 

0.0000 

-28.0000 

YBASE 

TGT.  BDT •  CDDRD*  INS. 

-70.0000 

YTOP 

TOP  APPR.  CODRD.  INS. 

9999.9999 

RELT 

PROB.  OF  REL.  FLIGHT 

.9900 

RELF 

PROB.  OF  REL.  FUZE 

.9280 
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TABLE  5.1  SAMPLE  OUTPUT  DATA  (CONTINUED) 


250  0  0 


1.7100 

0.0000 

0.0000 

19.8000 

15.6000 

0.0000 

0.0000 

0.0000 

0.0000 

.0357 

1.0161 

2.7038 

5.5923 

.2277 

.2277 

.6136 

.6872 

.7800 

1.4635 

1.4635 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0. 0000 

TIME  DATA 

(SECONDS) 

1ST  P.D,  SUBS.  RDS . 

FLT 

TIME  OF  FLIGHT 

1. 7100 

TF1,TFS 

FIXED  TIME 

0.0000 

0.0000 

XM1»  XM2 

MEDIAN  VARIABLE  TIME 

19.8000 

15.6000 

AMT1*  AMT2 

MINIMUM  TIME 

0.0000 

0.  0000 

STDl*  STD2 

VARIABILITY  FACTOR 

0. 0000 

0.0000 

ACCURACY  DATA  (METERS 

HOR. 

VERT. 

XB*  YB 

FIXED  BIAS 

.0357 

1.0161 

SIGXB*  SIGYB 

VAR.  BIAS  STD.  DEV. 

2.7038 

5.5923 

SIGXL/SIGYL 

LAY  ERROR  STD.  DEV. 

.2277 

.2277 

SIGXR,SIGYR 

RD.-TO-RD.STD.  DEV. 

.6136 

.6872 

PROBS 

PR OB.  SENS.  MISS.  RD. 

.  7800 

SIGXSjSIGYS 

SENS.  ERROR  STD.  DEV. 

1.4635 

1.4635 

PROJECTILE  985  VS.  VEHICLE  13*0  087 


TASK  SHEET  IDENTIFICATION  R15 
PROJECTILE  CODE  985 
TARGET  CODE  1340 
RANGE  IN  METERS  250 


X  Y 


0  DEG  ELEV*  0  DEG  AZIM 
M  F  M/F  K  EC 


0  DEG  ELE V#  180  DEG  AZIM 
M  F  M/F  K  EC 


X  Y 


0  DEG  ELEV,  30  DEG  AZIM 
M  F  M/F  K  EC 


0  DEG  ELEV,  150  DEG  AZIM 
M  F  M/F  K  EC 


X  Y 


0  DEG  EL EV»  60  DEG  AZIM  0  DEG  ELEV,  120  DEG  AZIM 

M  F  M/F  K  EC  M  F  M/F  K.  EC 
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X 

1BE 

5LE 

RDS 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

1 

2 

3 

4 

5 

6 

7 

8 

9 


TABLE  5.1  SAMPLE  OUTPUT  DATA  (CONTINUED) 


0  DEG  ELEV#  90  DEG  AZIM 


0  DEG  ELEV#  270  DEG  AZIM 


Y  M  F  M/F  K  EC 


M  F  M/F  K  EC 


OF  SAMPLES  *  10000 

-  1  2 


PRDB.  (EXP 


.0310000 

.1578000 

.2747000 

.3877000 

.4793000 

.5568000 

.6249000 

.6790000 

.7284000 

.7667000 

.8001000 

.8315000 

.85?8000 

.8787000 

.8978000 

.0299000 
.1499000 
.2701000 
.3‘,03000 
.4575000 
.5344000 
.5982000 
.6579000 
.7050000 
.7473000 
.7826000 
i 8 143 000 
.8402000 
.8646000 
.8826000 

.0828622 

.4197333 

.7544151 

1.0759812 

1.3500448 

1.5879577 


CAS.) 

. 0300000 

.0310000 

.0212000 

.0352000 

.1529000 

.1585000 

.1060000 

.1762000 

.2677000 

.2754000 

.1969000 

.3043000 

.3763000 

.3882000 

.2884000 

.4238000 

.4681000 

.4796000 

.3695000 

.5178000 

.5450000 

.5571000 

.4407000 

.5951000 

.6131000 

.6255000 

.5057000 

.6646000 

.6683000 

.67O8000 

.5600000 

.7178000 

.7178000 

.7292000 

.6114000 

.7655000 

.7566000 

.7674000 

.6515000 

.8028000 

.7900000 

.8008000 

.6888000 

.8348000 

.8217000 

. 8320000 

.7246000 

.8625000 

.8494000 

.8582000 

.7598000 

.8859000 

.8715000 

.8792000 

.7879000 

.9036000 

.8916000 

.8983000 

.8124000 

.9199000 

.0289000 

.0299000 

.0215000 

.0368000 

.1429000 

.1499000 

.1057000 

.1762000 

.2601000 

.2701000 

.1072000 

.3101000 

.3581000 

.3703000 

.2785000 

.4222000 

.4447000 

.4576000 

.3508000 

.5145000 

.5194000 

.5348000 

.4215000 

.5913000 

.5838000 

.5987000 

.4826000 

.6565000 

.6431000 

,6582000 

.5387000 

.7137000 

.6903000 

.7053000 

.5864000 

.7576000 

.7331000 

.7474000 

.6313000 

.7979000 

.7694000 

.7827000 

.6691000 

.8312000 

.8028000 

.8144000 

.7039000 

.8602000 

.8287000 

.8403000 

.7356000 

.8826000 

.8539000 

.8647000 

.7673000 

.9027000 

.8716000 

.8827000 

.7903000 

.9191000 

.0928149 

.4444624 

.8158528 

1.1396685 

1.4204448 

1.6770897 

1.8020048 

1.9762646 

2.1348499 

2.2621322 

2.3776576 

2.4867633 

2.5822915 


1.8936590 

2.0894381 

2.2512219 

2.3986‘>33 

2.5282113 

2.6408755 

2.7402884 
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TABLE  5,1  SAMPLE  OUTPUT  DATA  (CONTINUED) 


14  2.6566855  2.8284182 

15  2.7242920  2.8998553 


AVG.  NO.  OF  RDS • 


1 

2 

3 

4 

5 

6 

7 

8 
o 

10 


7.62 

7.80 

7.61 

o.91 

7.00 

8.00 

8.28 

8.00 

10.48 

-\03 


PROBABILITY  VERSUS  TIME 


0.00000 

0.00000 

0.00000 

0.00000 

.00050 

.00150 

.00250 

.00630 

•00880 

.012  80 

.01730 

.02320 

.02820 

.03500 

.04290 

.05290 

.06290 

.07420 

.08430 

. 0965C 

.10700 

,11?30 

.12940 

.14250 

.15470 

. 16°40 

.18300 

.19420 

.20820 

.2202C 

.23260 

.24510 

.25840 

.27030 

.28290 

.29720 

.31110 

.32260 

.33390 

•  3  44  5C 

.35540 

.36590 

•  3  7790 

.39030 

.40020 

.41080 

.42180 

.43100 

.44080 

•  4  51 1C 

.46000 

.55140 

.47050 

.47990 

.48920 

,40040 

.50900 

.51750 

.52550 

.*3400 

• 5426C 

0.00000 

0.00000 

0.00000 

0.00000 

.00050 

.00150 

.00240 

.00620 

.00870 

• 0126C 

.01700 

.02270 

.02760 

.03440 

.04190 

.05180 

.06140 

.07240 

.08220 

•  094 1C 

.10430 

.11390 

.12540 

.13810 

.15010 

.16460 

.17740 

.18860 

.20230 

•  21 36C 

.22600 

.23820 

.25090 

.26310 

.27550 

.28060 

.30320 

.31410 

.32540 

•  335  7C 

.34600 

.35670 

.36820 

.38050 

.39000 

.40030 

.41070 

.41980 

.42950 

•  4  39  3C 

.44800 

.54050 

.45830 

.46750 

.47680 

• 486°0 

.49660 

.50520 

.51390 

.52290 

• 5316C 

0.00000 

0.00000 

0.00000 

0.00000 

.00050 

.00150 

.00250 

.00630 

.00880 

•  012  8C 

.01730 

.02320 

.02820 

.03500 

.04310 

.05310 

.06310 

.07440 

.08450 

•  096  8C 

.10750 

.11770 

. 12090 

.14300 

.15530 

.17010 

.18370 

.19480 

.20880 

•  2  207C 

.23310 

.24560 

.25800 

.27080 

.28340 

.29770 

.31150 

.32290 

.33420 

.  3447C 

.35560 

.36610 

.37810 

.39050 

.40040 

.41100 

.42200 

.43140 

.44120 

•  451 5C 

.46040 

.55220 

.47090 

.48030 

.4  8960 

.49980 

. 50°20 

.51790 

.52600 

.53450 

• 5432C 

0.00000 

0.00000 

0.00000 

0.00000 

.00030 

.00090 

.00170 

.00430 

.00610 

• 0086C 

.01180 

.01560 

.01960 

.02410 

.02920 

.03550 

.04240 

.05050 

.05750 

•  0662C 

.07340 

.08030 

.08780 

,09690 

.10670 

.11740 

.12650 

.13500 

.14570 

• 1545C 

.16410 

.17360 

.18510 

.19490 

.20570 

.21700 

.22820 

.23880 

.24730 

•  2  56  5C 

.26650 

.27550 

.28570 

.29670 

.30550 

.313^0 

.32300 

.33110 

.34040 

•  3  49  4C 

.35710 

.44050 

.36700 

.37560 

.38270 

.39210 

.40130 

. 40°00 

.41760 

.42660 

•  4  334C 

0.00000 

0.00000 

0.00000 

0.00000 

.00050 

.00160 

.002_’0 

.00700 

.01000 

•  0 146C 

. 01930 

.02560 

.03150 

.03870 

.04800 

.05930 

.07070 

.08330 

.09540 

.  1 090C 

.12150 

.13240 

.14530 

.15980 

.17320 

.18990 

.20440 

.21680 

.23120 

. 244 1C 

.25710 

.27050 

.28400 

.29790 

.31110 

.32680 

.34120 

.35340 

.36570 

• 3768C 

.38800 

.39900 

.41220 

.42480 

.43500 

.44570 

.45780 

.46810 

.47780 

• 4884C 

.49720 

.58990 

.50800 

.51770 

.52680 

.53650 

.54590 

.55510 

.56320 

.57210 

• 5810C 

0.00000 

0.00000 

0.00000 

.00010 

.00030 

.00110 

.00220 

.00470 

.00740 

. 0098C 
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TABLE  5,1  SAMPLE  OUTPUT  DATA  (CONTINUED) 


.01390 

.01920 

.02660 

.03320 

.04170 

.05010 

.06110 

.07220 

.08250 

.09640 

.10670 

.11650 

.13000 

.14090 

.15260 

.16360 

.17480 

.18730 

.20130 

.21280 

.22430 

.23760 

.248^0 

.26090 

.27360 

.28650 

.29900 

.31060 

.32160 

.33400 

.34450 

.35500 

.36580 

.37610 

.38480 

.39500 

,40670 

.41690 

.42640 

.43560 

.144420 

.53010 

.45290 

.46170 

.47050 

.47990 

.48Q00 

.49690 

.50610 

.51360 

.52190 

0.00000 

0.00000 

0.00000 

.00010 

.00030 

.00110 

.00220 

.00440 

.00710 

.00940 

.01340 

.01840 

.02560 

.03170 

.03980 

.04790 

.05860 

.06940 

.07950 

.09260 

.10220 

.11160 

.12440 

.13460 

.14580 

.15670 

.16780 

.17950 

.19290 

.20470 

.21530 

.22850 

.23940 

.25110 

.26340 

.27590 

.28720 

.29820 

.30910 

.32140 

.33200 

.34240 

.35340 

.36360 

.37240 

.38310 

.39430 

.40450 

.41370 

.42310 

.43200 

.51750 

.44100 

.44990 

.45770 

.46710 

.47610 

.48410 

.49310 

.50050 

.50870 

0.00000 

0.00000 

0.00000 

.00010 

.00030 

.00110 

.00220 

.00470 

.00740 

.00980 

.01390 

.01920 

•C2660 

.03320 

.04170 

.05010 

.06110 

.07220 

.08250 

.09640 

.10670 

.11650 

.13000 

.14090 

.15260 

.16360 

.17480 

.18730 

.20130 

.21280 

.22430 

.23760 

.24870 

.26090 

.27360 

.28650 

.29900 

.31060 

.32160 

.33400 

.34460 

.35510 

.36590 

.37620 

.38490 

.39510 

.40680 

.41700 

.42650 

.43570 

.44440 

.53040 

.45310 

.46190 

.47070 

.48020 

.48930 

.49720 

.50640 

.51390 

.52220 

0.00000 

0.00000 

0.00000 

.00010 

.00020 

.00100 

.00200 

.00370 

.00590 

.00720 

.01000 

.01400 

.01960 

.02450 

.03070 

.03620 

.04320 

.05160 

.05970 

.06890 

.07660 

.08350 

.09290 

.10180 

.10900 

.11710 

.12620 

.13570 

.14570 

.155  30 

.16360 

.17590 

.18390 

.19300 

.20200 

.21110 

.22050 

.22960 

.23750 

.24740 

.25710 

.26660 

.27650 

.28510 

.29190 

.30120 

.31110 

.32020 

.32820 

.33560 

.34400 

.42200 

.35230 

.36040 

.36780 

.37470 

.38230 

.39100 

.39910 

.40670 

.41440 

0.00000 

0.00000 

0.00000 

.00010 

.00040 

.00130 

.00280 

.00560 

.00920 

.01250 

.01730 

.02430 

.03270 

.04040 

.05000 

.05980 

.07240 

.08530 

.09720 

.11330 

.12530 

.13680 

.15200 

.16510 

.1^790 

.19050 

.20340 

.21680 

.23260 

.24540 

.25920 

.27330 

.28620 

.29970 

.31360 

.32670 

.34030 

.35250 

.36460 

.37870 

.39020 

.40150 

.41400 

.42510 

.43420 

.44450 

.45780 

.46840 

.47820 

.48760 

.49690 

.58580 

.50600 

.51550 

.52440 

.53430 

.54370 

.55060 

.56050 

.56940 

.57740 

ETC.  .  .  .  FOR  OTHER  ANCLES 
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Average  numbers  of  rounds  expected  to  be  needed  to  hit  or  kill 
targets  consist  of  two  sets  of  five  numbers  each.  The  two  sets  correspond 
to  the  two  target  orientations  respectively.  The  ordering  within  each 
set,  according  to  the  type  of  kill  or  whether  only  a  target  hit  is 
considered,  is  the  same  as  for  the  ten  groups  of  hit  probabilities  and 
kill  probabilities  discussed  in  the  preceding  paragraph. 

Ten  groups  of  cumulative  hit  probabilities  and  kill  probabilities 
versus  engagement  time  follow.  Five  groups  apply,  as  before,  to  each 
target  orientation.  The  correspondence  with  kill  type  or  whether  only  a 
target  hit  is  of  concern  also  remains  as  before.  Each  group  contains  61 
numbers.  These  are  cumulative  probabilities  for  engagement  times  of  0, 

2,  4,  etc.,  118,  120  seconds. 

It  is  recognized  that  the  formats  used  to  print  output  proba¬ 
bilities  and  expected  personnel  casualties  contain  several  unnecessary 
decimal  places.  However,  all  results  have  been  consistently  rounded  to 
two  decimal  places  for  the  tank  effectiveness  pamphlets  published  by 
the  JMEM-SS  Methodology  and  Evaluations  Working  Group. 

The  choice  of  10,000  samples  that  was  made  results  in  the 
computed  output  quantities,  mostly  probabilities,  being  practically 
always  accurate  to  within  0.005.  This  degree  of  accuracy  is  basically 
equivalent  to  that  achievable  by  deterministic  programs. 
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6.  COMPUTER  ASPECTS 


The  direct  fire  program  is  coded  in  FORTRAN  IV. 


Production  runs  have  been  made,  with  the  program,  on  the  CDC 
7600  computer  at  Aberdeen  Proving  Ground,  Maryland  and  the  CDC  6600 
computer  at  Oklahoma  State  University  Field  Office,  Eglin  Air  Force  Base, 
Florida.  Also,  it  has  recently  been  learned  that  the  Oklahoma  State 
University  Field  Office  will  use  a  newer  CYBER  176  computer  for  future 
production  runs. 

Running  times  can  vary  according  to  various  factors,  e.g.,  the 
computer  used  and  the  program  options  exercised.  A  basic  run,  i.e. ,  one 
primarily  intended  to  simulate  10,000  engagements  for  several  target 
ranges  and  four  pairs  of  target  orientation  angles  at  each  range,  would 
typically  require  about  1  minute  of  computer  time  for  each  range  on 
either  the  CDC  7600  or  CYBER  176. 
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APPENDIX  A 


PROGRAM  LISTING  AND  DETAILED  EXPLANATIONS 
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PROGRAM  LISTING  AND  DETAILED  EXPLANATIONS 


The  principal  objective  of  this  appendix  is  to  provide  a  listing 
of  FORTRAN  statements  and  detailed  explanations  for  the  main  program. 

For  completeness,  however,  subroutine  NRAN  and  function  CNORM  are  also 
listed  and  briefly  described. 


Both  the  listing  of  program  statements  and  the  corresponding 
detailed  explanations  cover  several  pages.  Listing  and  explanation  pages 
have  been  intermixed  so  that  program  statements  and  accompanying  explanations 
are  reasonably  near  together. 


Some  input  quantities  have  already  been  sufficiently  explained 
in  Section  4.  Additional  explanations  related  to  inputs  are  included 
only  as  necessary  to  supplement  that  section. 


Main  program  lines  are  numbered  from  1  through  1385.  These 
lines  were  grouped  as  seemed  most  convenient  to  facilitate  the  explanations. 


Certain  program  statements  require  no  explanation, 
statements,  used  more  than  once,  are  covered  only  when  they 


Other 

first  occur. 


The  listing  and  explanations  for  the  main  program,  subroutine 
NRAN,  and  function  CNORM  follow. 
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PROGRAM  LISTING 


1 

PROGRAM  MEWGD  (  INPUT,  OUTPUT,  TAPE5 

- - - 

-INPUT,  TAPE6- 

OUTPUT,  TAPE2  ) 

'  COMMON  A( 80 ) 

*AFAIL (15) 

* AKIL (61*10) 

, AX ( 201 ) 

A  *  AY ( 51 ) 

*C(10) 

, IDCODEt 10) 

*  IK ( 50*2  00 ) 

5 

B  *  IP (10  > 

*ITAPE (80) 

, MPK  C  3000*8 ) 

*  NXTEMP ( 4 ) 

C  , NYTEMP (4  ) 

*  PK( 10 ) 

*  PKRLFT( 10 ) 

*  XTEMPL ( 4 ) 

D  *XTEMPH(4 ) 

* IMODXY { 3000*4) 

»YTEMPL  C  4) 

*  ZC 19, 12) 

E  *  XLEFTF (4  ) 

*XRGHTF(4) 

,XLEFTD(4) 

*  XRGHTD( 4) 

COMMON  KDMSTC ( 5 ) 

* B  X C  5) 

*  BY  (5  ) 

*  SI GX( 3, 5 ) 

10 

A  *  SIGY( 3*  5 ) 

*NMINUS ( 50*  2*5) 

»NPLUS  C  50*  2*  5 ) 

*COORD ( 2 ) 

B  *  TOTAL (2) 

*NPOS( 50) 

, NNEG( 50) 

»  NCMPOS ( 50 1 

C  *NCMNEG(50) 

*NRFLC  T(50) 

,BIASXY(2,5) 

, S IGMX Y( 3,2*5) 

D  *  PH<3»  5) 

,  PKILL  (4*  3,  5  ) 

COMMON  BBX( 4) 

*BB Y( 4 ) 

*  SS IGX( 4 ) 

»SS IGY { 4 ) 

15 

A  *  PPH(  5 ) 

*  PKSHOT (4*5) 

*  PPKHIT ( 4*  5) 

COMMON  KTRGTC (4 ) 

*  NCTRY (120*  4 ) 

,NCTRXL ( 120*4 ) 

»  NCTRXR( 120, 4 ) 

COMMON  XXB(3*6)*YYB(3*6)*SSIGXB(3*6) 

, SSIGYB (3,6) 

DATA  MASK11  / 

3777  B  / 

DATA  MASK16  / 

177777  B  / 

20 

WRITE  {  6,1010  ) 

1010  FORMAT  (  /  /  ) 

READ  (  5,1110  ) 

I  DC  ODE 

1110  FORMAT  (  10A1  ) 

WRITE  C  6*1112  ) 

I  DC  ODE 

25 

1112  FORMAT  (  1H1*  10X 

* 10A1  ) 

'  WRITE  (  6*1010  ) 

READ  (  5*1120 


30 


35 


A 

B 

C 


)  NRSKP  ,NCASES,NRFRST,NTGT 
*NADJST,NRFB$  ,NDROP  , NHIT 


*  NRDTYP, NRD1  »NRDS 
, NPRHIT*  NDTRM1*  ND”RM2 


*NDTRM3,  NDTRM4»NDTRM5,NHTKLL,NTCRD1*NTCRD2»  NTC  RD3 
*  NTCRD4 


1120 


FORMAT 
WRITE  C 


(  1415  ) 
6*1122 


A 

B 

C 


1122 


40 


45 


50 


55 


)  NRSKP  *NCASES*NRFRST*NTGT  *  NRDTYP*  NRD1  *  NRD S 
*NADJST*NRFBS  *NDROP  »NHIT  , NPRHIT, NDTRM1, ND'RM2 
*  NDTRM3* NDTRM4*  NDTRM5*  NHTKLL* NTCRD1* NTCRD2*  NTCRD3 
*NTCRD4 

FORMAT  C  10X* 1415  ) 

WRITE  (  6*1020  ) 

1020  FORMAT  (  ) 

IF  (  NRSKP  .EQ.  9999  )  NRSKP  *  0 
WRITE  C  6*1130  )  NRSKP  *NCASES*NRFRST»NTGT 
A  *NADJST*NRFBS  *NDROP  *  NHIT 

WRITE  C  6*1131  )  NDTRM3*  NDTRM4*  NDTRM5*  NHTKLL*  NTCRDl, NTCRD 2* NTC  RD3 


»NRDTYP*  NRD1  *  NRDS 
*NPRHIT*  ND  TRM 1*  NDTRM2 


1130  FORMAT 
A 
B 
C 
D 
E 
F 
G 
H 

J 

K 

L 

M 

N 


*  NTCRD4 

(  10X*  39H  NRSKP 
10X* 27H  NCASES 
10X*  41H  NRFRST 
10X*  32H  NTGT 
10X* 39H  NRDTYP 
10X*  3BH  NRD1 
10X*  31H  NRDS 
10X*  38H  NADJST 
10X* 41H  NRFBS 
10X*  39H  NDROP 
10X* 23H  NHIT 
10X*  36H  NPRHIT 


10X*  34H 
10X*  34H 


NDTRM1 

NDTRM2 


VULN.  TAPE  FILES  TO  SKIP*9X,I5  / 
NO.  OF  CASES*21X*I5  / 

FIRST  RANGE  FOR  VULN.  DATA, 7X, 15 
VULN.  TAPE  FORMAT* 16X* 1 5  / 

VULN.  DATA  DEP.  ON  RANGE*9X* 15  / 

SINGLE  ROUND  ENGAGEMENT*10X* 15  / 

MAX.  NO.  OF  RDS. * 17X* 15  / 

SP.  PROC.  FOR  MISS.  RD.*10X*I5  / 

RANGEFINDER  OR  B ATTLES I GHT, 7X,  1  5 
STANDARD  DROP  ADJUSTMENT*9X* 15  / 

ADJ.  HIT* 25X* 15  / 

HIT  PROB.  OUTPUT  0NLY*12X, 15  / 

DETERMINISTIC2  HIT1*14X*I5  / 

DET2  NO  K  KILL  HIT1*14X*I5  ) 
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EXPLANATIONS 


REFERENCES  TO  THE  LEFT 
OF  EXPLANATIONS  ARE  TO 
PROGRAM  LINE  NUMBERS 


1  TAPE  5  IS  ASSOCIATED  WITH  INPUT  CARDS  AND  TAPE  2  WITH  TAPE  OR  DISC  CONTAINING 

TARGET  VULNERABILITY  INPUT  DATA. 


3- 

17 

18— 

19 

ro 

o 

1 

165 

20- 

21 

22- 

25 

2’’- 

36 

37- 

38 

THESE  ARE  REALLY  DIMENSION  STATEMENTS.  OKLAHOMA  STATE  UNIVERSITY  FIELD  OFFICE 
AT  EGLIN  AIR  FORCE  BASE  HAS  FOUND  THAT  SYSTEMATIC  USE  OF  COMMON  STATEMENTS 
ENHANCES  COMPUTER  EFFICIENCY. 

THESE  DATA  STATEMENTS  ARE  RELATED  TO  RETRIEVAL  OF  KILL  PROBABILITY  AND 
PERSONNEL  CASUALTY  INFORMATION  FROM  FILES  CONTAINING  SUCH  INFORMATION  IN 
PACKED  FORM. 

PORTION  OF  PROGRAM  PRECEDING  PROCESSING  FOR  ANY  PARTICULAR  RANGE. 

SKIP  3  LINES. 

READ  AND  WRITE  INPUT  CARD  1. 

READ  AND  WRITE  INPUT  CARDS  2  AND  3. 

SKIP  1  LINE. 


39  NRSKP  =  9999  HAS  BEEN  TEMPORARILY  PUNCHED  ON  INPUT  CARDS  AS  CONSPICUOUS  ALERT 

THAT  PROPER  ENTRY  NEEDS  TO  BE  DETERMINED  LATER.  THIS  IS  SIMPLY  MATTER  OF 
CONVENIENCE.  TEMPORARY  VALUE  9^99  CAN  BE  RETAINED  AS  INPUT  IF  0  PROVES  TO  BE 
CORRECT  ENTRY#  BECAUSE  OF  RESETTING  DONE  HERE. 


40-  65  REWRITE  INPUTS  ON  CARDS  2  AND  3  WITH  ABBREVIATED  INDICATION  OF  WHAT  EACH 

NUMBER  REPRESENTS. 

66-  69  IF  NRD1  =  1#  NRDS  IS  IGNORED  AMD  NOT  USED.  OTHERWISE#  MAXIMUM  NUMBER  OF  ROUNDS 

TO  BE  FIRED  PER  ENGAGEMENT#  READ  AS  INPUT#  MAY  BE  OVERRIDDEN  EITHER  BY  LOWER 
OR  UPPER  BOUND.  REASON  FOR  LOWER  BOUND  OF  10  ROUNDS  IS  INDICATED  IN 
CONNECTION  WITH  PROGRAM  LINES  1333  THROUGH  1340.  UPPER  BOUND  OF  19  ROUNDS 
HAS  BEEN  ARBITRARILY  SELECTED  SIMPLY  TO  LIMIT  SCOPE  OF  PROGRAM  CALCULATIONS. 
RDS  IS  REAL  FORM. 


71-  77  IF  AT  LEAST  ONE  OF  CONTROLS  MDTRM1  THROUGH  NDTRM5  IS  NOT  0#  THEIR  SUM  NSMDTR 

IS  ALSO  NOT  0  AND  ALL  FIVE  QUANTITIES  ARE  ENTERED  IN  ARRAY  KDMSTC . 


78-  85  IF  AT  LEAST  ONE  OF  CONTROLS  MTCRD1  THROUGH  NT CP.D 4  IS  NOT  0#  THEIP.  SUM  NSTCRD 

IS  ALSO  NOT  0#  ALL  FOUR  QUANTITIES  ARE  ENTERED  IN  ARRAY  KTRGTC#  AND  PROGRAM 
LINES  86  THROUGH  95  ARE  SKIPPED  BECAUSE  TAPE  OR  DISC  WITH  VULNERABILITY  DATA 
IS  NOT  INVOLVED. 

87-  Q5  SKIP  FILES  AS  NECESSARY  ON  VULNERABILITY  DATA  TAPE  OR  DISC.  NSKP  FILES  SKIPPED 
ARE  NRSK°  FILES  MOT  APPLYING  TO  P.OUND/TARGET  COMBINATION  OF  CONCERN  AND 
POSSIBLY  ADDITIONAL  FILES  ASSOCIATED  WITH  RANGES  LESS  THAN  FIRST  ONE  TO  BE 
CONSIDERED.  STATEMENT  INVOLVING  GOING  TO  311  IF  E0F(IU)«1.0  PREVENTS  MINOR 
INCONSISTENCIES  IN  FILE  MARK  LOCATION  THAT  ARE  OCCASIONALLY  ENCOUNTERED  FROM 
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PROGRAM  LISTING  (CONTINUED) 


60 

65 

1131  1 
A 
B 
C 
D 
E 
F 
G 

70 

1132 

75 

1135 

80 

85 

11*5 

90 

95 

311 

rzoo 

100 

i 

1210 

i 

1212 

IF 

IF 

IF 

RDS 


IT  (  10X*25H 
10X*34H 
10X*36H 
10X/39H 
10X*32H 
10X*32H 
10X*32H 
10X*32H 
NRDl  .EQ.  1 
NRDS  .LT.  10 
NRDS  .GT.  19 
-  NRDS 


NDTRM3 
NDTRM4 
NDTRM5 
NHTKLL 
NTCRD1 
NTCRD2 
NTCRD3 
NTCRD4 
)  GO  TO 
)  NRDS 
)  NRDS 


1132 
•  10 
-  19 


DET2  MISS1*  23X* I  5  / 

DET2  OBSERVED  MISS1,14X,I5  / 
DET2  UNOBSERVED  MISS1*12X,  15  / 

HIT  AND  KILL  PROB.  INPUT*9X*I5 
CARDS  FOR  TARGET1*16X*I 5  / 

CARDS  FOR  TARGET2* 16X* 1 5  ! 

CARDS  FOR  TARGET3*16X* I  5  I 
CARDS  FOR  TARGET4* 16X* 15  ) 


CONTINUE 
NSMDTR  -  NDTRM1 
IF  (  NSMDTR  .EQ 
KDMSTC (1) 

KDMSTC (2) 

KDMSTC ( 3 ) 

KDMSTC ( 4 ) 

KDMSTC ( 5) 
CONTINUE 
NSTCRD  - 


+  NDTRM2  + 
0  )  GO  TO 

-  NDTRM1 

-  NDTRM2 

-  NDTRM3 
*  NDTRM4 

-  NDTRM5 

NTCRD1  +  NTCRD2  + 


NDTPM3 

1135 


+  NDTRM4  +  NDTRM5 


IF  (  NSTCRD  .EQ.  0 
KTRGTC ( 1)  -  NTCRD1 
KTRGTC ( 2 )  -  NTCRD2 
KTRGTC (3)  -  NTCRD3 
KTRGTC (4)  -  NTCRD4 
GO  TO  1200 
CONTINUE 
REWIND  2 
IU  -  2 
NS KP  - 
IF  ( 

IF  ( 

CALL 
READ 
IF  ( 

CONTINUE 
CONTINUE 
WRITE  (  6*1010 
READ  (  5*1210 


)  GO  TO 


NTCRD3 

1145 


♦  NTCRD4 


NRSKP 

NRDTYP  .EQ.  1  ) 
NSKP  .EQ.  0  )  GO 
SKIPFE  ( IU*NSKP*0) 
(  IU*  3110  )  A 
EOF(IU)  .EQ.  1.0  ) 


NSKP  -  NSKP 
TO  1200 


♦  NRFRST / 500 


GO  TO  311 


105 


) 

)  PASSN 
*RELF 

FORMAT  (  7F10.4  ) 

WRITE  (  6*1212  )  PASSN 
l  *  RELF 

FORMAT  (  10X, 7F 10. 4  ) 

WRITE  (  6,1020  ) 

W  -  WCELL 

IF  (  WCELL  .EQ.  4.0  )  GO  TO  1215 
W  -  W  /  25.4 


* WCELL  *XC 

*  YC 

, YBASE 

,  YTOP 

,RELT 

*  WCELL  *XC 

*  YC 

, YBASE 

*  YTOP 

*  RELT 

XC  «  XC  /  25.4 
YC  *  YC  /  25.4 

HO  YBASE  *  YBASE  /  25.4 

IF  (  YTOP  .LT.  9999.0  )  YTOP  -  YTOP  /  25.4 
——  5715  CONTINUE 

- "  WRITE  <  6*1220  )  , 

1220  FORMAT  (  65X*16H  HOR.  VERT.  ) 
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EXPLANATIONS  (CONTINUED) 


STOPPING  RUN. 

98-  103 

READ  AMD  WRITE  IMPUT  CARDS  4  AMD  5. 

10?-  Ill 

VJ  IS  NEEDED  BECAUSE  WCELL  NEEDS  TO  BE  PRESERVED  AS  READ  FROM  INPUT  CARO,  IF 
WCELL  =  100.0/  VI/  XC/  YC/  YBASE/  AND  POSSIBLY  YTOP  ARE  ALL  CONVERTED  FROM 
MILLIMETERS  TO  INCHES.  YTOP  =  <>999.0  OR  AMY  LARGER  NUMBER  INDICATES  YTOP  IS 
UNMEEDED  AND  TO  BE  IGNORED. 

113-  124 

REWRITE  INPUTS  ON  CARDS  4  AND  5  WITH  ABBREVIATED  INDICATION  OF  WHAT  EACH 

NUMBER  REPRESENTS. 

125 

H  IS  HALF  OF'  W  AND  CONSEQUENTLY  REPRESENTS  PERPENDICULAR  DISTANCE  FROM  CENTER 

OF  AMY  TARGET  CELL  TO  SIDES. 

126-  127 

VERTICAL  COORDINATE  YAIM  OF  INTENDED  AIMPOINT  CAN  CORRESPOND  EITHER  TO 
APPROXIMATE  CENTER  OF  MASS  OR  TO  BASE  OF  TARGET. 

128-  134 

VALUES  2/  3/  4,  AND  5  FOR  NADJST  HAVE  BEEN  USED  TO  DATE  ONLY  IN  CONNECTION 

WITH  SPECIAL  STUDY  NADE  FOR  US  ARMY  ARMOR  SCHOOL/  FORT  KNOX  IN  197e, 
REQUIREMENTS  OF  THAT  EFFORT  ARE  NOT  OF  GENERAL  INTEREST  AND  NEED  NOT  BE 
DISCUSSED  IN  DETAIL  IN  THIS  REPORT.  SINCE  5  IS  LARGEST  VALUE  USED  TO  DATE 

FOR  NADJST/  CONSIDERATION  IS  RESTRICTED  TO  NADJST  EQUALLING  0  OR  1  IN 
REMAINDER  OF  THESE  EXPLANATIONS. 

136-  158 

NSTCP.D  IS  NORMALLY  0  BECAUSE  MOST  RUMS  INVOLVE  TAPE  OR  DISC  CONTAINING  TARGET 
VULNERABILITY  DATA.  EXCEPTIONALLY/  INPUT  CARDS  CONTAINING  TARGET  SHAPE  DATA 
NEED  TO  BE  READ.  IF  MSTCRD  IS  NOT  0/  AT  LEAST  ONE  OF  CONTROL  QUANTITIES  IN 
KTRGTC  ARRAY  MUST  ALSO  DIFFER  FROM  0.  ITGT  VALUES  1/  2/  3/  AND  4  CORRESPOND 
TO  TARGET  ORIENTATION  ANGLES  OF  0/  30/  60/  AND  90  DEGREES  RESPECTIVELY.  ALL 
COORDINATES  ARE  IN  SAME  UNIT  AS  WCELL.  VERTICAL  COORDINATE  OF  CELL  CENTERS 
FOR  EACH  HORIZONTAL  STRIP  OF  TARGET  CELLS  IS  ENTERED  IN  NCTRY  ARRAY. 
HORIZONTAL  COORDINATES  OF  CENTERS  OF  CELLS  AT  LEFT  AMD  RIGHT  ENDS  OF  STRIP 
APE  STORED  IN  NCTRXL  AND  NCTRXR  ARRAYS.  LAST  COORDINATE  IN  EACH  SET  IS 
IDENTIFIED  BY  USE  OF  999«  AS  NEXT  ENTRY. 

160-  165 

INITIAL  SETTINGS. 

166-1382 

CYCLE  FOR  EACH  RANGE  CONSIDERED. 

166 

EXCEPT  WHEN  SET  TO  INITIAL  INTPL  SETTING/  CONTROL  QUANTITY  JNTPL  PRESERVES 
VALUE  OF  IMTPL  USED  IN  CALCULATION  CYCLE  FOR  PREVIOUS  RANGE  WHEN  RANGE 
INTERPOLATION  OF  TARGET  VULNERABILITY  DATA  MAY  HAVE  BEEN  INVOLVED.  INPUT 
VALUE  987  FOR  INTPL  IS  SPECIAL  SETTING  ENTERED  ONLY  ON  INPUT  CARD  R1  FOR 
FIRST  RANGE  CONSIDERED  IN  RUM  WHEN  MOVING  TARGET  OR  MOVING  FIRING  WEAPON  IS 
INVOLVED. 

167-  168 

BEGIN  PRINTING  FOR  EACH  RANGE  ON  NEW  PAGE. 

169-  170 

READ  AND  WRITE  INPUT  CARD  Rl. 

171-  1*’2 

CONTROL  QUANTITY  ISTMOV  IS  SET  TO  1  AND  COUNTER  MSET  IS  GIVEN  INITIAL  SETTING 
IF  RUN  INVOLVES  MOVING  TARGET  OR  MOVING  FIRING  WEAPON. 

174 

PROGRAM  LIMES  175  THROUGH  255  ARE  NORMALLY  APPLICABLE  BUT  ARE  SKIPPED  FOR 
MOVING  TARGET  OR  MOVING  FIRING  WEAPON  RUN. 

175-  245 

READ  AND  WRITE  INPUT  CARDS  P.2  THROUGH  R8.  REWRITE  INPUTS  ON  THESE  CARDS  WITH 
INDICATION  OF  WHAT  EACH  NUMBER  REPRESENTS.  OPTW  AND  NHTKLL  CONTROL  REWRITING 
OR  SKIPPING  OF  SEVERAL  QUANTITIES  THAT  ARE  OFTEN  UNNEEDED  FOR  RUN  AND  SET 

TO  0  (OR  LEFT  BLANK)  ON  INPUT  CARDS. 
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PROGRAM  LISTING  (CONTINUED) 


115 


120 


WRITE 

WRITE 


6# 1020 
6*1230 


125 


1230  FORMAT 
A 
B 
C 
D 
E 

_ F 

H  »  W  / 


PASSN 
*  RELT 
(  10X#  36H  PASSN 
10X#34H  W#W 
10X#36H  XC#YC 
10X# 36H  YBASE 
10X#  36H  YTOP 
10X#35H  RELT 
10X#33H  RELF 
2.0 


*  W 

#RELF 


#  W 


xc 


YC 


* YBASE  > YTOP 


TGT. 

CELL 

TGT. 

TGT. 


PERS.  EXCL.  CREW#  7X#F10.4  / 

DIM.  IN  INCHES#19X#2F10.4  / 

CTR.  COORD.  INS  •  # 17X#  2F10. 4  / 

BOT.  COORD.  INS •  #  27X#  F10.4  / 


TOP  APPR.  COORD.  INS.# 27X, F10. 4 
PROB.  OF  REL.  FLIGHT#8X» F10.4  / 
PROB.  OF  REL.  FUZE#10X#F10.4  ) 


130 

135 

1300 

140 

145 

i 

1 

i 

1330 

150 

i 

\ 

1332 

155 

1 

i 

1320 

1310 

YAIM  -  YC 
IF  (  NRFBS 
IF  (  NADJST 
XLEFTFl 1)  - 
XLEFTF ( 4 )  - 
XLEFTD(l)  - 
XRGHTF (1)  - 
XRGHTF (4)  - 
XRGHTD (1 )  - 
CONTINUE 
IF  (  NSTCRD  .EQ. 
WRITE  (  6#  1010  ) 

DO  1310  ITGT  -  1# 4 
IF  (  KTRGTC  tITGT ) 
WRITE  (  6# 1020  ) 

DO  1320  J  -  1# 30 
K  -  4  *  t  J— 1 ) 

READ  <  5# 1330 


.EQ.  2  )  YAIM  -  YBASE 
.LT.  2  .OR.  NADJST  .GT. 
-1650.0  /  25.4 
-3050.0  /  25.4 
-1350.0  /  25.4 
1650.0  /  25.4 
/  25.4 
/  25.4 


5  )  GO  TO  1300 


3550.0 

1450.0 


0  )  GO  TO  2000 


EQ.  0  )  GO  TO  2000 


FORMAT 
WRITE  ( 


(  1215  ) 
6#1332 


)  NCTRY(K+1# ITGT )#NCTRXL(K+1# ITGT )#NCTRXR<  K+l# ITGT ) 
#NCTRY<K+2#ITGT),NCTRXL (K+2#ITGT)#NCTRXR(K+2#I~GT) 
# NCTRY ( K+3#  ITGT)#NCTRXL(K  +  3#ITGT)#NCTRXR(K+3»I1‘GT  ) 
,NCTRY(K+4#ITGT)#NCTRXL(K+4#ITGT)#NCTRXR<K+4#I~GT  ) 


t 


)  NCTRY <K+1# ITGT)#NCTRXL<K+1#ITGT)#NCTRXR(K+1#ITGT) 
#  NCTR Y(K+2# ITGT)#NCTRXL<K  +  2# ITGT) # NCTRXRI K+2# I'GT ) 
#NCTRY( K+3# ITGT) #NCTRXL (K+3# ITGT) #NCTRXR( K+3# I TGT ) 
#NCTRY( K+4# ITGT)#NCTRXL(K+4#ITGT)#NCTRXR(K+4#I“GT ) 
10X#12I5  ) 


.EQ. 

.EQ. 

.EQ. 

.EQ. 


9999 

9999 

9999 

9999 


.OR. 
•  0  R. 
.OR. 

)  GO 


TO  1310 


160 


165 


TUTS 


2020 


170 


FORMAT 

IF  (  NCTRYIK+l# ITGT) 

A  NCTRYtK+2# ITGT ) 

B  NCTRY(K+3# ITGT) 

C  NCTRY( K+4# ITGT) 

CONTINUE 
CONTINUE 
CONTINUE 
NRANGE  -  1 
NEJECT  -  0 
JNTPL  -  0 
INTPL  *  0 
ISPLIT  -  0 
ISTMOV  -  0 

IF  (  NRDTYP  .EQ.  1  .AND.  INTPL  .NE.  987  )  JNTPL  -  INTPL 
WRITE  (  6#2020  ) 

FORMAT  (  1H1  ) 

READ  (  5# 1120  ) 

WRITE  (  6# 1122  ) 

IF  (  NRANGE  .EQ. 


I  RANGE#  I  NT  PL#  I  MI  LS 
IRANGE#INTPL#IMILS 
1  .AND.  INTPL  .EQ. 


987  )  ISTMOV  -  1 
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175 

180 

185 

190 

195 

200 

205 

210 

215 

220 

225 


PROGRAM  LISTING  (CONTINUED) 


IF  (  ISTMOV  .EQ.  1  )  MSET  -  1 
WRITE  (  6,1010  ) 

IF  (  ISTMOV  .EQ.  1  )  GO  TO  2200 


READ 

(  5,1210 

)  FLT 

,  TF1 

,  TFS 

, XM1  , XM2  ,  AMT1 

,  AMT2 

A 

,  STD1 

,  STD2 

,  XB 

, YB  , SIGXB  , S I6YB 

, SIGXL 

B 

, SIGYL  , SIGXR 

,SIGYR 

, PROBS  , SIGXS  , SIGYS 

,  PGH 

C 

,  PGS 

,  PGCH 

,  PGCS 

, SGHX  , SGHY  , SGSX 

,  SGSY 

D 

,CDRX 

,CDRY 

,HSX 

, HSY  , XBH  , YBH 

, SI GXBH 

E 

, SIGYBH,XBL 

,  YBL 

,SIGXBL,SIGYBL,XBS 

,  YBS 

F 

WRITE 

(  6,1212 

, SIGXBS, SIGYBS 
)  FLT  , TF1 

,TFS 

, XM1  , XM  2  , AMT1 

,  AMT2 

A 

,  STD1 

,  STD2 

,  XB 

, YB  , SIGXB  , SIGYB 

, SIGXL 

B 

, SIGYL  , SIGXR 

, SI GYR 

, PROBS  , SIGXS  , SIGYS 

,  PGH 

C 

,PGS 

,  PGCH 

,  PGCS 

, SGHX  , SGHY  , SGSX 

,  SGS  Y 

D 

,  C  DRX 

,  CDRY 

,HS  X 

, HSY  , XBH  , YBH 

, SIGXBH 

E 

, SIGYBH,XBL 

,  YBL 

,SIGXBL,SIGYBL,XBS 

,  YBS 

( 


( 


( 


WRITE 
WRITE 

2110  FORMAT 
A 

WRITE  ( 
WRITE  ( 

A 

2120  FORMAT 
A 
B 
C 
D 

WRITE  ( 

IF  (  IMILS 


6,1010  ) 
6,2110  ) 

63  X,20H 
/  64X,19H 

6,1020  ) 
6,2120  )  FIT 
,  STD1 

(  10X, 29H  FIT 
10X, 25H 
10X,35H 
10X, 27H 
10X, 33H 
6,1020  ) 


, SIGXBS, SIGYBS 


TIME  DATA  (SECONDS), 
1ST  RD.  SUBS.  RDS.  ) 


,  TF1 
,  STD2 


TF1, TFS 
XM1,  XM2 
AMT1, AMT2 
STD1,STD2 


, TFS  , XM1  , XM2  , AMT1  , AMT2 

TIME  OF  FLIGHT,14X,F10.4  / 

FIXED  TIME,28X,2F10. A  / 

MEDIAN  VARIABLE  TIME, 18X, 2F1 0. A 
MINIMUM  TIME,26X,2F10«4  / 

VARIABILITY  FACTOR, 20X, 2F10. 4  ) 


1  cn  Tn 


WRITE  (  6,2130  ) 

2130  FORMAT  (  62X,23H  ACCURACY  DATA  (METERS), 

A  /  66X, 16H  HOR.  VERT.  ) 

GO  TO  2135 


2125  WRITE  (  6,2140  ) 

2140  FORMAT  (  63 X, 2 1H  ACCURACY  DATA  (MILS), 

A  /  66 X, 16H  HOR.  VERT.  ) 

2135  WRITE  (  6,1020  ) 

WRITE  (  6,2150  )  XB  , YB  ,SIGXB  ,SIGYB  ,SIGXL  ,SIGYL  ,SIGXR 


A  ,S IGYR  , PROBS 

2150  FORMAT  (  10X,25H  XB, YB 

A  10X,34H  SIGXB, SIGYB 

B  10X, 34H  S IGXL, SIGYL 

C  10X, 34H  Si GXR, SI GYR 

D  10X, 36H  PROBS 

E  10X, 36H  S IGXS, SIGYS 

OPTW  -  PGH  ♦  PGS  ♦  PGCH  +  PGC 
A  ♦  CDRY  ♦  HSX  ♦  HSY 


, S IGXS  , SIGYS 
FIXED  BI AS, 2 8 X, 2 F10. 4  / 

VAR.  BIAS  STD.  DEV.,19X, 2F10.4  / 

LAY  ERROR  STD.  DEV., 19X, 2F10.4  / 

RD.-TO-RD.STD.  DEV.,19X,2F10.4  / 

PROB.  SENS.  MISS.  RD.,6X,F10.4  / 

SENS.  ERROR  STD.  DEV., 17X, 2F 10. 4  ) 

+  SGHX  ♦  SGHY  ♦  SGSX  +  SGSY  ♦  CDRX 


IF  (OPTW  .EQ.  0.0  )  GO  TO  2155 


WRITE  ( 

6,1010  ) 

WRITE  ( 

6,2160  )  PGH 

,PGS 

A 

,  SGSY 

,  CDRX 

2160  FORMAT 

(  10X,38H  PGH 

A 

10X, 37H  PGS 

B 

10X, 35H  PGCH 

C 

10X, 36H  PGCS 

, PGCH  , PGCS  , SGHX  , SGHY  ,SGSX 
,CDRY  ,HSX  , HS Y 
GUNNER  SENS.  PROB.  HIGH,5X,F10. 4 
GNR.  SENS.  PROB.  SHORT, 6X, F10.4 
G/C  SENS.  PROB.  HIGH,8X,F10.4  / 

G/C  SENS.  PROB.  SHORT, 7X, F 10. 4  / 
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PROGRAM  LISTING  (CONTINUED) 


230 


235 


240 


245 


D 

E 

F 

G 


10X, 35H 
10X,36H 
10X,33H 
10X,33H 


SGHX,SGHY 

SGSX,SGSY 

CDRX,CDRY 

HSX,HSY 


GNR.  SENS.  S.D.,  HIGH*  18X,2F10.4 
GNR.  SENS.  S.D.  SHORT, 17X, 2F10. 4 
CMDR.-TO-GNR.  S.D ., 20X, 2F10. 4  / 

HIT  ADJ.  STD.  DEV., 20X, 2 FI 0.4  ) 


2155 


CONTINUE 
IF  {  NHTKLL 


,EQ .  0  )  GOTO  2165 


WRITE 

WRITE 

WRITE 

WRITE 


6,1020 

6,2140 

6,1020 

6,2170 


XBH  ,  YBH 
, S I GYBL, XBS  , YBS 
10X,23H  XBH, YBH  BIAS 

SIGXBH, SIGYBH  STD. 
XBL, YBL  BIAS 

SIGXBL,SIGYBL  STD. 
XBS, YBS  BIAS 

S IGXBS, SIGYBS  STD. 


, SIGXBH, SIGYBH, XBL  , YBL  ,SIGXBL 


250 


2170  FORMAT  { 

A 
B 
C 
D 
E 

2165  CONTINUE 

IF  (  INTPL  .EQ.  0  )  GO  TO  2175 
INCRNG  -  500 

IRNG2  *  IRANGE/INCRNG  ♦  1 
DELRNG  -  IRNG2*INCRNG  -  IRANGE 


10X, 28H 
10X, 29H 
10X, 34H 
10X, 31H 
10X, 36H 


,SIGXBS, SIGYBS 
HIT,30X,2F10.4  / 

DEV.  HIT , 25X, 2F10 . 4  / 

LOST  MISS,24X,2F10.4  / 

DEV.  LOST  MISS,19X,2F10.t  I 
SENSED  MISS, 22X,2F 10.4  i 
DEV.  SENSED  MISS, 17X, 2F 10. 4 


2175 


2206 

260 

2220 

265 

2222 

2210 

270 

2230 

2240 

275 

2300 

2305 

TOT 


285 


RATIO  *  DELRNG  /  500.0 
CONTINUE 

IF  (  STD1  .EQ.  0.0  )  STD1  ■  0.4983 
IF  {  STD2  .EQ.  0.0  )  STD2  *  0.4983 
GO  TO  2300 

READ  (  5,1210  )  FLT  , TF1  ,TFS 
WRITE  (  6,1212  )  FLT  ,TF1  ,TFS 
STD1  -  0.0 
STD2  *  0.0 

WRITE  (  6,1212  )  STD1  ,STD2 
DO  2210  M  *  1,3 
DO  2210  N  »  1,6 

READ  (  5,2220  )  XXB(M,N)  ,YYB(M,N) 
FORMAT  {  20X, 4F10. 4  ) 

WRITE  (  6,2222  )  XXB(M,N)  , YYB ( M, N ) 
FORMAT  (  30X,4F10.4  ) 

CONTINUE 
SIGXL  «  0.0 

WRITE  (  6,2230  )  SIGXL 
FORMAT  (  70X,F10.4  ) 

READ  (  5,2240  )  SIGYL  ,SIGXR  ,SIGYR 

FORMAT  (  6F10.4  /  /  /  /  ) 

WRITE  (  6,1212  )  SIGYL  »SIGXR  ,SIGYR 

CONVRT  »  1.0 
PI  *  3.14159 
R  *  IRANGE 

IF  (  IMILS  .EQ.  0  )  GO  TO  2305 

CONVRT  «  R*PI  /  3200.0 

CONVRT  ■  CONVRT  *  39.37 

IF  (  ISTMOV  .EQ.  1  )  GO  TO  3000 

XB  *  XB  *  CONVRT 

YB  *  YB  *  CONVRT 

SIGXB  ■  SIGXB  *  CONVRT 

SIGYB  *  SIGYB  *  CONVRT 

SIGXL  *  SIGXL  *  CONVRT 


*  XM1 
,XM1 


,  XM2 
,  XM2 


,  AMT1 
,  AMT1 


,  AMT2 
,  AMT  2 


,SSIGXB (M,N ) 
,SSIGXB(M,N ) 


, SSI GYB { M»  N ) 
,SSIGYB { M,N ) 


, PROBS  , SIGXS  , S I GYS 
, PROBS  , SIGXS  , S IGYS 
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EXPLANATIONS  (CONTINUED) 


247-  251 

IF  needed,  range  interpolation  factor  ratio  is  DETERMINED.  FOR  EXAMPLE*  RANGES 
OF  600  AND  2350  METERS  WOULD  CORRESPOND  TO  0.8  AND  0.3  RESPECTIVELY.  PRODUCT 
OF  RATIO  AND  500  IS  DIFFERENCE  IN  METERS  BETWEEN  RANGE  IPANGE  AND  NEXT 

HIGHER  MULTIPLE  OF  500  METERS. 

253-  254 

FOR  CONVENIENCE*  ZEROS  CAN  BE  USED  FOR  STD1  AND  STD2  ON  INPUT  CARD  R3  INSTEAD 
OF  0.4983  WHICH  IS  GENERALLY  APPLICABLE.  IF  THIS  HAS  BEEN  DONE*  RESET  STD1 
AND  STD2. 

255 

SKIP  PROGRAM  LINES  256  THROUGH  273. 

256-  273 

READ  AMD  WRITE  INPUT  CARDS  CONTAINING  FIRST  ROUND  BIAS  AND  STANDARD  DEVIATION 
DATA  NEEDED  FOP.  MOVING  TARGET  OR  MOVING  FIRING  WEAPON  RUN.  SINCE  TIMES  ARE 
NOT  INVOLVED*  FIRST  CARD  CAN  CONTAIN  ONLY  ZEROS  AND  STD1  AND  STD2  ARE  SET 

TO  0.  DO  2210  LOOP  PROCESSES  SET  OF  18  INPUT  CARDS  AND  STORES  DATA  IN  XXB* 
YYB*  SSIGXB*  AND  SSIGYB  ARRAYS.  SINCE  STANDARD  DEVIATIONS  ENTERED  IN  LATTER 
TWO  ARRAYS  ALREADY  ACCOUNT  FOR  RANDOM  ERRORS  AS  WELL  AS  OTHER  ERROR 
CONTRIBUTORS*  SIGXL  IS  SET  TO  0.  LAST  CARD  CAN  CONTAIN  ONLY  ZEROS  BECAUSE 
QUANTITIES  RELATED  TO  SUBSEQUENT  ROUND  SENSING  AMD  ADJUSTMENT  AS  WELL  AS  TO 
RANDOM  ERRORS  ARE  UNNEEDED.  THE  18  SETS  OF  INPUT  BIASES  AND  STANDARD 
DEVIATIONS  MAY_  INCLUDE  DUMMY  SETS  FOP.  WHICH  CALCULATIONS  ARE  NOT  TO  BE  MADE. 
SUCH  DUMMY  SETS  ARE  IDENTIFIED  BY  HORIZONTAL  BIAS  SETTING  OF  999.9999. 

274-  27? 

CONVERSION  FACTOR  CONVRT  IS  SET  TO  APPROPRIATE  VALUE  FOR  CHANGING  METERS  OR 
MILS  TO  INCHES* 

280 

SKIP  PROGRAM  LINES  281  THROUGH  328  FnR  MOVING  TARGET  OR  MOVING  FIRING  WEAPON 
RUN . 

281-  298 

CONVERT  FROM  METERS  OR  MILS  t0  INCHES. 

299 

SKIP  PROGRAM  LINES  300  THROUGH  328  IF  ASSOCIATED  QUANTITIES  ARE  UNNEEDED. 

300-  312 

CONVERT  FROM  MILS  TO  INCHES.  IF  IMILS  EQUALS  0*  PREVIOUS  SETTING  OF  CONVPT  HAS 
NOT  PROVIDED  FOR  CONVERSION  FROM  MILS  TO  METERS. 

313-  328 

ENTER  BIASES  AND  STANDARD  DEVIATIONS  IN  BBX*  BBY*  SSIGX*  AND  SSIGY  ARRAYS. 

330 

SKIP  PROGRAM  LINES  331  THROUGH  462  IF  TARGET  SHAPE  DATA  FROM  CARDS  ARE  NEEDED. 

331-  462 

CYCLE  FOR  READING  AMD  PROCESSING  TARGET  VULNERABILITY  DATA  FROM  TAPE  OR  DISC. 

331-  334 

DO  NOT  READ  AND  PROCESS  VULNERABILITY  DATA  FROM  TAPE  OR  DISC  WHEN  NEEDED 

INFORMATION  HAS  ALREADY  BEEN  READ  AMD  PROCESSED.  CONDITION  INVOLVING  NFANGE* 
IMTpL*  AND  IS°L IT  COVERS  INSTANCES  WHERE  FIRST  RANGE  OF  CONCERN  FOR  RANGE 
INTERPOLATION  OF  VULNERABILITY  DATA  HAS  BEEN  INVOLVED  IN  PREVIOUS  RANGE 

CYCLE.  JNTPL  EQUALLING  1  INDICATES  THAT  SECOND  RANGE  USED  IN  INTERPOLATION 

OF  IMMEDIATELY  PRECEDING  RANGE  CYCLE  IS  RANGE  OF  CONCERN  FOR  CURRENT  CYCLE. 
VULNERABILITY  DATA  THAT  DO  NOT  VARY  WITH  RANGE  ARE  READ  AND  PROCESSED  IN 
FIRST  RANGE  CYCLE  ONLY  AMD  CAM  THEM  BE  REUSED. 

335 

CONDITION  IS  ASSOCIATED  WITH  RANGE  INTERPOLATION  OF  TARGET  VULNERABILITY  DATA. 
ISPLIT  EQUALLING  2  INDICATES  THAT  CALCULATIONS  HAVE  ALREADY  BEEN  COMPLETED 
FOR  LESSER  OF  TWO  RANGES  BRACKETING  ACTUAL  RANGE  AND  PARTICULAR  PAIR  OF 
FORWARD  AMD  REVERSE  ORIENTATION  ANGLES.  IF  FORWARD  ANGLE  OF  0  DEGREES  IS 
INVOLVED*  CONTINUE  WITH  PROGRAM  LINES  336  THROUGH  374  BEFORE  OBTAINING 
VULNERABILITY  DATA  FOR  RANGE  BRACKETING  THAT  OF  TARGET  ON  HIGH  SIDE  AND  FOR 
ANGLES  OF  0  AND  180  DEGREES.  IF  FORWARD  ANGLE  IS  30*  60*  OR  90  DEGREES*  SKIP 
TO  PROGRAM  STATEMENT  320C. 
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PROGRAM  LISTING  (CONTINUED) 


290 


295 


300 


305 


310 


315 


320 


325 


SIGYL  - 

SIGYL 

*  CONVRT 

SIGXR  - 

SIGXR 

*  CONVRT 

SIGYR  - 

SIGYR 

*  CONVRT 

SIGXS  - 

SIGXS 

*  CONVRT 

SIGYS  - 

SIGYS 

*  CONVRT 

SGHX  - 

SGHX  * 

CONVRT 

SGHY  - 

SGHY  * 

CONVRT 

SGSX  - 

SGSX  * 

CONVRT 

SGSY  - 

SGSY  * 

CONVRT 

CDRX  - 

CDRX  * 

CONVRT 

CDRY  - 

CDRY  * 

CONVRT 

HSX 
HSY 
IF  ( 
IF  ( 
XBH 
YBH 


■  HSX  *  CONVRT 
«  HSY  *  CONVRT 

NHTKLL  .EQ.  0  )  GOTO  3000 

IMILS  .EQ.  0  )  CONVRT  ■  CONVRT  *  R+PI/3200.0 

■  XBH  *  CONVRT 
•  YBH  *  CONVRT 

SIGXBH  ■  SIGXBH  *  CONVRT 
SIGYBH  =■  SIGYBH  *  CONVRT 
XBL  -  XBL  *  CONVRT 
YBL  -  YBL  *  CONVRT 
SIGXBL  *  SIGXBL  *  CONVRT 
SIGYBL  -  SIGYBL  *  CONVRT 
XBS  -  XBS  *  CONVRT 
YBS  -  YBS  *  CONVRT 


330 


SIGXBS  - 
SIGYBS  - 
BBX(l)  - 
BBY(l)  - 
BBX (2)  * 
BBY  (2 )  - 
BBX  (3 )  - 
BBY (3 )  * 
BBX (4 )  = 
BBY(4)  - 
SSIGX(l) 
SSIGY(l) 
SSIGXI2 ) 
SSIGYI2 ) 
SSIGXt  3 ) 
SSIGY(3) 
SSIGX<4> 
SS IGY (4 ) 
CONTINUE 
IF 
IF 


SIGXBS  *  CONVRT 

SIGYBS  *  CONVRT 

XB 

YB 

XBH 

YBH 

XBL 

YBL 

XBS 

YBS 

-  SORT (SI GXB**2+SIGXL**2+SIGXR**2 ) 

*  SQRT(SIGYB**2  +  SIGYL**2«-SIGYR**2) 

-  SIGXBH 

-  SIGYBH 

*  SIGXBL 

•  SIGYBL 

•  SIGXBS 

-  SIGYBS 


)  GO  TO  3800 
.AND.  INTPL  . 


ITT 


1  )  GO  TO  4010 
.  0  .AND.  NRANGE 
.AND.  NANGLE 


.GT. 

.NE. 


340 


{  NSTCRD  .GT.  0 

{  NRANGE  .GT.  1 

A  4010 
IF  (  JNTPL  .EQ. 

IF  (  NRDTYP  .EQ. 

IF  {  ISPLIT  .EQ. 

-  WRITE  (  6/1010  ) 

IF  (  NEJECT  .EQ. 

3110  FORMAT  (  80A1  ) 

IF  (  EOF(IU)  .EQ.  1.0  )  GO  TO  301 
301  CONTINUE 

IF  (  ISTMOV  .EQ.  1  .AND.  NRANGE  .NE. 
IF  (  EOF(IU)  .EQ.  1.0  )  GO  TO  312 


EQ.  1  .AND.  ISPLIT  .EQ.  0  )  GO  TO 


GO 

GO 


TO 

TO 


4000 

3200 


1  )  READ  (  IU# 3110  )  A 


1  )  READ  (  I U/ 3110  )  A 
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PROGRAM  LISTING  (CONTINUED) 


345 

302 

3112 

350 

355 

3120 

3130 

360 

3150 

3140 

365 

3170 

3160 

370 

3180 

1 

i 

375 

3200 

380 

303 

385 

304 

3210 

390 

GO  TO  3120 
)  GO  TO  3120 
)  GO  TO  3120 


312  CONTINUE 

READ  (  IU/3110  )  ITAPE 
IF  {  EOF { IU)  .EQ.  1.0  )  GO  TO  302 
CONTINUE 

WRITE  {  6/3112  )  ITAPE 
FORMAT  {  10X/80A1  ) 

DO  3120  K  -  1/7 A 
IF  {  ITAPE (K )  .NE.  IDCODE(A)  ) 

IF  {  ITAPEtK+1)  .NE.  IDC0DEI5) 

IF  {  IT  APE (K+2)  .NE.  IDC0DEI6) 

M  *  K+3 
GO  TO  3130 
CONTINUE 
GO  TO  3140 
DO  3150  L  -  M/77 
IF  {  ITAPEU)  .NE.  IDCODE  (7 )  ) 

IF  (  ITAPEtL+l )  .NE.  IDC0DEI8) 

IF  {  ITAPEIL+2)  .NE.  IDC0DE<9) 

IF  {  ITAPE (L+3)  .NE.  IDCODEUO) 

GO  TO  3160 
CONTINUE 

WRITE  {  6/3170  ) 

FORMAT  {  1H0/  10X/32H  PROJECTILE/TARGET  CODE  MISMATCH  ) 
GO  TO  9900 
CONTINUE 

WRITE  (  6/3180  )  {  IDCODE  { I )/  1*1/  3  l/{  IDCODE  (I )/ 1-4/6 
k  /(  IDCODE (I)/I*7/10  ) / IRANGE 

FORMAT  {  24X/26H  TASK  SHEET  IDENTIFICATI0N/10X/3A1  / 

A  24X/16H  PROJECTILE  C0DE/20X/3A1  / 

B  24X/12H  TARGET  C0DE/23X/4A1  / 

C  v  24X/16H  RANGE  IN  METERS/ 18X/I5  ) 


GO  TO  3150 
)  GO  TO  3150 
)  GO  TO  3150 
)  GO  TO  3150 


2  )  IC 


IANGLE 


1.0  )  GO  TO  303 


395 


3400 


3410 


3405 


24X/16H 
24X/12H 
24X/16H 
IANGLE  -  1 
IF  {  ISPLIT  .NE 
IB  -  IC  *  2 
IA  •  IB  -  1 
WRITE  {  6/1010  ) 

READ  {  IU/3110  )  i 
IF  {  EOF { IU)  .EQ. 

CONTINUE 
NSKIP  -  3 
WRITE  {  6/1020  ) 

DO  3210  L  -  1/ NSKIP 

READ  {  IU/3110  )  A 

IF  {  EOF { IU)  .EQ.  1.0  )  GO  TO  304 

CONTINUE 

WRITE  (  6/3112  )  A 
CONTINUE 

YTEMPL { IC )  -  9999.0 
YTEMPH  »  -9999.0 
XTEMPLtlC)  -  9999.0 
XTEMPH(IC)  -  -9999.0 
M  -  1 

IF  {  NTGT  .GT.  1  )  GO  TO  3405 
READ  {  IU/ 3410  )  X/Y/t  PKRLFTt J) / J-l/4  )/{ 
FORMAT  {  1X/F6.1/F7.1/5X/4F7.3/6X/4F7.3  ) 
GO  TO  305 

IF  {  NTGT  .GT.  2)  GO  TO  3415 


PKRLFTU  J/J-6/9  ) 
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EXPLANATIONS  (CONTINUED) 


337-  343 

SPECIAL  STATEMENTS  NEEDED  TO  READ  AND  IGNORE  BLANK  LINE  OR  EJECT  PAGE  LINE 

PERIODICALLY  USED  ON  TAPE  OR  DISC. 

344-  373 

READ  NEXT  LIME  OF  INFORMATION  FROM  TAPE  OR  DISC.  THIS  LINE  SHOULD  CONTAIN 
PP.OJECTILE  CODE*  RANGE*  AND  TARGET  CODE.  VERIFY  AGREEMENT  OF  PROJECTILE  AMD 
TARGET  CODES  WITH  THOSE  IN  IDCODE.  IF  DISCREPANCY  IS  FOUND*  PRINT  ERROR 
INDICATION  AND  GO  TO  PROGRAM  STATEMENT  9900  TO  STOP  RUN.  OTHERWISE  WRITE 
SPECIFIED  REFERENCE  INFORMATION. 

374 

INITIAL  SETTING. 

375-  460 

CYCLE  FOR  I ANGLE  EQUALLING  1*  2*  3*  OR  4.  THESE  VALUES  CORRESPOND  TO  TARGET 
ORTENTATION  ANGLES  OF  0*  30*  60*  AMD  90  DEGREES  RESPECTIVELY  TOGETHER  WITH 
ASSOCIATED  REVERSE  ANGLES. 

375-  377 

INITIAL  SETTINGS. 

379-  35° 

READ  BLANK  LINE  FROM  TAPE  OR  DISC.  ALSO  READ  AND  WRITE  THREE  LINES  OF  HEADER 
INFORMATION  SUCH  AS  0  DEG  ELEV*  30  DEG  AZIM  ETC... 

3°0-  394 

INITIAL  SETTINGS. 

395-  447 

CYCLE  FOP.  EACH  TARGET  CELL. 

395-  405 

READ  CELL  CENTER  COORDINATES  AND  RELATED  VULNERABILITY  DATA.  FORMAT  PREVIOUSLY 
USED  TO  STORE  NUMBERS  ON  TAPE  OR  DISC  IS  IDENTIFIED  BY  NTGT .  LABEL  CELL 
COORDINATES  AS  X  AMD  Y.  WCELL  INDICATES  WHETHER  X  AND  Y  ARE  IN  INCHES  OR 
MILLIMETERS.  STORE  VULNERABILITY  DATA  IN  PKRLFT  ARRAY*  USING  PKRLFT(l)  AND 
PKRLFT(  6)  FOP.  M  (MOBILITY)  KILL  PROBABILITIES*  PKP.LFT(2)  AND  (7)  FOR  F 
(FIREPOWER)  KILL  PROBABILITIES*  (3)  AND  (8)  FOP  M  OR  F  (MOBILITY  AND/OR 
FIREPOWER)  KILL  PROBABILITIES*  (4)  AND  (9)  FOR  K  (COMPLETE  DESTRUCTION)  KILL 
PROBABILITIES  AND  IF  NECESSARY  (5)  AND  (10)  FOR  PASSENGER  PERSONNEL 
CASUALTIES.  PKRLFT(l)  THROUGH  (5)  ARE  ASSOCIATED  WITH  TARGET  ORIENTATION 
ANGLES  OF  0*  30*  60*  AND  90  DEGREES*  AND  PKRLFT(6)  THROUGH  (10)  WITH 
CORRESPONDING  REVERSE  ANGLES. 

406-  411 

NUMBER  READ  AS  X  COORDINATE  OF  CELL  CENTER  MAY  EXCEPTIONALLY  BE  99°. 9  USED  ON 
TAPE  OR  DISC  TO  INDICATE  ALL  DATA  FOR  PARTICULAR  SET  OF  FORWARD  AND  REVERSE 
ANGLES  HAVE  BEEN  READ.  IN  THAT  CASE*  MULTIPLICATION  OF  X  BY  10.0  SETS  I END 

Tf)  Q99Q,  ADDITION  OF  10000.001  TO  X  AND  Y  ELIMINATES  ALL  MINUS  SIGNS. 

INTEGER  QUANTITIES  IMODX  AND  IMODY  ARE  READILY  COMBINED  INTO  SINGLE  ENTRY 
STORABLE  IN  IMODXY  ARRAY.  WHEN  X  EQUALS  °°°.°*  IMODX  IS  SET  TO  1999°  AND 
COMBINED  WITH  Y  COORDINATE  0  BEFORE  STORAGE  OF  RESULT  IN  IMODXY  ARRAY.  AFTER 
END  INDICATOR  IS  THUS  STORED*  RUN  CONTINUES  AT  PROGRAM  STATEMENT  370C. 

412—  427 

IGNORE  THESE  SPECIAL  STATEMENTS  FOR  NADJST  EQUALLING  2  THROUGH  5. 

428-  430 

IF  CELL  COORDINATES  ARE  IN  MILLIMETERS*  CONVERT  TO  INCHES. 

431-  434 

UPDATE*  AS  NECESSARY*  MINIMUM  AND  MAXIMUM  COORDINATES  RELATED  TO  PECTANGLE 
ENCLOSING  TARGET.  IF  PREVIOUSLY  ESTABLISHED  VALUE  OF  YTEMPL(IC)  IS  GREATER 
THAN  Y*  RESET  YTEMPL(IC)  TO  Y.  SIMILARLY*  RESET  YTEMPH  IF  LESS  THAN  Y* 
XTEMPL(IC)  IF  GREATER  THAN  X*  AND  XTEMPH(IC)  IF  LESS  THAN  X. 

435 

IF  NPRHIT  EQUALS  1*  SKIP  PROGRAM  LIMES  436  THROUGH  447  SINCE  THEY  INVOLVE 
UNHEEDED  TARGET  VULNERABILITY  DATA. 

436-  447 

QUANTITIES  IN  PKRLFT  ARRAY  CONTAIN  THREE  DECIMAL  PLACES.  CONVERT  TO  CONVENIENT 
INTEGER  FORM  AND  ENTER  IN  IP  ARRAY.  PACK  FIRST  FIVE  IP  VALUES  INTO  SINGLE 
NUMBER  AND  STORE  RESULT  IN  MPK  ARRAY.  TARGET  CELL  INVOLVED  IS  IDENTIFIED 

BY  M.  IA  EQUALS  1*  3*  5*  OR  7  FOR  TARGET  ORIENTATION  ANGLE  OF  0*  30*  60*  OR 
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PROGRAM  LISTING  (CONTINUED) 


400 


405 


410 


READ  (  IU/3420  )  X, Y/PKRLFT 
3420  FORMAT  (  1X/2F6.0/2 13X/4F6.2/F7.2)  ) 

GO  TO  305 

3415  READ  (  IU/3430  )  X# Y/PKRLFT 
3430  FORMAT  (  IX/ F6. 1,F7. 1, 5X, 2( 4F7.3, F6. 1 ) 
305  IF  (  EOF ! IU)  .EQ.  1.0  )  GO  TO  3500 
3500  IEND  -  X  *  10.0 

IMODX  -  X  +  10000.001 
IF  (  IEND  .EQ.  9990  ) 

IMODY  -  Y  +  10000.001 
IMODXY (M/IC )  ■  IMODXnOOOOO  ♦  IMODY 


IMODX  ■  19999 


415 


420 


425 


IF 

1 

IEND  .EQ.  9999  )  GO 

TO  3700 

TO 

IF 

( 

NADJST 

.EQ. 

2  .AND. 

IANGLE  .EQ. 

2 

)  GO 

3600 

IF 

I 

NADJST 

.EQ. 

3  .AND. 

IANGLE  .EQ. 

2 

)  GO 

TO 

3600 

IF 

( 

NADJST 

.EQ. 

4  .AND. 

IANGLE  .EQ. 

2 

)  GO 

TO 

3600 

IF 

( 

NADJST 

.EQ. 

5  .AND. 

IANGLE  .EQ. 

2 

)  GO 

TO 

3600 

IF 

I 

NADJST 

.EQ. 

2  .AND. 

IANGLE  .EQ. 

3 

)  GO 

TO 

3600 

IF 

( 

NADJST 

.EQ. 

3  .AND. 

IANGLE  .EQ. 

3 

)  GO 

TO 

3600 

IF 

1 

NADJST 

.EQ. 

4  .AND. 

IANGLE  .EQ. 

3 

)  GO 

TO 

3600 

IF 

( 

NADJST 

.EQ. 

5  .AND. 

IANGLE  .EQ. 

3 

)  GO 

TO 

3600 

) 

GO 

IF 

1 

NADJST 

.EQ. 

2  .AND. 

IANGLE  .EQ. 

A 

.AND 

.  YBASE 

.GT. 

o 

• 

o 

T 

A 

i 

0  3600 

.GT. 

) 

GO 

IF 

1 

NADJST 

.EQ. 

3  .AND. 

IANGLE  .EQ. 

A 

.AND 

.  YBASE 

o 

• 

o 

T 

A 

0  3600 

.GT. 

GO 

IF 

1 

NADJST 

• 

m 

& 

• 

4  .AND. 

IANGLE  .EQ. 

A 

.AND 

.  YBASE 

o 

• 

o 

) 

T 

A 

0  3600 

.GT. 

) 

GO 

IF 

! 

NADJST 

.EQ. 

5  .AND. 

IANGLE  .EQ. 

A 

.AND 

•  YBASE 

o 

• 

o 

T 

430 


3505 


435 


^  0  3600 

IF  (  WCELL  .EQ.  4.0  )  GO  TO  3505 
X  -  X  /  25.4 
Y  -  Y  /  25.4 

YTEMPL(IC)  -  AMIN1 { YTEMPL ( IC )/ Y) 
YTEMPH  »  AMAX1 ( YTEMPH/ Y ) 

XTEMPL  t  IC )  -  AMINUXTEMPL  ( IC )/ X ) 
XTEMPHIIC)  *  AMAX1 { XTEMPH! I C)»  X ) 

- IF  (  NPRHIT  .EQ.  1  )  GO  TO  3600 

DO  3510  K  -  1/10 

3510  IP (K)  -  1000. 0*PKRLFT (K )  ♦  0.1 


440 


445 


4  50 


"JEM 


3605 

3610 


T5T 


~mo 


MPKtM/IA)  -  SHIFT (IP(l)/49)  .OR. 

SHIFT(IP(2)/38)  .OR. 

SHI  FT ( IP( 3) / 27)  .OP. 

SHIFT! IP(4)/16)  .OR. 

IP  t  51 

MPKIM/IB)  *  SHIFT! IP( 6) /49)  .OR. 

SHIFT i I P( 7) / 38 )  .OR. 

SHIFT! IP!8)/27)  .OR. 

SHIFT i IP i 9) ,16)  .OR. 

IP! 10) 

CONTINUE 
M  »  M  +  1 

IF  (  M  .GT.  3000  )  GO  TO  3605 
GO  TO  3400 
WRITE  I  6/3610  ) 

FORMAT  (  1H0,  10X/22H  TARGET  MATRIX  TOO  BIG  ) 

GO  TO  9900 

CONTINUE 

NXTEMP ! IC )  »  !XTEMPHIIC)-XTEMPL!IC))/W  ♦  1.001 
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EXPLANATIONS  (CONTINUED) 


90  DEGREES  RESPECTIVELY.  SIMILARLY  PACK  LAST  FIVE  IP  VALUES  FOR  STORAGE  IN 
MPK  ARRAY.  IB  IS  ASSOCIATED  WITH  REVERSE  ANGLES  AND  EQUALS  1  MORE  THAN 
CORRESPONDING  IA  VALUE. 

448-  454 

RESET  CELL  COUNTER  M.  IF  M  THEN  EQUALS  3001,  3000  CELLS  HAVE  ALREADY  BEEN 
PROCESSED  AND  MAXIMUM  ALLOWED  29«>9  CELLS  HAS  BEEN  EXCEEDED.  IF  TARGET 

CONTAINS  EXACTLY  2999  CELLS/  CELL  3000  IS  ASSOCIATED  WITH  099.9  ENTRY.  IF 
LIMIT  HAS  BEEN  VIOLATED,  PRINT  ERROR  INDICATION  AND  GO  TO  PROGRAM  STATEMENT 
9900  TO  STOP  RUM.  OTHERWISE  RECYCLE  FOP.  NEXT  CELL. 

456-  457 

CALCULATE  NUMBER  OF  CELLS  ALONG  EACH  SIDE  OF  RECTANGLE  THAT  IS  JUST  LARGE 
ENOUGH  TO  ENCLOSE  TARGET  FOR  FORWARD  AND  REVERSE  ORIENTATIONS  CORRESPONDING 
TO  IC.  RECTANGLE  HAS  NXTEMP(IC)  CELLS  IN  HORIZONTAL  DIRECTION  AND  NYTEMP(IC) 
CELLS  IN  VERTICAL  DIRECTION.  NUMBERS  OF  CELLS  BASED  ON  DIFFERENCES  BETWEEN 
MINIMUM  AND  MAXIMUM  COORDINATES  OF  CELL  CENTERS  ARE  INCREASED  BY  1.000  TO 
ACCOUNT  FOR  HALF  CELLS  AT  EACH  EDGE.  ADDITION  OF  .001  COVERS  ROUNDING 
INACCURACIES. 

458 

WHEN  RANGE  INTERPOLATION  OF  VULNERABILITY  DATA  IS  INVOLVED  AND  TARGET  RANGE  IS 
NOT  FIRST  ONE  CONSIDERED  IN  RUN,  CONTINUE  AT  PROGRAM  STATEMENT  4100  FOR  EACH 
IANGLE  VALUE. 

459 

CONTINUE  AT  PROGRAM  STATEMENT  4000  IF  IANGLE  VALUES  1  THROUGH  4  HAVE  ALL  BEEN 
COVERED. 

460 

IF  THIRD  OF  FOUR  INTERPOLATION  STEPS  CORRESPONDING  TO  ISPLIT  IS  OF  CURRENT 
CONCERN,  CONTINUE  AT  PROGRAM  STATEMENT  4020  FOR  PARTICULAR  IANGLE  VALUE 

UNDER  CONSIDERATION. 

461-  462 

SINCE  IANGLE  VALUES  1  THROUGH  4  HAVE  NOT  ALL  BEEN  PROCESSED,  REPEAT  FOP  NEXT 
VALUE. 

463-  505 

CYCLE  FOR  PROCESSING  TARGET  SHAPE  DATA  READ  FROM  CARDS. 

464 

INITIAL  SETTING. 

465-  471 

INITIAL  SETTINGS.  JCELL  IS  INTEGER  FORM  OF  WCELL. 

472-  473 

IT  VALUES  1>  2,  3,  AND  4  CORRESPOND  TO  TARGET  ORIENTATION  ANGLES  OF  0,  30,  60, 
AMD  90  DEGREES  RESPECTIVELY.  COUNTER  K  INDICATES  WHICH  STRIP  OF  TARGET  CELLS 
IS  OF  CURRENT  CONCERN.  STRIP  CAM  EXCEPTIONALLY  CONSIST  OF  ONLY  ONE  CELL.  SET 
JY  TO  EQUAL  Y  COORDINATE  COMMON  T0  CENTERS  OF  ALL  CELLS  IN  STRIP.  SET  JX  TO 
EQUAL  X  COORDINATE  OF  CENTER  OF  CELL  AT  LEFT  EDGE  OF  STRIP. 

474-  477 

ADDITION  OF  10000  TO  JX  AND  JY  ELIMINATES  ALL  MINUS  SIGNS.  COMBINE  IMODX  AMD 
IMODY  AND  STORE  RESULTING  QUANTITY  IN  IMODXY  ARRAY.  WHEN  JX  EQUALLING  Q9Q9 
INDICATES  THAT  ALL  TARGET  CELLS  FOR  ANGLE  CORRESPONDING  TO  IT  HAVE  BEEN 
PROCESSED,  RUN  CONTINUES  AT  PROGRAM  STATEMENT  3900. 

4^8-  482 

IF  CELL  COORDINATES  ARE  IN  MILLIMETERS,  CONVERT  TO  INCHES. 

483-  486 

UPDATE,  AS  NECESSARY,  MINIMUM  AND  MAXIMUM  COORDINATES  RELATED  TO  RECTANGLE 
ENCLOSING  TARGET.  IF  PREVIOUSLY  ESTABLISHED  VALUE  OF  YTEMPL(IT)  IS  GREATER 
THAN  Y,  RESET  YTEMPL(IT)  TO  Y.  SIMILARLY,  RESET  YTEMPH  IF  LESS  THAN  Y, 

XTEMPH IT1  IF  GREATER  THAN  X,  AND  XTEMPH(IT)  IF  LESS  THAN  X. 

487 

RESET  CELL  COUNTER. 

488 

CELL  COUNTER  SHOULD  NEVER  EXCEED  3000.  IF  LIMIT  HAS  BEEN  VIOLATED,  GO  TO 
PROGRAM  STATEMENT  3845. 
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PROGRAM  LISTING  (CONTINUED) 


NYTEMP(IC)  ■  (YTEMPH-YTEMPLtlC)  )/W  «■  1.001 

-  IF  (  INTPL  .EQ.  1  .AND.  NRANGE  .GT.  1  )  GO  TO  *100 

-  IF  (  IANGLE  .EQ.  4  )  GO  TO  4000 

460  IF  {  ISPLIT  .EQ.  2  )  GO  TO  4020 

-  IANGLE  -  IANGLE  +  1 

GO  TO  3200 
3800  CONTINUE 
IT  *  1 

465  3805  YTEMPL(IT)  -  9999.0 

YTEMPH  -  -9999.0 
XTEMPL ( IT )  -  9999.0 
XTEMPH ( IT )  -  -9999.0 
M  -  1 

470  K  -  1 

JCELL  -  WCELL 

3810  JY  ■  NCTRY(K,IT) 

JX  -  NCTRXL (K, IT) 

3820  IMODX  -  JX  +  10000 

475  IMODY  -  JY  +  10000 

IMODXY (M, IT)  -  IM0DX4100000  +  IMODY 
IF  (  JX  .EQ.  9999  )  GO  TO  3900 
X  -  JX 
Y  »  JY 

480  IF  (  WCELL  .EQ.  4.0  )  GO  TO  3825 

X  -  X  /  25.4 

_  Y  -  Y  /  25.4 

3825  YTEMPL (IT)  -  AHIN1 (YTEMPL ( IT), Y) 

YTEMPH  -  AMAX1 ( YTEMPH, Y) 

485  XTEMPL ( IT)  *  AMIN1 ( XTEMPL ( I T), X) 

XTEMPH(IT)  -  AMAX1(XTEMPH(IT),X) 

M  -  M  ♦  1 

IF  (  M  .GT.  3000  )  GO  TO  3845 
IF  (  JX  .EQ.  NCTRXR (K,IT)  )  GO  TO  3835 
490  JX  *  JX  +  JCELL 

GO  TO  3820 
3835  K  »  K  +  1 

GO  TO  3810  , 

384?  WRITE  (  6,3610  ) 

495  GO  TO  9900 

_ 3900  CONTINUE 

MCELLS  -  M  -  1 
WRITE  (  6,1020  ) 

WRITE  (  6,1122  )  MCELLS 

TOO  NXTEMP ( IT )  *  ( XTEMPH ( IT) -XTEMPL (IT) ) /W  +  1.001 

NYTEMP (IT)  -  ( YTEMPH-YTEMPL ( IT) )/W  ♦  1.001 
IF  (  IT  .EQ.  4  )  GO  TO  4000 
IT  -  IT  +  1 

IF  C  KTRGTC ( IT )  .EQ.  0  )  GO  TO  4000 
505  GO  TO  3805 

4000  CONTINUE 

IF  (  NSTCRD  .GT.  0  )  GO  TO  4200 
4010  NANGLE  -  1 

4TZ3  NALPHA  -  NANGLE  /  2 

510  NBETA  -  NALPHA  *  2 

-  IF  (  NBETA  .EQ.  NANGLE  )  GO  TO  4105 

IC  -  NALPHA  +  1 

E  IF  (  ISTMOV  .NE.  1  )  GO  TO  4025 
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489-  403 

MCTRXR(K*IT)  IS  X  COORDINATE  OF  CENTER  OF  CELL  AT  RIGHT  EDGE  OF  STRIP.  WHEN 

JX  EQUALS  THIS  COORDINATE*  STRIP  HAS  BEEN  COMPLETELY  PROCESSED.  IN  THAT 

CASE*  RESET  K  To  CORRESPOND  TO  NEXT  STRIP  OF  TARGET  CELLS  AND  CONTINUE  AT 
PROGRAM  STATEMENT  3810.  IN  OTHER  INSTANCES*  RESET  JX  TO  EQUAL  X  COORDINATE 

OF  NEXT  CELL  IN  STRIP  AMD  BEGIN  PROCESSING  FOR  THIS  CELL  AT  PROGRAM 

STATEMENT  3820. 

494-  495 

PRINT  ERROR  INDICATION  AND  GO  TO  PROGRAM  STATEMENT  Q900  TO  STOP  RUN, 

49-_  499 

MCELLS  IS  TOTAL  NUMBER  OF  TARGET  CELLS  FOR  ANGLE  ASSOCIATED  WITH  IT.  CELL 
COUNTER  M  IS  REDUCED  BY  1  BECAUSE  END  INDICATOR  HAS  BEEN  INCLUDED.  WRITE 
MCELLS  FOR  POSSIBLE  USE  IN  CHECKING  RUN  OUTPUT. 

500-  501 

CALCULATE  NUMBER  OF  CELLS  ALONG  EACH  SIDE  OF  RECTANGLE  THAT  IS  JUST  LARGE 

ENOUGH  TO  ENCLOSE  TARGET  FOR  ORIENTATION  CORRESPONDING  TO  IT.  EXPLANATIONS 
FOR  PROGRAM  LINES  456  AND  457  APPLY  HERE  ALSO  IF  IC  IS  REPLACED  BY  IT. 

502 

CONTINUE  AT  PROGRAM  STATEMENT  4000  IF  IT  VALUES  1  THROUGH  4  HAVE  ALL  BE Ef 
COVERED. 

503-  505 

RESET  it  FOR  NEXT  ORIENTATION  ANGLE.  KTRGTC (IT)  CAN  EQUAL  0  ONLY  IF  RUN 

INVOLVES  FEWER  THAN  FOUR  ANGLES.  IN  THAT  CASE*  CONTINUE  AT  PROGRAM  STATEMENT 
4000.  OTHERWISE  BEGIN  PROCESSING  FOR  NEW  ORIENTATION  AT  PROGRAM  STATEMENT 
3805. 

r 

0 

in 

SKIP  PROGRAM  LINES  508  THROUGH  532  IF  RUN  DOES  NOT  INVOLVE  VULNERABILITY  DATA 
FROM  TAPE  OR  DISC. 

508 

INITIAL  SETTING. 

509-  510 

NANGLE  EQUALS  1*  3*  5*  OR  7  FOR  TARGET  ORIENTATION  ANGLE  OF  0*  30*  60*  OF:  90 
DEGREES  RESPECTIVELY  AND  2*  4*  6*  OR  8  FOR  CORRESPONDING  REVERSE  ANGLE, 
MALPHA  EQUALS  0  WHEN  MANGLE  IS  1  OR  1  WHEN  NANGLE  IS  2.  CORRESPONDING  VALUES 
OF  NBETA  ARE  0  OR  2  RESPECTIVELY  FOR  MANGLE  EQUALLING  1  OR  2.  NBETA  ALWAYS 
DIFFERS  FROM  NANGLE  FOR  FORWARD  ORIENTATION  ANGLES  AMD  EQUALS  NANGLE  FOR 
REVERSE  ANGLES. 

511 

SKIP  PROGRAM  LIMES  512  THROUGH  542  WHEN  NANGLE  CORRESPONDS  TO  REVERSE  ANGLE. 

512 

IC  IS  SET  TO  1*  2*  3*  OR  4  RESPECTIVELY  FOR  NANGLE  EQUALLING  1*  3*  5*  OR  7. 

513 

PROGRAM  LINES  514  THROUGH  523  APPLV  ONLY  FOR  MOVING  TARGET  OP.  MOVING  FIRING 
WEAPON  RUN  AMD  ARE  SKIPPED  OTHERWISE. 

514 

IF  DUMMY  SET  FOR  WHICH  HO  CALCULATIONS  ARE  NEEDED  IS  INVOLVED*  SKIP  PROGRAM 
LIMES  515  THROUGH  1374. 

515-  523 

INITIAL  SETTINGS. 

524—  528 

IGNORE  STATEMENTS  APPLICABLE  ONLY  FOP.  NADJST  EQUALLING  2  THROUGH  5. 

530-  531 

INITIAL  SETTINGS. 

532 

SKIP  PROGRAM  LINES  533  THROUGH  537  APPLICABLE  ONLY  WHEN  F.UN  INVOLVES  TARGET 
SHAPE  DATA  FROM  CARDS. 

533-  537 

SPECIAL  SETTINGS  ARE  NEEDED  WHEN  TARGET  SHAPE  DATA  FROM  CARDS  ARE  USED. 

534 

INITIAL  SETTING. 

535-  537 

SET  IC*  NANGLE*  AND  NBETA  SO  THAT  PROGRAM  STATEMENT  4100  AND  SUBSEQUENT 
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PROGRAM  LISTING  (CONTINUED) 


- IF  (  XXB(IC,MSET)  .EQ.  999.9999  )  GO  TO  9305 

315  BBX(l)  -  XXB ( ICjMSET)  *  CONVRT 

BBY(l)  •  YYB<IC,MSET)  *  CONVRT 
SSIGX(l)  -  SSIGXBt IC#  MSET)  *  CONVRT 
SSIGY(l)  ■  SSIGYBt IC/MSET)  *  CONVRT 
DO  4040  N  *  2*4 

520  BBX(N)  -  0.0 

BBY(N)  ■  0.0 
SSIGX(N)  -  0.0 
_ 4040  SSIGYtN )  •  0.0 


4025 

IF 

( 

NADJST  .LT.  2  .OR.  NADJST 

.GT.  5 

)  GO  TO  4045 

525 

IF 

( 

IC  .EQ.  2  .OR.  IC  .EQ. 

3  ) 

NANGLE 

-  NANGLE  +  1 

IF 

( 

IC  .EQ.  2  .OR.  IC  .EQ. 

3  ) 

GO  TO 

9100 

IF 

( 

IC  .EQ.  4  .AND.  YBASE 

.GT. 

0.0  ) 

NANGLE  -  NANGLE  ♦ 

IF 

( 

IC  .EQ.  4  .AND.  YBASE 

.GT. 

0.0  ) 

GO  TO  9100 

4045 

CONTINUE 

536 

IB 

m 

IC  *  2 

IA 

• 

IB  -  1 

-  GO  TO  4100 

5  CONTINUE 

- NTRGTC  -  1 

T35  42l0  IC  -  NTRGTC 

NANGLE  •  0 
NBETA  •  9 

4lOO  NX  -  NXTEMP(IC) 
NY  -  NYTEMP  < IC ) 
540  CNX  -  NX 


CNY  «  NY 

XCORNR  -  XTEMPL  < IC ) 

- 4T(J5  IF  (  NBETA  .EQ.  NANGLE  )  XCORNR  -  -XTEMPHIIC) 

”  DO  4110  I  -  1# NX 
545  U  -  I  -  1 

4110  AX ( I )  -  XCORNR  +  U*W 
DO  4120  J  •  1#  NY 

Y  *  J  -  1 

4120  AY( J )  -  YTEMPL < IC )  +  V*W 
550  DO  4130  I  -  1,  NX 

DO  4130  J  -  1,NY 
4130  XK( J,I)  ■  0 

-  DO  4140  M  -  1,3000 

IMODX  -  IMODXYt  M,IC ) /100000 
IMODY  -  IMODXY ( M, IC )  -  IM0DX+100000 
X  -  IMODX  -  10000 
I  END  -  X 

IF  (  IEND  .EQ.  9999  )  GO  TO  4145 
IF  (  NBETA  .EQ.  NANGLE  )  X  «  -X 

Y  -  IMODY  -  10000 
IF  (  WCELL  .EQ.  4.0  )  GO  TO  4155 
X  -  X  /  25.4 

Y  -  Y  /  25.4 

4155  I  ■  (X-XCORNR)/ W  +  1.001 

J  -  (Y-YTEMPL (IC) ) /W  ♦  1.001 
4140  IK( J, I )  -  M 
TT75  CONTINUE 

-  IF  (  NHTKLL  .EQ.  9  )  GOTO  7400 

IF  (  NSTCRD  .GT.  0  )  GO  TO  4300 
570 _  IF  (  INTPL.  EQ.  0  .OR.  ISPLIT  .EQ.  0  )  GO  TO  4300 


555 


560 


565 
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538-  541 

542 

543 

544-  549 

550-  552 
553-  566 

568 

560-  5->0 

571-  5-4 

5*,6 

577-  578 
580 

581-  585 

587 

588 

589-  597 


INSTRUCTIONS  APPLY  FOR  TARGET  SHAPE  CARD  INPUT  AS  WELL  AS  FOR  VULNERABILITY 
DATA  INPUT  FROM  TAPE  OR  DISC.  SELECTION  OF  ARBITRARY  BUT  DIFFERENT  VALUES 
FOR  NANGLE  AND  NBETA  ENSURES  CALCULATIONS  .NORMALLY  ASSOCIATED  WITH  REVERSE 
ANGLES  ARE  NOT  ATTEMPTED. 

NX  AND  MY  ARE  SET  TO  EQUAL  NUMBERS  OF  CELLS  ON  EDGES  OF  SMALLEST  RECTANGLE 
ENCLOSING  TARGET  FOR  IC  VALUE  OF  CURRENT  CONCERN.  CNX  AND  CNY  ARE  REAL 
FORMS. 

XCORNR  IS  SMALLEST  X  COORDINATE  CORRESPONDING  TO  CENTER  OF  AT  LEAST  ONE  TARGET 
CELL.  FOR  EXAMPLE,  IF  PREVIOUS  CALCULATIONS  HAVE  RESULTED  IN  XTEMPL(IC) 

BEING  -204  AMD  XTEMPH(IC)  128  FOR  IC  VALUE  OF  INTEREST,  XCORNR  EQUALS  -204. 

XCORNR  NEEDS  TO  BE  RESET  FOR  REVERSE  ANGLE.  FOR  EXAMPLE,  IF  XTEMPL(IC)  EQUALS 
-204  AND  XTEMPH (IC )  128,  CORRECT  SETTING  FOR  XCORNR  IS  NOW  -128. 

AX  AND  AY  ARRAYS  MEED  TO  CONTAIN  X  AND  Y  COORDINATES  RESPECTIVELY  OF  CENTERS 
OF  CELLS  IN  SMALLEST  RECTANGLE  ENCLOSING  TARGET.  ARRAYS  MAY  POSSIBLY  INCLUDE 
SOME  VALUES  FOR  WHICH  THERE  IS  NO  CORRESPONDING  TAPGET  CELL. 

INITIAL  SETTINGS. 

COORDINATES  OF  CENTER  OF  EACH  TARGET  CELL  ARE  RETRIEVED  AND  CONVERTED  TO  PAIR 
OF  INDICES  J  AND  I  THAT  IDENTIFY  CELL  LOCATION  IN  RECTANGLE  ENCLOSING 
TARGET.  ONCE  PROCESSING  IS  DONE  FOR  ALL  CELLS,  IK  ARRAY  CONTAINS  1)  ZEROS 
FOR  VALUES  OF  J  AND  I  ASSOCIATED  WITH  CELLS  OUTSIDE  TARGET,  AND  2)  CELL 
COUNTER  INDEX  M  FOR  J  AND  I  PAIRS  CORRESPONDING  TO  TARGET  CELLS.  USE  OF 
INDICES  J  AND  I  IN  IK  ARRAY  PARALLELS  EARLIER  USE  OF  THESE  SAME  INDICES  IN 
AY  AND  AX  ARRAYS  RESPECTIVELY. 

SKIP  PROGRAM  LINES  569  THROUGH  1132  IF  MHTKLL  INDICATES  RUN  DOES  NOT  INVOLVE 
SIMULATED  FIRING  ENGAGEMENTS. 

SKIP  PROGRAM  LINES  571  THROUGH  574  UNLESS  ISPLIT  EQUALLING  1,  2,  OR  3 
INDICATES  FIRST,  SECOND,  OR  THIRD  OF  FOUR  SETS  OF  SIMULATED  ENGAGEMENTS 
INVOLVED  IN  RANGE  INTERPOLATION  HAVE  BEEN  COMPLETED. 

ISPLIT  EQUALLING  1  INDICATES  FIRST  SET  OF  SIMULATED  ENGAGEMENTS  INVOLVED  IN 
RANGE  INTERPOLATION  HAS  BEEN  COMPLETED.  FIRST  SET  IS  FOR  FORWARD  ORIENTATION 
ANGLE  OF  CONCERN.  RESETTING  OF  ISAMP  TO  0  BY  PROGRAM  STATEMENT  5000  IS 
NECESSARY  BEFORE  SIMULATION  OF  ENGAGEMENTS  FOR  REVERSE  ANGLE  CAN  BEGIN.  IF 
ISPLIT  EQUALS  2  OR  3,  NEMD1  ENGAGEMENTS  HAVE  ALREADY  BEEN  SIMULATED  FOR  BOTH 
FORWARD  AND  REVERSE  ANGLES.  IF  NECESSARY,  RESET  ISAMP  TO  NEND1  BEFORE  NEXT 
SET  OF  ENGAGEMENTS  IS  INITIATED. 

BASIC  ADJUSTMENT  PROCEDURE  DOES  NOT  INVOLVE  PROGRAM  LINES  577  AND  578. 

INITIAL  SETTINGS. 

MSMDTR  EQUALLING  0  INDICATES  PROGRAM  LIMES  581  THROUGH  585  ARE  TO  BE  SKIPPED. 

INITIAL  SETTINGS. 

PROGRAM  LINES  588  THROUGH  643  ARE  OF  POSSIBLE  CONCERN  FOP.  FORWARD  ORIENTATION 
ANGLE  OF  TARGET  BUT  NOT  FOR  REVERSE  ANGLE. 

PROGRAM  LINES  589  THROUGH  597  ARE  NEVER  INVOLVED  AFTER  PROCESSING  FOR  FIRST 
PAIR  OF  TARGET  ORIENTATION  ANGLES  HAS  BEEN  COMPLETED. 

ENGAGEMENT  IS  TO  BE  REPEATEDLY  SIMULATED  TOTAL  OF  NS  AMP  TIMES  FOR  EACH  SET  OF 
CONDITIONS.  EACH  SIMULATED  ENGAGEMENT  IS  CALLED  SAMPLE  ENGAGEMENT  OR  SIMPLY 
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575 


4533 


4310 

- 4135 

$56 


585  4320 

- 4*0(5 


590 

4410 

595 

4415 

4425 

4500 

600 

605 


610 


4510 


615 


620 
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4520 

4530 

4540 

4550 

4560 


IF  (  ISPLIT  .EQ.  1  )  GO  TO  5000 

IF  {  ISPLIT  .EQ.  2  )  GO  TO  5100 

ISAHP  -  NEND1 
GO  TO  5100 
CONTINUE 

IF  {  NADJST  .EQ.  0  )  GO  TO  4305 
DO  4310  I  -  1,NRDS 
AFAIL ( I )  »  0.0 
CONTINUE 

IF  {  NSHDTR  .EQ.  0  )  GO  TO  4400 

DO  4320  N  -  1*5 

DO  4320  I  -  1,2 

DO  4320  J  -  1,50 

NHINUS ( J, I, N )  -  0 

NPLUS { J,I,N )  -  0 

CONTINUE 

IF  {  NBETA  .EQ.  NANGLE  )  GO  TO  5000 

IF  C  NANGLE  .GT.  1  )  GO  TO  4500 

NSAHP  -  10000 

SAMP  -  NSAHP 

WRITE  {  6,4410}  NSAHP 

FORHAT  {  /  /,  10X, 21H  NUHBER  OF  SAHPLES  -  ,16  ) 

IF  {  NSTCRD  .GT.  0  )  GO  TO  4415 
IF  (  INTPL  .EQ.  1  )  GO  TO  4425 
NEND1  -  NSAHP 
GO  TO  4500 

IF  (  NANGLE  .EQ.  1  )  NEND1  -  RATIO  *  SAHP 
CONTINUE 

IF  {  NSTCRD  .GT.  0  )  GO  TO  4600 

IF  (  NRANGE  .GT.  1  }  GO  TO  4600 

IF  {  NANGLE  .GT.  1  )  GO  TO  4600 

DROP  -  0.0 

IF  {  NDROP  .EQ.  0  }  GO  TO  4600 
IF  {  NDROP  .GT.  1  }  GO  TO  4510 
DROP  -  (YTOP-YBASE)  *  0.5 
GO  TO  4600 

IF  (  NADJST  ,NE.  2  .AND,  NADJST  .NE.  4  )  GO  TO  4600 
PI  «  3.14159 
R  -  IRANGE 

CNHTER  -  R*PI  /  3200.0 
CNINCH  -  CNHTER  *  39.37 
IF  (  NDROP  .GT.  2  )  GO  TO  4520 
DROP  -  0.5  *  CNINCH 
GO  TO  4600 

IF  {  NDROP  .GT.  3  )  GO  TO  4530 
DROP  »  1.0  *  CNINCH 
GO  TO  4600 

IF  {  NDROP  .GT.  4  )  GO  TO  4540 
DROP  «  2.0  ♦  CNINCH 
GO  TO  4600 

IF  {  NDROP  .GT.  5  )  GO  TO  4550 
DFLCTN  «  1.0  *  CNINCH 
GO  TO  4600 

IF  {  NDROP  .GT.  6  )  GO  TO  4560 
DFLCTN  -  2.5  *  CNINCH 
GO  TO  4600 

DFLCTN  -  5.0  *  CNINCH 
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SAMPLE  IN  THESE  EXPLANATIONS.  SAMP  IS  REAL  FORM.  NEND1  REPRESENTS  NUMBER  OF 
ENGAGEMENTS  THAT  REQUIRE  CONSIDERATION  AS  PART  OF  FIRST,  AND  POSSIBLY  ONLY, 
GROUP  PROCESSED.  NEND1  IS  SET  EITHER  TO  NSAMP  OR,  IF  RANGE  INTERPOLATION  OF 
TARGET  VULNERABILITY  DATA  IS  INVOLVED,  TO  APPROPRIATE  FRACTION  OF  NSAMP .  IN 
LATTER  CASE  PROGRAM  STATEMENT  6015  SUBSEQUENTLY  CAUSES  REMAINING  PORTION  OF 
NSAMP  ENGAGEMENTS  TO  BE  ACCOUNTED  FOR. 

599-  601 

PROGRAM  LINES  602  THROUGH  627  CAM  BE  ,0F  CONCFRN  ONLY  WHEN  FIRST  RANGE  A  -ID 

FIRST  PAIR  OF  TARGET  ORIENTATION  ANGLES  ARE  BEING  PROCESSED.  CONDITION 
INVOLVING  NSTCRD  IS  INCLUDED  BECAUSE  SPECIAL  SETTINGS  OF  NANGLE  AND  N3ETA 
ESTABLISHED  BY  PROGRAM  LINES  536  AND  537  WHEN  RUN  INVOLVES  TARGET  SHA 3E  DATA 
FROM  CARDS  WOULD  OTHERWISE  CAUSE  PROBLEM. 

602 

INITIAL  SETTING. 

603 

SKIP  PROGRAM  LINES  604  THROUGH  627  IF  NDROP  IS  0. 

604-  627 

RESET  STANDARD  VERTICAL  ADJUSTMENT  DROP  TO  HALF  OF  TARGET  HEIGHT  FOR  NDROP 
EQUALLING  1.  NDROP  VALUES  2  THROUGH  6  ARE  ASSOCIATED  WITH  NADJST  EQUALLING 

2  THROUGH  5.  IGNORE  CORRESPONDING  PROGRAM  STATEMENTS. 

629-  634 

INITIAL  SETTINGS. 

635-  643 

CONDITIONS  OF  INTEREST  INCLUDE  PAIP  OF  TARGET  ORIENTATION  ANGLES  REPRESENTED 

BY  MANGLE  AND  NP.VRSE  OR  PARTICULAR  TARGET  ASSOCIATED  WITH  NTRGTC. 

644 

INITIAL  SETTING  BEFORE  ANV  SAMPLE  ENGAGEMENTS  FOR  PARTICULAR  CONDITIONS. 

645-  892 

CYCLE  FOR  EACH  SAMPLE  ENGAGEMENT  SIMULATED. 

645-  653 

INITIAL  SETTINGS  FOR  PARTICULAR  SAMPLE  ENGAGEMENT  BEING  SIMULATED. 

654—  656 

R1  AND  R2  ARE  RANDOM  NORMAL  DEVIATES  DRAWN  FROM  NORMAL  DISTRIBUTION  WITH  MEAN 
OF  0  AMD  STANDARD  DEVIATION  OF  1.  CALCULATED  VALUES  OF  XC1  AND  YC1  ARE 
COORDINATES  OF  POINT  IN  TARGET  COORDINATE  SYSTEM  WHERE  FIRST  ROUND  WOULD 
IMPACT  IF  THERE  WERE  MO  RANDOM  ERRORS.  XC1  IS  SUM  OF  1)  HORIZONTAL 

COORDINATE  XC  OF  INTENDED  AIMPOIMT,  2)  HORIZONTAL  FIXED  BIAS  XB,  AND  3) 
HORIZONTAL  VARIABLE  BIAS  DETERMINED  FOR  PARTICULAR  FIRING  ENGAGEMENT  OF 
CONCERN  BY  APPLYING  R1  AS  ADJUSTMENT  FACTOR  TO  VARIABLE  BIAS  STANDARD 
DEVIATION  SIGXB .  SIMILARLY  FOR  YC1. 

657-  658 

INITIAL  SETTINGS. 

65q-  891 

CYCLE  FOR  EACH  ROUND  FIRED. 

659-  660 

INITIAL  SETTINGS. 

661-  669 

COORDINATES  XI  AND  Y1  OF  ROUND  IMPACT  POINT  IN  PLANE  OF  TARGET  ARE  DETEF MIMED . 
FAC  IS  CONTROL  THAT  CAN  EQUAL  ONLY  0  OR  1.  FAC  IS  INITIALLY  SET  TO  1  TO 
INDICATE  THAT  NO  ADJUSTMENT  BASED  ON  SENSING  HAS  YET  OCCURRED,  AND  RESET  TO 

0  AFTER  FIRST  SUCH  ADJUSTMENT.  LAY  ERRORS  ARE  EXCLUDED  FROM  FURTHER 
CONSIDERATION  WHEN  ADJUSTMENT  BASED  ON  SENSING  HAS  BEEN  MADE,  BECAUSE 
ASSOCIATED  SENSING  ERRORS  ARE  DEFINED  TO  INCLUDE  LAY  ERRORS  AS  WELL  AS  ROUND 
OBSERVATION  ERRORS.  HORIZONTAL  LAY  ERROR  X2,  IF  APPLICABLE,  IS  GOTTEN  BY 
APPLYING  R1  TO  HORIZONTAL  LAY^  ERROR  STANDARD  DEVIATION  SIGXL.  SIMILAPLY  FOR 
VERTICAL  LAY  ERROR  Y2.  HORIZONTAL  RANDOM  ERROR  XRE  ACCOUNTS  FOP  LAY  ERROR, 

IF  NECESSARY,  AMD  ALWAYS  INCLUDES  HORIZONTAL  ROUND-TO-ROUND  ERROR 

DETERMINED,  FOR  SPECIFIC  ROUND  BEING  CONSIDERED,  BY  APPLYING  R1  TO 

HORIZONTAL  ROUND-TO-ROUND  STANDARD  DEVIATION  SIGXR.  SIMILARLY  FOR  VERTICAL 
RANDOM  ERROR  YP.E.  COORDINATES  XI  AND  Y1  OF  ACTUAL  IMPACT  POINT  ARE  OBTAINED 
BY  ADJUSTING  XC1  AMD  YC1  TO  REFLECT  CONTRIBUTION  OF  RANDOM  ERRORS. 
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PROGRAM  LISTING  (CONTINUED) 


630 

4610 

4620 

635 

4630 

640 

4640 

4650 

5000 

TS5 - 5T5T5 


650 

5110 

655 

660 

5200 

665 

5205 

670 

675 

5215 

5220 

660 

CONTINUE 

DO  4610  J  •  1,10 
DO  4610  I  »  1,61 
AKIM  I,  J  )  •  0.0 
DO  4620  J  -■  1,12 
DO  4620  I  •  1,NRDS 
Z(I,J)  -  0.0 

IF  {  NSTCRD  .GT.  0  )  GO  TO  4640 
NRVRSE  •  NANGLE  ♦  1 
WRITE  (  6,4630  )  NANGLE, NRVRSE 
FORMAT  (  /  /,  10X, 10H  ANGLES  -  ,217  ) 

WRITE  (  6,1020  ) 

GO  TO  5000 

WRITE  (  6,4650  )  NTRGTC 

FORMAT  C  /  /,  1  OX, 10H  TARGET  ■  ,  17  ) 

WRITE  (  6,1020  ) 

ISAMP  •  0 

ISAMP  •  ISAMP  +  1 

T  -  FLT 

SUMCASE  -  0.0 

XC1  -  0.0 

YC1  -  0.0 

IHW  -  0 

NF AIL  -  0 

DO  5110  J  -  1,10 

C(J)  -  0.0 

CALL  NRAN  (Rl, R2  ) 

XC1  -  XC1  +  XC  ♦  XB  +  R1*SIGXB 
YC1  -  YC1  +  YAIM  +  YB  ♦  R2*SIGYB 
FAC  -  1.0 
J  *  0 
J  -  J  +  1 
NSHORT  -  0 

IF  {  FAC  .EQ.  0.0  )  GO  TO  5205 
CALL  NRAN(R1,R2 ) 

X2  -  R1  *  SIGXL 
Y2  -  R2  *  SIGYL 
CALL  NRAN(R1,R2) 

XRE  «  FAC*X2  +  R1*SIGXR 
YRE  ■  FAC*Y2  +  R2*SIGYR 
XI  «  XC1  +  XRE 
Y1  -  YC1  ♦  YRE 
CALL  NRAN(T1,T2DUM) 

IF  (  J  .GT.  1  )  GO  TO  5215 
HWM  -  XM1  *  EXP (T1*STD1 ) 

IF  (  HWM  .LT.  AMT1  )  HWM  -  AMT1 
T  *  T  +  TF1  ♦  HWM 
GO  TO  5220 

HWM  ■  XM2  *  EXP (T1*STD2 ) 

IF  (  HWM  .LT.  AMT2  )  HWM  -  AMT2 

T  -  T  +  TFS  +  HWM 

CONTINUE 

J10  *  2.0  ♦  T/2.0 

IF  (  NADJST  .LT.  2  .OR.  NADJST  .GT.  5  )  GO  TO  5225 
IF  {  WCELL  .EQ.  100.0  .AND.  YBASE  .GT.  0.0  )  GO  TO  5235 
IF  {  XLEFTF(IC)  .LE»  XI  .AND.  XI  .LE.  XRGHTF { IC )  )  GO  TO  5245 
IF  (  NADJST  .EQ.  2  .OR.  NADJST  .EQ.  3  )  GO  TO  5665 


65 


EXPLANATIONS  (CONTINUED) 


6?0-  678 

cm  tap  cm  put  TTMF  T  t  S  RESET  TO  EQUAL  TIME  AT  WHICH  ROUND  UNDER  CONSIDERATION 
REACHES  TARGET  RANGE.  SECOND  RANDOM  NORMAL  DEVIATE  T2DUM  IS  UNNEEDED  AND 

TENOR  ED  DISTINCTION  IS  MADE  BETWEEN  FIRST  ROUND  AND  SUBSEQUENT  ROUND-. 
FIRSTEROUND  firing  time,  measured  from  beginning  of  engagement,  INCLUDE. 
nYFn  TTMF  PflMPflNENT  TFl  THAT  MAY  EQUAL  0«  VARIABLE  COMPONENT*  WHICH  IS 

IJrinG  TINE  IF  FIXED  TIME  EQUALS  0,  IS  GOTTEN  FROM  LOGARITHMICONORMAL 
DISTRIBUTION  WITH  MEDIAN  TIME  XMl  AND  VARIABILITY  FACTOR  ST1.  HWM  REPRESENTS 
SPECIFIC  VARIABLE  TIME  GENERATED  FROM  DISTRIBUTION  OF  TIMES  BY  APPLICATION 

OF  T1  IF  INITIAL  VALUE  OF  HWM  IS  EVER  LESS  THAN  MINIMUM  TIME  AMT1,  LATTER 
VALUE  OVERRIDES.  UPDATED  VALUE  OF  T  IS  TIME  AT  WHICH  FIRST  ROUND  HITS  OR 
MISSES  TARGET*  BECAUSE  FLIGHT  TIME  FLT  HAS  ALREADY  BEEN  ACC°U^E^ 

INITIAL  SETTING.  NOTE  THAT  PROGRAM  DOES  NOT  PROCESS  FLIGHT  TIME  AND  F  RING 
TIME  OF  FIRST  ROUND  IN  ORDER  MATCHING  CHRONOLOGICAL  SEQUENCE  OF  BAT^^ 

PVFNTS  STNCE  firing  of  round  obviously  precedes  its  flight  toward  target. 
FIrSI  tIhE  FOR  EACH  SUBSEQUENT  ^OUND  IS  ESTABLISHED  BASICALLY  LIKE  FIR 

Dnimn  ptptMP  ttMF  EACH  SUBSEQUENT  ROUND  FIRING  TIME  IS  CONSIDERED  MEASURED 
fgsr^g1!?  WHICH  PREVIOUS  ROUND  IS  FIRED.  SUBSEQUENT  ROUND  FIRING  TIME 
a rrm iwtc  pnp  FTYED  TIME  COMPONENT  THAT  MAY  AGAIN  BE  0«  VARIABLE  TIME 
cSpSnENT>  WHICH  HAY  BE  TOTAL  FIRING  TIME,  IS  GOTTEN  FROM  LOGARITHMICONORMAL 
DISTRIBUTION  WITH  MEDIAN  TIME  XM2  AND  VARIABILITY  FACTOR  ST2.  HWM  REPRESENTS 
SPECIFIC1 VARIABLE  TIME  BASED  ON  Tl.  IF  NECESSARY,  INITIAL  VALUE  OF  HWH  IS 
OVERRIDDEN  BY  MINIMUM  TIME  AMT2.  UPDATED  VALUE  OF  T  IS  TIME  AT  WHICH 
^iirs'coiifnt  ROUND  OF  INTEREST  HITS  0P.  MISSES  TARGET.  ANOTHER  OBSERVATION  IS 
uSfuS  suppose!  Sr  EXAMPLE,  THAT  FOURTH  ROUND  IS  BEING  PROCESSED.  VALUE  OF 

T  PRIOR  TO  UPDATING  IS  TIME  AT  WHICH  THIRD  ROUND  HITS  OR  MISSES  TARGET. 
STRICTLY  SPEAKING,  ONE  NEEDS  TO  SUBTRACT  PROJECTILE  FLIGHT  TIME  TO  GET  T.ME 

AT  uh t r H  THIRD  ROUND  IS  FIRED,  ADD  FOURTH  ROUND  FIRING  TIME,  AND  THEN  ADD 

TIHE  OF  FLIGHT  FOR  FOURTH  ROUND.  FIRST  AND  THIRD  OF  THESE  STEPS  CLEAR.Y 

CANCEL  EACH  OTHER  OUT  AND  ARE  CONSEQUENTLY  NOT  EXPLICITLY  REQUIRED. 

680 

ttuc  T  is  CONVERTED  TO  INDEX  J10  ASSOCIATED  WITH  AKIL  ARRAY.  J10  EQUALLING  1 

IS  RECEPVED  FOP  0  SECONDS,  P.ATHEP.  THAN  FOR  ANY  INTERVAL.  SUBSEQUENT  INTEGER 
VALUES  OF  J 10  CORRESPOND  TO  TIME  INTERVALS  OF  0  TO  2  SECONDS  (EXCLUDING  BOTH 
0  AND  2  SECONDS),  2  TO  4  SECONDS  (EXCLUDING  4  SECONDS),  ETC  .  .  ., 
RESPECTIVELY.  ASSUME,  FOR  EXAMPLE,  THAT  T  EQUALS  47.8  SECONDS.  THEN,  T/2 
EQUALS  23.9,  AND  23.9  +  2  EQUALS  25.9.  SINCE  J10  IS  INTEGER  QUANT  IT Y>  IT  IS 
SET  TO  25. 

681-  694 

IGNORE  THESE  SPECIAL  STATEMENTS  FOP.  NADJST  EQUALLING  2  THROUGH  5. 

696-  710 

DETERMINE  WHETHER  ROUND  HITS  TARGET. 

696-  701 

&NY  ROUND  FOR  WHICH  VERTICAL  COORDINATE  OF  IMPACT  POINT  IS  LESS  THAN  VERTICAL 
*^COORO IN ATE^YBASE  DEFINING  TARGET  BASE  HISSES  TARGET.  WHEN (NADJST  EOU.LS  1. 
MSHDRT  TS  SET  TO  1  TO  INDICATE  THAT  PARTICULAR  MISSING  ROUND  UNDER 
CONSIDERATION  HAS  BEEN  SHORT  OF  TARGET.  FIRE  ADJUSTMENT  JJ°CEDURE  AS-OCUTED 
WITH  MAD J ST  EQUALLING  0  DOES  NOT  INVOLVE  DISTINGUISHING  WHETHER  MISSING 

ROUND  IS  SHORT  OR  NOT. 

7  03 

IGNORE. 

704 

HORIZONTAL  COORDINATE  OF  IMPACT  POINT  IS  SMALLER  THAN  SMALLEST  X  COORD .NATE 
ASSOCIATED  WITH  VERTICAL  RECTANGLE  ENCLOSING  TARGET.  ROUND  MISSFS  TO  LEFT. 

705 

ROUND  MISSES  TO  RIGHT. 

706 

ROUND  IS  TOO  LOW  TO  HIT  TARGET.  ALTHOUGH  COMPARISON  INVOLVING  YBASE  HAS  BEEN 
MADE  EARLIER,  IT  IS  POSSIBLE  YBASE  AND  AY(D-H  DIFFER  SLIGHTLY  DUE  TO 
ROUNDING. 
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PROGRAM  LISTING  (CONTINUED) 


4 


685 


690 


695 


700 


755 


710 


715 


720 


725 


730 


775 


7*0 


5235 


5777 


5255 


52*5 


5280 


IF 

( 

NADJST 

IF 

( 

XI  .LT. 

IF 

( 

XI  .GT. 

GO 

TC 

1  5800 

IF 

( 

XLEFTD  ( 

IF 

( 

NADJST 

IF 

( 

NADJST 

IF 

( 

XI  .LT. 

IF 

( 

XI  .GT. 

GO 

TC 

1  5800 

CONTINUE 

IF 

( 

NADJST 

IF 

( 

Yl  .LT. 

GO 

TO 

1  52*5 

IF 

( 

Yl  .GE. 

.EQ.  5  )  GO 
XLEFTF(IC) 
XRGHTF(IC) 

IC)  .LE.  XI 
.EQ.  2  .OR. 
.EQ.  5  )  GO 
XLEFTD(IC) 
XRGHTD(IC) 


.EQ.  1  )  GO 
YBASE  )  GO 

YBASE  )  GO 


TO  5625 
)  XC1  -  XC1  ♦ 

)  XC1  -  XC1  - 

.AND.  XI  .LE. 
NADJST  .EQ.  3 
TO  5625 
)  XC1  «  XC1  + 

)  XC1  «  XC1  - 


TO  5255 
TO  5600 

TO  52*5 


DFLCTN 

DFLCTN 

XRGHTD(IC)  )  GO  TO  52*5 
)  GO  TO  5665 

DFLCTN 

DFLCTN 


NSHORT  •  1 
GO  TO  5600 
CONTINUE 

IF  I  NADJST  .EQ.  *  .OR.  NADJST  .EQ.  5  )  GO  TO  5*00 

IF  I  XI  .LT.  AX(1)-H  )  GO  TO  5600 

IF  {  XI  ,GT,  AX (NX) +H  )  GO  TO  5600 

IF  (  Y1  .LT.  AY ( 1)— H  )  GO  TO  5600 

IF  {  Y1  .GT.  AY (NY ) +H  )  GO  TO  5600 

JV  -  MINK  { Yl— AY(1  )+H)/W  ♦  1.0*CNY) 

IH  -  MINK  (X1-AX(  1 ) +H) /W  +  1.0/CNX) 

IF  {  IK(JV*IH)  .LE.  0  )  GO  TO  5600 
IF  I  J  .GT.  1  )  GO  TO  5265 
IF  (  INTPL  .EQ.  1  )  GO  TO  5265 

IF  (  NDTRM1  .EQ.  0  .AND.  NDTRM2  .EQ.  0  )  GO  TO  5265 
COORD(l)  -  XI  -  XRE  -  XC 
C00RD(2)  -  Yl  -  YRE  -  YC 
DO  5270  I  ■  1*2 

IF  (  COORD(I)  .GE.  0.0  )  GO  TO  5275 
NCOORD  -  -COORD  (D/20.0  ♦  1.0 


IF  (  NCOORD  .GT.  50  )  NCOORD  -  50 
DO  5280  N  ■  1#  2 

IF  (  KDMSTC(N)  .EQ.  0  )  GO  TO  5280 

NMINUS  (NCOORD* Is N)  *  NMINUS(NCOORD*I*N>  +  1 
CONTINUE 


GO  TO  5270 

5275  NCOORD  -  COORD( I) /20.0  +  1.0 

IF  I  NCOORD  .GT.  50  )  NCOORD  -  50 
DO  5290  N  -  1*2 

IF  (  KDMSTC(N)  .EQ.  0  )  GO  TO  5290 
N PLUS (NCOORD# I#  N)  =  NPLUS (NCOORD* I*N)  +  1 
5290  CONTINUE 


5270  CONTINUE 
7775  CONTINUE 

_  IF  t  RANF(DUM)  .GT.  RELT  )  GO  TO  5700 

_  IF  (  NHIT  .EQ.  0  )  GO  TO  5*00 

CALL  NRAN(R1*R2) 

XC1  -  XC1  -  XI  +  XC  -  R1*H$X 
YC1  -  YC1  -  Yl  +  YC  -  R2*HSY 

_  FAC  -  0.0 

5*00  CONTINUE 

IF  (IHW  .GT.  0  )  GO  TO  5*05 
IHW  -  1 
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EXPLANATIONS  (CONTINUED) 


70? 

708-  710 


"’ll-  713 

714-  715 

71'6-  731 


733 

',3‘t 

735-  738 


ROUND  MISSES  HIGH. 

ROUND  HITS  WITHIN  RECTANGLE  ENCLOSING  TARGET.  IMPACT  POINT  IS  WITHIN  CELL 
WHOSE  CENTER  HAS  Y  AND  X  COORDINATES  EQUAL  TO  JV  AND  IH  RESPECTIVELY. 

FORMULAS  FOR  JV  AND  IH  SYSTEMATICALLY  RESOLVE  AMBIGUITIES  THAT  ARISE, 
EXCEPTIONALLY,  IF  IMPACT  POINT  IS  COMMON  TO  TWO  OR  FOUR  CELLS.  FOLLOWING 
ILLUSTRATION  SHOULD  CLARIFY  FORMULA  FOP.  JV.  1)  Y1-(AY(1)-H),  WHICH  EQUALS 
Yl-AY ( 1) +H,  IS  DISTANCE  IN  INCHES  FROM  IMPACT  POINT  TO  BOTTOM  EDGE  OF 
RECTANGLE.  LET  THIS  DISTANCE  EQUAL  22.415  INCHES.  2)  DIVISION  BY  CELL 
DIMENSION  W  CONVERTS  FROM  INCHES  TO  CORRESPONDING  NUMBER  OF  CELLS.  LET  W 
EQUAL  3.937  INCHES.  THEN  22.415  INCHES  EQUAL  5.693  CELLS.  3)  FOR  EXAMPLE 
CHOSEN,  IMPACT  POINT  OBVIOUSLY  LIES  IN  SIXTH  ROW  OF  CELLS.  SINCE  JV  IS 
INTEGER  QUANTITY,  ANY  DECIMAL  FRACTION  SUCH  AS  .693  IS  EVENTUALLY  DROPPED 
FROM  CALCULATION  AND  ADDITION  OF  1  TO  GET  6.693  IS  NECESSARY.  INTEGER  VALUE 
6  FOR  JV  IS  OBTAINED  BY  DROPPING  .693.  COMPARISON  OF  6.693  TO  C NY  ENSURES 
THAT  MAXIMUM  POSSIBLE  VALUE  NY  IS  NEVER  EXCEEDED  DUE  TO  ROUNDING.  SIMILAR 
EXPLANATIONS  APPLY  TO  IH.  PREVIOUS  PROCESSING  OF  TARGET  DESCRIPTION  AND 
VULNERABILITY  DATA  HAS  SET  IK  ARRAY  VALUES  TO  BE  POSITIVE  IF  AND  ONLY  IF 
CORRESPONDING  CELLS  ARE  INCLUDED  IN  TARGET.  ROUND  MISSES  IF  IT  IMPACTS  ON 
CELL  WHICH,  ALTHOUGH  IN  RECTANGLE  ENCLOSING  TARGET,  IS  NOT  PART  OF  TARGET. 

PROGRAM  LINES  714  THROUGH  731  INVOLVE  SPECIAL  CALCULATIONS  ASSOCIATED  WITH 
AT  LEAST  ONE  OF  CONTROLS  NDTRM1  AND  NDTRM2  EQUALLING  1.  THESE  CALCULATIONS 
CANNOT  BE  ATTEMPTED  IF  ROUND  INDEX  J  INDICATES  THAT  SOME  SUBSEQUENT  ROUND 
RATHER  THAN  FIRST  ROUND  IS  BEING  PROCESSED  OR  IF  RANGE  INTERPOLATION  OF 
TARGET  VULNERABILITY  DATA  IS  NEEDED. 

COORD (1)  AND  COORD( 2)  RESPECTIVELY  ARE  ESSENTIALLY  HORIZONTAL  AND  VERTICAL 
MISS  DISTANCES  IN  TARGET  PLANE  WITH  REFERENCE  TO  APPROXIMATE  CENTER  OF  MASS 
OF  TARGET.  TO  VERIFY  THIS,  SEE  HOW  XC1  AND  YC1  ARE  CALCULATED  BY  PROGRAM 
LINES  655  AND  656  AND  HOW  XI  AND  Y1  ARE  SET  BY  LINFS  668  AND  669.  COORD(l) 

IS  DIFFERENCE  BETWEEN  X  COORDINATE  OF  APPROXIMATE  CENTER  OF  MASS  AND  X 
COORDINATE  OF  POINT  THAT  WOULD  BE  HIT  REPEATEDLY  IF  THERE  WERE  NO  RANDOM 
ERRORS  AND  WEAPON  FIRED  SEVERAL  ROUNDS  WITHOUT  AMY  ADJUSTMENT  BETWEEN 
ROUNDS.  SIMILARLY,  C00RD(2)  IS  DIFFERENCE  OF  CORRESPONDING  Y  COORDINATES. 

INDEX  I  VALUES  1  AND  2  ARE  ASSOCIATED  WITH  PROCESSING  OF  HORIZONTAL  AND 

VERTICAL  MISS  DISTANCES  RESPECTIVELY.  PROGRAM  LINES  718  THROUGH  723  APPLY  TO 
NEGATIVE  COORD (I)  VALUES,  WHILE  LINES  725  THROUGH  730  COVER  POSITIVE  VALUES 
OR  0.  QUANTITIES  IN  NMINUS  AND  NPLUS  ARPAYS  ARE  COUNTERS  INDICATING  NUMBERS 
OF  ENGAGEMENTS  FOR  WHICH  MISS  DISTANCE  FALLS  WITHIN  ANY  20-INCH  INTERVAL. 
TOTAL  SPREAD  COVERED  IS  BASICALLY  FROM  -1000.0  TO  1000.0  (EXCLUDING  -1000.0 
AND  1000,0)  INCHES.  ANY  MISS  DISTANCE  OUTSIDE  THIS  SPREAD  WOULD  BE  COUNTED 
AS  BEING  WITHIN  INTERVAL  NUMBER  50  ASSOCIATED  WITH  NMINUS  OP  NPLUS  ARRAY.  N 
EQUALLING  1  CORRESPONDS  TO  ALL  FIRST  ROUNDS  THAT  HIT  TARGET.  CALCULATIONS 
FOR  N  EQUAL  TO  2  ARE  IDENTICAL  AT  THIS  STAGE  TO  THOSE  FOR  N  EQUAL  TO  1  BUT 
EVENTUALLY  APPLY  TO  FIRST  ROUND  HITS  THAT  DO  NOT  RESULT  IN  K  (COMPLETE 
DESTRUCTION)  KILL  OF  TARGET. 

RANDOM  NUMBER  RAMF(DUM)  IS  COMPARED  TO  PROBABILITY  P.ELT  THAT  ROUND  HAS  FLOWN 
RELIABLY  TO  TARGET  RANGE.  DEFECTS  CAUSING  ERRATIC  FLIGHT  OF  ROUND  AP.E 
CONSIDERED  TO  MAKE  IT  IMPOSSIBLE  FOR  ROUND  TO  HIT  AMD  FOR  FIRE  ADJUSTMENT 
TO  BE  ATTEMPTED.  PROCESSING  CONTINUES  AT  PROGRAM  STATEMENT  5700  AFTER 
DETERMINATION  OF  DEFECTIVE  TRAJECTORY. 

SKIP  PROGRAM  LINES  735  THROUGH  738  IF  MHIT  CONTROL  INDICATES  AIMPOINT  IS  NOT 
TO  BE  ADJUSTED  AFTER  TARGET  HIT. 

ADJUSTMENT  OF  AIMPOINT  AFTER  TARGET  HIT  IS  CONSIDERED  BASED  ON  GUNNER  SENSING 
OF  IMPACT.  UPDATED  VALUES  OF  XC1  AND  YCl  AP.E  COORDINATES  OF  POINT  IN  TARGET 
COORDINATE  SYSTEM  WHERE  NEXT  ROUND  WOULD  IMPACT  IF  THERE  WERE  NO  RANDOM 
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740-  741 

742-  748 


749 

’’SO 

751-  784 

751-  753 

754-  772 

773-  784 


EXPLANATIONS  (CONTINUED) 


ERRORS.  XC1  IS  OBTAINED  BY  COMBINING  1)  HORIZONTAL  COORDINATE  XC  OF  INTENDED 
AINPOINT*  2)  DIFFERENCE  XC1-X1  WHICH  INVOLVES  XC1  AND  XI  VALUES  FROM  EARLIER 
CALCULATIONS  AND  WHICH*  AS  CAN  BE  SEEN  FROM  PROGRAM  LINE  668*  CAUSES 
SPECIFIC  HORIZONTAL  RANDOM  ERRORS  THAT  AFFECTED  IMPACT  POINT  OF  HITTING 
ROUND  TO  BE  REMOVED*  AND  3)  ADJUSTMENT  ERROR  OBTAINED  BY  APPLYING  R1  TO 
HORIZONTAL  STANDARD  DEVIATION  HSX.  TERM  ACCOUNTING  FOR  ADJUSTMENT  ERROR 

involves  minus  sign  because*  if  hitting  round  is  judged  to  be  off  target 

CENTER  IN  ONE  DIRECTION*  WEAPON  IS  MOVED  IN  OPPOSITE  DIRECTION.  FORMULA  FOR 
YC1  OBVIOUSLY  PARALLELS  THAT  FOR  XC1.  FAC  IS  RESET  SO  THAT  LAY  ERRORS* 
CONSIDERED  INCLUDED  AS  PART  OF  ADJUSTMENT  ERRORS,  WILL  NO  LONGER  BE  TREATED 
EXPLICITLY. 

IF  IHW  EQUALS  0*  RESET  TO  1.  IHW  CAM  BE  ONLY  0  OR  1*  LATTER  INDICATING  THAT 
TARGET  HAS  BEEN  HIT  AT  LEAST  ONCE  IN  PARTICULAR  SAMPLE  ENGAGEMENT  UNDER 
CONSIDERATION.  ONCE  IHW  IS  SET  TO  1*  IT  REMAINS  UNCHANGED  FOR  REMAINDER  OF 
ENGAGEMENT. 

UPDATE  Z  AND  AKIL  ARRAYS  TO  ACCOUNT  FOR  TARGET  HIT  BEING  PROCESSED.  SET  INDEX 
KH  TO  5  FOR  0*  30*  60*  OP.  90  DEGREES  OR  TO  10  FOR  REVERSE  ANGLE  MATCHED  WITH 
EACH  FIRST  ANGLE.  EACH  Z(K*5)  OR  Z(K*10)  VALUE  REPRESENTS  HOW  MANY  TIMES 
FIRST  TARGET  HIT  OCCURRED  EITHER  ON  ROUND  K  OP.  ON  EARLIER  POUND  IN  SAMPLE 
ENGAGEMENTS  PROCESSED  SO  FAR.  INDEX  J10  ASSOCIATED  WITH  AKIL  ARRAY  IS 
RESTRICTED  AT  THIS  POINT  TO  MAXIMUM  VALUE  OF  61  TO  PRECLUDE  FURTHER 
CONSIDERATION  OF  HITS  THAT  OCCUR  AT  TIMES  EQUALLING  OR  EXCEEDING  LIMIT  OF 
120  SECONDS  (2  MINUTES).  AKIL  ARRAY  IS  UPDATED  TO  ACCOUNT  FOR  TIME  AT  WHICH 
HIT  UNDER  CONSIDERATION  OCCURS.  EACH  AKIL(K*5)  0P.  AKIL(K*10)  VALUE 
REPRESENTS  HOW  MANY  TIMES  TARGET  HAS  BEEN  HIT  AT  LEAST  ONCE  PRIOR  TO  TIME 
EQUALLING  2(K-1)  SECONDS  IN  SAMPLE  ENGAGEMENTS  PROCESSED  SO  FAR.  ANY 
PARTICULAR  ENGAGEMENT  DOES  NOT  CONTRIBUTE  TO  AKIL(K*5)  0P.  AKIL(K*10)  VALUE 
IF  TARGET  IS  HOT  HIT  AT  ALL  OR  IF  IT  IS  FIRST  HIT  AT  TIME  EQUALLING  OR 
EXCEEDING  2(K-1)  SECONDS. 

SKIP  PROGRAM  LINES  750  THROUGH  891  IF  NPRHIT  INDICATES  TARGET  HIT  ONLY  IS  OF 
CONCERN  BUT  KILLS  ARE  NOT.  ENGAGEMENT  IS  OVER  SINCE  FIRST  HIT  ON  TARGET  HAS 
BEEN  ACHIEVED. 

IT  IS  CONSIDERED  THAT  TARGET  CANNOT  BE  KILLED  IF*  FOR  ROUNDS  WITH  FUZE*  LATTER 
DOES  NOT  PERFORM  RELIABLY  WHEN  TARGET  IS  HIT.  INPUT  PROBABILITY  RELF  IS 
SIMPLY  SET  EQUAL  TO  1.0000  FOR  ANY  POUND  TYPE  THAT  DOES  NOT  INVOLVE  FUZING. 
PROCESSING  CONTINUES  AT  PROGRAM  STATEMENT  5700  AFTER  DETERMINATION  OF 
DEFECTIVE  FUZE  FUNCTIONING. 

KILLING  EFFECTS  THAT  RESULT  FROM  TARGET  BEING  HIT  BY  PARTICULAR  ROUND  OF 
CONCERN*  EXCEPT  FOR  PASSENGER  PERSONNEL  CASUALTIES  CONSIDERED  LATER,  ARE 
DETERMINED. 

INITIAL  SETTINGS.  HW2  IS  USED  FOR  TEMPORARY  STORAGE  OF  RANDOM  NUMBER  NEEDED- 
SUBSEQUENTLY  TO  DETERMINE  WHETHER  OR  NOT  TARGET  IS  KILLED  AS  RESULT  OF  BEING 
HIT.  K22  IS  USED  FOR  TEMPORARY  STORAGE  OF  IK(JV*IH)  VALUE  ASSOCIATED  WITH 
TARGET  CELL  HIT.  N  IS  INVOLVED  IN  UNPACKING  OF  ARRAY  MPK  QUANTITIES. 

UNPACK  MPK  ARRAY  QUANTITY  THAT  CORRESPONDS  TO  TARGET  CELL  HIT  AND  ASSOCIATED 
TARGET  ORIENTATION  TO  RETRIEVE  PROBABILITIES  OF  M  KILL*  F  KILL*  M  OR  F  KILL* 
AND  K  KILL  PER  HIT  AND*  IF  NECESSARY*  EXPECTED  CASUALTIES  PER  HIT.  STORE 
THESE  IN  PK  ARRAY*  USING  PK(1)  OR  PK(6)  FOR  M  KILL  PROBABILITY*  PK(2)  OR 
PK( 7)  FOR  F  KILL  PROBABILITY*  PK(3)  OP.  PK(8)  FOR  M  OR  F  KILL  PROBABILITY* 

PK( 4)  OR  PK ( 9)  FOR  K  KILL  PROBABILITY*  AND  PK(5)  OR  PK(10)  FOR  CASUALTIES. 

PK ( 1 }  THROUGH  PK(5)  ARE  USED  FOR  TARGET  ORIENTATION  ANGLE  OF  0*  30,  60*  OR 
90  DEGREES*  AMD  PK(6)  THROUGH  PK(10)  FOR  CORRESPONDING  REVERSE  ANGLE. 

CONSIDER.  KILL  CRITERIA  IN  TURN.  INDEX  K21  EQUALS  1  OR  6  FOR  M  KILL*  2  OR  7  FOR 
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PROGRAM  LISTING  (CONTINUED) 


I 


-  KH  »  5 

IF  (  NBETA  .EQ.  NANGLE  )  KH  »  10 
DO  5410  K  »  J,NRDS 

745  5410  Z (K,KH)  *  Z(K,KH)  +  1.0 

IF  {  J10  .GT.  61  )  GO  TO  5405 
DO  5420  K  »  J 10/61 
5420  AKIL<K,KH)  »  AKIL(K,KH)  +  1.0 
5405  IF  (  NPRHIT  .EQ.  1  )  GO  TO  6000 
750  IF  {  RANF(DUM)  .GT.  RELF  )  GO  TO  5700 

-  HW2  •  RAN  F { DUM ) 


755 

760 

5435 

5430 

765 

5445 

770 

5450 

775 

780 

5490 

5480 

5460 

785 


K22  »  IK(JV,IH) 

N  *  -49 

DO  5430  K  *  1,4 

IF  {  NBETA  .EQ.  NANGLE  )  GO  TO  5435 
1ST  -  SHIFT{MPK(K22,IA),N)  .AND.  MASK11 
ST  *  1ST 

PK  (K )  *  0.001  *  ST 
GO  TO  5430 

IST1  *  SHIFT(MPK(K22,IB),N)  .AND.  MASK11 
ST1  *  IST1 

PKtK+5)  *  0.001  *  ST1 
N  -  N  +  11 

IF  {  NBETA  .EQ.  NANGLE  )  GO  TO  5445 
1ST  *  MPKIK22, IA)  .AND.  MASK16 
ST  *  1ST 

PK(5)  *  0.001  *  ST 
GO  TO  5450 

IST1  *  MPK(K22, IB)  .AND.  MASK16 
ST1  -  I ST1 

PKtlO)  »  0.001  *  ST  1 
CONTINUE 

DO  5460  KILL  *  1,4 
K21  •  KILL 

IF  (  NBETA  .EQ.  NANGLE  )  K21  *  K21+5 

IF  (  C (K21)  .GT.  0.0  )  GO  TO  5460 

IF  (  PK (K21 )  .LT.  HW2  )  GO  TO  5460 

C (K21)  *  1.0 

DO  5490  K  *  J, NRDS 

Z (K,K21 )  *  Z (K,K21)  +  1.0 

IF  (  J10  .GT.  61  )  GO  TO  5460 

DO  5480  K  *  J10,61 

AKIL (K,K21 )  ■  AKILtK,K21)  +  1.0 

CONTINUE 

IF  (  J  .GT.  1  )  GO  TO  5505 
IF  (  INTPL  .EQ.  1  )  GO  TO  5505 


IF  {  NDTRM2  .EQ .  0  )  GO  TO  5505 
KKILL  =  4 

IF  (  NBETA  .EQ.  NANGLE  )  KKILL  *  KKILL  +  5 
790  IF  (  C { KKILL )  .EQ.  0.0  )  GO  TO  5505 

N  *  2 

COORD ( 1)  -  XI  -  XRE  -  XC 
COORDd)  *  Y1  -  YRE  -  YC 
DO  5510  I  ■  1,2 

795  IF  {  COORDd)  .GE.  0.0  )  GO  TO  5515 

NCOORD  =*  -COORDd) /20.0  +  1.0 
IF  (  NCOORD  .GT.  50  )  NCOORD  *  50 
NHINUStNCOORD, I,N)  «  NMINUS (NCOORD, I, N )  -  1 
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EXPLANATIONS  (CONTINUED) 


76  5- 

768- 

-91- 

805- 

805- 

807- 


F  KILL,  3  OR  8  FOR  M  OR  F  KILL,  AND  4  OR  9  FOR  K  KILL.  VALUES  1  THROUGH  4 
APPLY  TO  TARGET  ORIENTATION  ANGLES  OF  0,  30,  60,  AND  90  DEGREES,  WHILE  6 
THROUGH  9  APPLY  FOR  ASSOCIATED  REVERSE  ANGLES.  EACH  QUANTITY  IN  C  ARRAY  CAN 
BE  ONLY  0  OR  1,  LATTER  VALUE  INDICATING  TARGET  HAS  BEEN  KILLED  FOR  ANGLE 
UNDER  CONSIDERATION  AND  KILL  CRITERION  CORRESPONDING  TO  K21.  NO  FURTHER 
PROCESSING  IS  TO  BE  DONE  FOP.  PARTICULAR  ANGLE  AND  CRITERION  WHEN  RELATED  C 
VALUE  ALREADY  EQUALS  1.  IF  TARGET  KILL  HAS  NOT  ALREADY  BEEN  ACHIEVED  BY 
EARLIER  ROUND,  FOR  SPECIFIC  COMBINATION  OF  ANGLE  AND  CRITERION,  RANDOM 
NUMBER  HW2  AND  APPLICABLE  KILL  PROBABILITY  PK(K21)  ARE  COMPARED  TO  DETERMINE 
WHETHER  TARGET  IS  KILLED.  IF  TARGET  KILL  IS  ACHIEVED,  RESET  C(K21)  AND 
UPDATE  Z  AND  AKIL  ARRAYS  ACCORDINGLY.  Z  AMD  AKIL  ARRAY  UPDATING  IS  VERY 
SIMILAR  TO  THAT  ALREADY  EXPLAINED  FOP.  PROGRAM  LINES  742  THROUGH  748.  EACH 
Z(K,K21)  VALUE  REPRESENTS  HOW  MANY  TIMES  TARGET  ORIENTED  AS  INDICATED  BY  K21 
HAS  BEEN  KILLED  ACCORDING  TO  CRITERION  ASSOCIATED  WITH  K21  EITHER  ON  ROUND  K 
OP.  ON  EARLIER  ROUND  IN  SAMPLE  ENGAGEMENTS  PROCESSED  SO  FAR.  EACH  AKIL(K,K21) 
VALUE  REPRESENTS  HOW  MANY  TIMES  TARGET  HAS  BEEN  KILLED  PRIOR  TO  TIME 
EQUALLING  2(K-1)  SECONDS. 

787  PROGRAM  LINES  788  THROUGH  803  INVOLVE  SPECIAL  CALCULATIONS  ASSOCIATED  WITH 

NDTRM2  EQUALLING  1.  THESE  CALCULATIONS  CANNOT  BE  ATTEMPTED  IF  ROUND  INDEX  J 
INDICATES  THAT  SOME  SUBSEQUENT  ROUND  RATHER  THAN  FIRST  ROUND  IS  BEING 
PROCESSED  OR  IT  RANGE  INTERPOLATION  OF  TARGET  VULNERABILITY  DATA  IS  NEEDED. 

-790  SET  KKILL  TO  EQUAL  4  FOR  FORWARD  ANGLE  OR  9  FOP.  REVEPSE  ANGLE.  IF  C(KKILL) 
INDICATES  K  KILL  OF  TARGET  HAS  NOT  BEEN  ACHIEVED,  SKIP  PROGRAM  LINES  791 
THROUGH  803  BECAUSE  HORIZONTAL  AMD  VERTICAL  MISS  DISTANCES  HAVE  ALREADY  BEEN 
ACCOUNTED  FOP.  BY  LINES  714  THROUGH  731. 

803  CALCULATIONS  APPLY  ONLY  IF  K  KILL  HAS  BEEN  INFLICTED  ON  TARGET.  SINCE  MISS 

DISTANCES  ALREADY  PROCESSED  BY  LINES  714  THROUGH  731  ARE  NOT  TO  BE  COUNTED, 
CORRESPONDING  COUNTERS  IN  MMINUS  ARRAY  AMD/OP.  NPLUS  ARRAY  ARE  REDUCED 
ACCORDINGLY. 

•  828  PROCESSING  FOR  HITTING  ROUND  CONTINUES  WITH  ACCOUNTING  FOR  PASSENGER  PERSONNEL 
CASUALTIES,  AS  NECESSARY,  AND  DETERMINATION  OF  WHETHER  ENGAGEMENT  HAS  BEEN 
COMPLETED. 

806  INITIAL  SETTINGS. 


828 


NONZERO  VALUE  OF  PASSN  INDICATES  PASSENGER  PERSONNEL  CASUALTIES  NEED  TO  BE 
CONSIDERED.  EXPECTED  CASUALTIES  PER  HIT  HAVE  ALREADY  BEEN  STORED  IN  PK(5) 

FOP.  TARGET  ORIENTATION  ANGLE  OF  0,  30,  60,  OP.  90  DEGREES,  OR  IN  PK(10)  FOR 
ASSOCIATED  REVERSE  ANGLE.  RATIO  CASLT/PASSN  INDICATES  WHAT  FRACTION  OF 
PASSENGER  PERSONNEL  CARRIED  IN  TARGET  VEHICLE  WOULD  BECOME  CASUALTIES  FROM 
EFFECTS  OF  ROUND  BEING  PROCESSED  IF  MO  CASUALTY  HAD  BEEN  CREDITED  TO  ANY 
PREVIOUS  ROUND.  DIFFERENCE  PASSN-SUMCASE,  WHICH  EQUALS  ORIGINAL  NUMBER  OF 
PASSENGER  PERSONNEL  PASSN  IF  TOTAL  CASUALTIES  SUMCASE  INFLICTED  BY  PREVIOUS 
ROUNDS  EQUAL  ZERO,  REPRESENTS  HOW  MANY  PASSENGER  PERSONNEL  HAVE  NOT  ALREADY 
BECOME  CASUALTIES.  UPDATE  SUMCASE  TO  ACCOUNT  FOR  ADDITIONAL  CASUALTIES, 
REASONABLY  ESTIMATED  AS  PRODUCT  OF  RATIO  CASLT/PASSN  AND  DIFFERENCE 
PASSN-SUMCASE,  THAT  RESULT  FROM  ROUND  OF  CURRENT  CONCERN.  UPDATE  Z  ARRAY  SO 
THAT  Z(J>11),  OR  Z(J,12)  FOR  REVERSE  ANGLE,  REPRESENTS  TOTAL  NUMBER  OF 
PASSENGER  PERSONNEL  CASUALTIES  INFLICTED,  IN  ALL  SAMPLE  ENGAGEMENTS 
PROCESSED  SO  FAR,  BY  FIRST  J  ROUNDS.  FIRST  THREE  CONDITIONAL  STATEMENTS  IN 
DO  5550  LOOP  ENSURE  CONSIDERATION  OF  PROPER  K  VALUES,  BE  THEY  1  THROUGH  4  OP. 
6  THROUGH  9.  IF  C(K)  EVER  EQUALS  0,  ENGAGEMENT  IS  STILL  INCOMPLETE  UNLESS 
CURRENT  PROCESSING  INVOLVES  LAST  OF  NRDS  ROUNDS  ALLOWED  AS  MAXIMUM.  IF  ALL 
FOUR  C(K)  VALUES  OF  CONCERN  EQUAL  1,  PROCEED  TO  NEXT  ENGAGEMENT  BUT  ONLY 
AFTER  UPDATING,  AS  NECESSARY,  CUMULATIVE  PERSONNEL  CASUALTIES  IN  Z  ARRAY  FOP. 
ROUNDS  WHICH  NEED  NOT  BE  FIRED  IN  CURRENT  ENGAGEMENT. 
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PROGRAM  LISTING  (CONTINUED) 


800 

5515 

5510 

5505 

TOT 


810 

815 

5535 

5545 

820 

5550 

825 

5565 

5560 

5600 

FTO 


835 

840 

5615 

845 

5610 

5TO5 

850 

5635 

5625 

8  55 

GO  TO  5510 

NCOORD  -  COORD( I ) /20.0  +  1.0 

IF  (  NCOORD  .GT.  50  )  NCOORD  »  50 

NPLU5 (NCOORD# I,N)  -  NPLUS (NCOORD, I, N)  -  1 

CONTINUE 

CONTINUE 

KCAS  -  5 

IF  (  NBETA  .EQ.  NANGLE  )  KCAS  -  10 
IF  (  PASSN  .IE.  0.0  )  GO  TO  5545 
CASLT  «  PK.(KCAS) 

SUMCASE  -  SUMCASE  +  (CASLT/ PASSN) *( PASSN-SUMC AS E > 
IF  (  NBETA  .EQ.  NANGLE  )  GO  TO  5535 
Z(J,11)  -  Z(J,11)  ♦  SUMCASE 
GO  TO  5545 

Z ( J, 12 )  *  Z ( J, 12)  ♦  SUMCASE 
DO  5550  K  -  1,10 

IF  (  K  .EQ.  5  .OR.  K  .EQ.  10  )  GO  TO  5550 

IF  (  K  .GT.  KCAS  )  GO  TO  5550 

IF  (  KCAS  .EQ.  10  .AND.  K  .LT.  5  )  GO  TO  5550 

IF  (  C(K)  .EQ.  0.0  )  GO  TO  5800 

CONTINUE 

IF  (  J  .EQ.  NRDS  )  GO  TO  6000 

JNXT  -  J  +  1 

DO  5560  K  -  JNXT, NRDS 

IF  (  NBETA  .EQ.  NANGLE  )  GO  TO  5565 

Z(K,11)  -  Z ( K, 1 1 }  ♦  SUMCASE 

GO  TO  5560 

Z (K, 12 )  -  Z ( K, 1 2 )  +  SUMCASE 

CONTINUE 

GO  TO  6000 

CONTINUE 

IF  (  J  .GT.  1  )  GO  TO  5605 
IF  (  INTPL  .EQ.  1  )  GO  TO  5605 
IF  (  NDTRM3  .EQ.  0  )  GO  TO  5605 
N  -  3 

COORD ( 1 )  -  XI  -  XRE  -  XC 
C00RD12)  *  Y1  -  YRE  -  YC 
DO  5610  I  -  1,2 

IF  (  COORD(I)  .GE.  0.0  )  GO  TO  5615 
NCOORD  -  -COORD(I) /20.0  +  1.0 
IF  (  NCOORD  .GT.  50  )  NCOORD  -  50 
NMINUS (NCOORD,I,N)  -  NMINUS (NCOORD,I,N)  +  1 
GO  TO  5610 

NCOORD  -  COORD( IJ/20.0  +  1.0 

IF  (  NCOORD  .GT.  50  )  NCOORD  •  50 

NPLUS (NCOORD, I,N)  -  NPLUS (NCOORD, I, N )  +  1 

CONTINUE 

CONTINUE 

IF  (  NADJST  .NE.  2  )  GO  TO  5625 
IF  (  Y1  .LT.  YBASE  )  GO  TO  5635 
YC1  -  YC1  -  DROP 
GO  TO  5800 
CONTINUE 

YC1  *  YC1  ♦  DROP 
GO  TO  5800 
CONTINUE 
CALL  NRAN (R1,R2 ) 
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EXPLANATIONS  (CONTINUED) 


830-  B84 
B30-  B32 

833-  B45 

B47-  B53 

855 

856 

857-  860 

B61—  881 

B61-  B64 

B65-  B72 

873-  B81 

BB2-  BB4 


CYCLE  FOP.  EACH  ROUND  THAT  IS  FIRED#  FLIES  RELIABLY#  AND  MISSES  TARGET. 

PROGRAM  LINES  B33  THROUGH  B45  INVOLVE  SPECIAL  CALCULATIONS  ASSOCIATED  WITH 
NDTRM3  EQUALLING  1.  THESE  CALCULATIONS  CANNOT  BE  ATTEMPTED  IF  ROUND  INDEX  J 
INDICATES  THAT  SQNE  SUBSEQUENT  ROUND  RATHER  THAN  FIRST  ROUND  IS  BEING 
PROCESSED  OR  IF  RANGE  INTERPOLATION  OF  TARGET  VULNERABILITY  DATA  IS  NEEDED. 

M  EQUALLING  3  CORRESPONDS  TO  ALL  FIRST  ROUNDS  THAT  MISS  TARGET.  MISS  DISTANCES 
FOR  SUCH  ROUNDS  ARE  USED  FOR  SETTING  COUNTERS  IN  NMINUS  AND  NPLUS  ARRAYS  IN 
SAME  MANNER  THAT  MISS  DISTANCES  FOP.  FIRST  ROUNDS  HITTING  TARGET  ARE 
PROCESSED  BY  PROGRAM  LINES  714  THROUGH  731. 

IGNORE  STATEMENTS  APPLICABLE  ONLY  FOR  NAD J  ST  EQUALLING  2. 

RANDOM  NORMAL  DEVIATES  ARE  OBTAINED  FOR  POSSIBLE  LATER  USE  IN  PROGRAM  LINES 
B5B  AMD  859  OR  LINES  B75  AND  B76. 

PROGRAM  LIMES  B57  THROUGH  860  APPLY  FOR  BASIC  ADJUSTMENT  PROCEDURE  AND  ARE 
SKIPPED  IF  NADJUST  IS  1. 

IF  NAD JST  EQUALS  0#  RANDOM  NUMBER  IS  COMPARED  TO  PROBABILITY  PROBS  OF  ROUND 
BEING  SENSED.  IF  R  1UND  IS  NOT  SENSED#  BURST-ON-TARGET  ADJUSTMENT  IS  NOT  MADE 
AND  PROCESSING  CONTINUES  AT  PROGRAM  STATEMENT  5700.  OTHERWISE  DETERMINE  ADJX 
AND  ADJY.  THESE  T1J0  QUANTITIES#  RELATED  TO  BUP.ST-ON-TARGET  ADJUSTMENT#  ARE 
EXPLAINED  IN  CONNECTION  WITH  PROGRAM  LINES  882  THROUGH  884  WHERE  PROCESSING 
CONTINUES. 

ADJX  AND  ADJY  ARE  DETERMINED  FOR  NADJST  EQUALLING  1.  TWO  QUANTITIES#  RELATED 
TO  BURST-ON-TARGET  ADJUSTMENT#  ARE  EXPLAINED  IN  CONNECTION  WITH  PROGRAM 
LINES  8B2  THROUGH  884. 

MHIGH  EQUALS  0  IF  HISSING  ROUND  IS  SHORT  OF  TARGET#  1  OTHERWISE.  HIGH  AND 
SHORT  ARE  REAL  FORMS.  PRSCG  IS  SET  TO  APPLICABLE  PROBABILITY  OF  SENSING  FOR 
GUNNER  AND/OR  COMMANDER. 

RANDOM  NUMBER  RANF(DUM)  IS  STORED  IN  RANDOM  FOR  REUSE  IN  PROGRAM  LINE  87B. 
COMPARE  RANDOM  TO  PROBABILITY  OF  ROUND  BEING  SENSED  BY  GUNNER  AND/OR 
COMMANDER.  PROCESSING  CONTINUES  AT  PROGRAM  STATEMENT  5655  IF  ROUND  IS 
SENSED.  OTHERWISE  NUMBER  NFAIL  OF  UNSENSED  ROUNDS  IS  UPDATED.  IF  NFAIL  THEN 
EQUALS  1#  DROP  AIMPOINT  BY  HALF  TARGET  HEIGHT  BEFORE  CONSIDERATION  OF  NEXT 
ROUND.  TWO  FAILURES  TO  SENSE  CAUSE  THIS  OUTCOME  TO  BE  RECORDED  IN  AF AIL 
ARRAY  AND  ENGAGEMENT  TO  BE  ENDED. 

SET  SGX  AND  SGY  TO  APPLICABLE  HORIZONTAL  AND  VERTICAL  SENSING  ERROR  STANDARD 
DEVIATIONS  FOR  GUNNER#  DETERMINE  CORRESPONDING  PARTICULAR  ERRORS  ADJX  AND 
ADJY#  AND  SET  PRSG  TO  APPLICABLE  PROBABILITY  OF  ROUND  BEING  SENSED  BY 
GUNNER.  COMPARE  RANDOM  TO  GUNNER  SENSING  PROBABILITY  TO  ESTABLISH  WHETHER 
GUNNER  OR  COMMANDER  SHOULD  BE  CREDITED  WITH  SENSING.  IF  GUNNER  SENSED  ROUND# 
ADJX  AND  ADJY  ARE  ALREADY  CORRECT.  OTHERWISE  ADJUSTMENTS  ARE  MADE  TO  ACCOUNT 
FOR  ADDITIONAL  ERRORS  ASSOCIATED  WITH  SENSING  OF  ROUND  BY  COMMANDER  WITH 
INFORMATION  TRANSMITTAL  TO  GUNNER. 

THESE  STATEMENTS  COMPLETE  ACCOUNTING  FOR  BURST-ON-TARGET  ADJUSTMENT  BASED  ON 
SENSING.  UPDATED  VALUFS  OF  XC1  AMD  YC1  ARE  COORDINATES  OF  POINT  IN  TARGET 
COORDINATE  SYSTEM  WHERE  NEXT  ROUND  WOULD  IMPACT  IF  THERE  WERE  NO  RANDOM 
ERRORS.  XC1‘  IS  OBTAINED  BY  COMBINING  1)  HORIZONTAL  COORDINATE  XC  OF  INTENDED 
AIMPOINT#  2)  DIFFERENCE  XC1-X1  WHICH  INVOLVES  XC1  AND  XI  VALUES  FROM  EARLIER 
CALCULATIONS  AND  WHICH#  AS  CAN  BE  SEEN  FROM  PROGRAM  LINE  668#  CAUSES 
SPECIFIC  HORIZONTAL  RANDOM  ERRORS  THAT  AFFECTED  IMPACT  POINT  OF  SENSED 
MISSING  ROUND  TO  BE  REMOVED#  AND  3)  ADJX  VALUE  RELATED  TO  BURST-ON-TARGET 
ADJUSTMENT.  WHEN  NADJST  EQUALS  0#  ADJX  HAS  BEEN  OBTAINED  BY  APPLYING  R1  TO 


73 


PROGRAM  LISTING  (CONTINUED) 


IF  (  NADJST  .EQ.  1  )  GO  TO  5645 

IF  (  RANF(DUM)  .GT.  PROBS  )  GO  TO  5700 

ADJX  -  R1  *  SIGXS 


ADJY  -  R2  *  SIGYS 

860 _  GO  TO  5650 

5533  NHIGH  »  1  -  NSHORT 

HIGH  -  NHIGH 
SHORT  -  NSHORT 

PRSCG  -  HIGH+PGCH  +  SHORT+PGCS 
T55  RANDOM  -  RANFIDUM) 

IF  (  RANDOM  .LE.  PRSCG  )  GO  TO  5655 
NFAIL  *  NFAIL  +  1 
IF  (  NFAIL  .EQ.  2  )  GO  TO  5665 
YC1  -  YC1  -  DROP 
870  GO  TO  5800 

5665  AFAIL(J)  -  AFAIL(J)  +  1.0 

_  GO  TO  6000 

5635  SGX  *  HIGH*SGHX  +  SHORT*SGSX 

SGY  *  HIGH*SGHY  +  SHORT*SGSY 
875  ADJX  -  R1  *  SGX 

ADJY  -  R2  *  SGY 

PRSG  -  HIGH*PGH  +  SHORT*PGS 

IF  (  RANDOM  .LE.  PRSG  )  GO  TO  5650 

CALL  NRAN (R3*R4 ) 

880  ADJX  -  ADJX  +  R3*CDRX 

_ ADJY  -  ADJY  +  R4*CDRY 

5653  XC1  *  XC1  -  XI  +  XC  -  ADJX 

YC1  -  YC1  -  Y1  +  YC  -  ADJY 

_  FAC  *  0.0 

885  5700  IF  (  PASSN  .LE.  0.0  )  GO  TO  5800 

IF  (  NADJST  .NE.  0  )  GO  TO  5800 
IF  (  NBETA  .EQ.  MANGLE  )  GO  TO  5735 
Z(J*11)  -  Z(J*11)  +  SUMCASE 
GO  TO  5800 

890 _ 5735  Z  ( J>12 )  *=  Z(J*12>  +  SUMCASE 

5333  IF  (  J  .LT.  NRDS  )  GO  TO  5200 

_ 6000  CONTINUE 

IF  (  ISAMP  .LT.  NEND1  )  GO  TO  5100 

_ _ IF  {  INTPL  .EQ.  0  )  GO  TO  6005 

333  IF  (  ISAMP  .GT.  NEND1  )  GO  TO  6015 

ISPLIT  -  ISPLIT  +  1 

_  GO  TO  9010 

6015  IF  (  ISAMP  .LT.  NS  AMP  )  GO  TO  5100 
ISPLIT  =  ISPLIT  +  1 

900 _  IF  {  ISPLIT  .EQ.  3  )  GO  TO  9010 

_ 6005  CONTINUE 

_  IF  (NADJST  .EQ.  0  )  GO  TO  7000 

IFAIL  =  AFAIL(l) 

WRITE  (  6*1020  ) 

905  WRITE  (  6*6020  )  1* IFAIL 

6020  FORMAT  (  10X*  21 10  ) 

SUM  *  0.0 

DO  6030  12  *  2*  NRDS 
SUM  -  SUM  +  AF AIK  12-1) 

910  IFAIL  *  AFAIKI2)  +  SUM 

WRITE  (  6*6020  )  I2*IFAIL 
6030  CONTINUE 
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EXPLANATIONS  (CONTINUED) 


885-  890 

891 

893-  901 

893-  8<H 

895-  897 

898-  900 

903-  912 

914 

915-113Z 


134-1135 


HORIZONTAL  SENSING  ERROR  STANDARD  DEVIATION  SIGXS.  WHEN  NADJST  EQUALS  1* 

ADJ X  HAS  BEEN  DETERMINED  IN  SIMILAR  FASHION  WITH  CONSIDERATION  OF  POSSIBLY 
DIFFERENT  ERRORS  FOR  GUNNER  AND  COMMANDER.  MINUS  SIGN  IS  USED  WITH  ADJX 
BECAUSE*  IF  ROUND  IS  SENSED  AS  BEING  OFF  IN  ONE  DIRECTION  (FOR  EXAMPLE* 
RIGHT  AND  DOWN)*  FIRING  WEAPON  IS  MOVED  IN  OPPOSITE  DIRECTION  (IN  THIS 
EXAMPLE*  LEFT  AND  UP).  TERMS  IN  FORMULA  FOR  YC1  PARALLEL  THOSE  IN  XC1 
FORMULA.  FAC  NEEDS  TO  BE  RESET  SO  THAT  LAY  ERRORS  WILL  NO  LONGER  BE  TREATED 
EXPLICITLY  SINCE  THEY  ARE  INCLUDED  IN  SENSING  ERRORS. 

WHEN  PASSN  HAS  NONZERO  VALUE  AND  ENGAGEMENT  INVOLVES  BASIC  ADJUSTMENT  METHOD* 
UPDATE  Z  ARRAY  TO  ACCOUNT  FOR  CUMULATIVE  PASSENGER  PERSONNEL  CASUALTIES 
CAUSED  BY  ROUNDS  FIRED  BEFORE  MISSING  ROUND  OF  CURRENT  CONCERN. 

PROCEED  WITH  NEXT  ROUND  UNLESS  MAXIMUM  NUMBER  OF  ROUNDS  ALLOWED  HAVE  ALREADY 
BEEN  ACCOUNTED  FOR. 


PROCEED  WITH  NEXT  SAMPLE  ENGAGEMENT  UNLESS  ALL  ENGAGEMENTS  NEEDED  AT  THIS 
STAGE  OF  CALCULATIONS  HAVE  BEEN  COMPLETED. 

SKIP  PROGRAM  LINES  8°5  THROUGH  900  IF  ISAMP  EQUALS  NEND1  AND  RUN  DOES  NOT 
INVOLVE  RANGE  INTERPOLATION  OF  VULNERABILITY  DATA. 

WHEN  ISAMP  EQUALS  NEND1*  RESET  ISPLIT  TO  INDICATE  FIRST  PART  OF  ENGAGEMENT 
SIMULATIONS  FOR  IC  VALUE  OF  INTEREST  AND  ASSOCIATED  FORWARD  OR  REVERSE  ANGLE 
HAS  BEEN  COMPLETED.  AFTER  BEING  RESET*  ISPLIT  EQUALS  1  FOR  FORWARD  ANGLE  OR 
2  FOR  REVERSE  ANGLE.  CONTINUE  PROCESSING  AT  PROGRAM  STATEMENT  9010.  PROGRAM 
LINE  895  CAUSES  NEXT  TWO  LINES  TO  BE  SKIPPED  WHEN  ISAMP  EXCEEDS  NENDl. 

PROCEED  WITH  NEXT  ENGAGEMENT  UNLESS  ALL  MS AMP  ENGAGEMENTS  HAVE  BEEN  COMPLETED. 
RESET  ISPLIT  TO  3  FOR  FORWARD  ANGLE  OR  4  FOR  REVERSE  ANGLE  TO  INDICATE 
SECOND  PART  OF  ENGAGEMENT  SIMULATIONS  FOR  IC  AND  RELATED  FORWARD  OR  REVERSE 
ANGLE  HAS  BEEN  DONE.  CONTINUE  AT  PROGRAM  STATEMENT  9010  IF  CALCULATIONS  FOR 
REVERSE  ANGLE  ARE  STILL  INCOMPLETE. 

IFAIL  ARRAY  CONTAINS  CUMULATED  VALUES  CORRESPONDING  TO  AFAIL  ARRAY.  IFAIL(I2) 
IS  NUMBER  OF  ENGAGEMENTS  IN  WHICH  SECOND  FAILURE  TO  SENSE  THAT  CAUSED  END  OF 
ENGAGEMENT  OCCURRED  ON  ROUND  LESS  THAN  OR  EQUAL  TO  12.  WRITE  IFAIL  ARRAY 
QUANTITIES. 

PROGRAM  LINES  Q15  THROUGH  1132  INVOLVE  SPECIAL  CALCULATIONS  ASSOCIATED  WITH  AT 
LEAST  ONE  OF  NDTRM1*  NDTRM2*  MDTRM3*  NDTRM4*  AND  NDTRM5  EQUALLING  1.  THESE 
CALCULATIONS  CANNOT  BE  ATTEMPTED  IF  PANGE  INTERPOLATION  OF  TARGET 
VULNERABILITY  DATA  HAS  BEEN  INVOLVED. 

THESE  STATEMENTS  ESSENTIALLY  CONSTITUTE  COMPUTER  ROUTINE  (ALTHOUGH  STATEMENTS 
HAVE  NOT  BEEN  STRUCTURED  AS  SEPARATE  ROUTINE  IN  THIS  PROGRAM)  FOP  PROCESSING 
HORIZONTAL  OR  VERTICAL  OFFSET  DISTANCES  BASED  ON  IMPACT  POINTS  OF  HITTING  OR 
MISSING  ROUNDS  IN  TARGET  PLANE.  OBJECTIVE  OF  PROCESSING  IS  TO  ESTIMATE  MEAN 
AND  OBTAIN  THREE  ESTIMATES  OF  STANDARD  DEVIATION  FOR  NORMAL  DISTRIBUTIONS 
TENTATIVELY  ASSUMED  TO  FIT  SUCH  OFFSET  DATA.  AFTER  COMPLETION  OF  COMPUTER 
RUN*  ANALYST  CAM  STUDY  THREE  SETS  OF  PARAMETERS  TO  JUDGE  HOW  CLOSELY  NORMAL 
DISTRIBUTION  APPLIES  AND  TO  MAKE  BEST  ESTIMATE  FOR  PARAMETERS  OF  INTEREST. 
PAPER  EXPLAINING  PROCESSING  OF  OFFSET  DISTANCES  WAS  PRESENTED  AT  1979  ARMY 
NUMERICAL  ANALYSIS  AMD  COMPUTERS  CONFERENCE.  THIS  PAPER  (WITH  MINOR  CHANGES) 
IS  INCLUDED  AS  APPENDIX  B  OF  THIS  REPORT.  N  EQUALLING  1,  2*  OR  3  IS  RELATED 
TO  HITTING  OR  MISSING  FIRST  ROUNDS  AS  ALREADY  EXPLAINED.  PROVISION  HAS  ALSO 
BEEN  MADE  FOR  N  VALUES  OF  4  AND  5  TO  BE  AVAILABLE  FOR  POSSIBLE  FUTURE  USE. 
CALCULATIONS  INVOLVING  ANY  PARTICULAR  N  VALUE  ARE  DONE  ONLY  IF  CONTROL 
QUANTITY  KDMSTC(N)  IS  NOT  0. 

PROGRAM  LINES  1136  THROUGH  1295  INVOLVE  DETERMINISTIC  CALCULATIONS  OF  HIT 
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TIT 

920 

925 

930 

935 

940 

945 

950 

955 

960 

965 


PROGRAM  LISTING  (CONTINUED) 


7000  CONTINUE 

_ __  IF  (  ISPLIT  .EQ.  4  )  GO  TO  7400 

N  *  0 

7005  IF  (  N  .EQ.  5  )  GO  TO  7400 
N  •  N  +  1 

IF  (  KDMSTC(N)  .EQ.  0  )  GO  TO  7005 

DO  7010  1-1,2 

AVRG  -  0.0 

TMINUS  -  0.0 

TPLUS  -  0.0 

KTIMES  -  0 

DO  7020  J  -  1,50 

TERM1  -  NMINUS! J,I,N) 

TERM2  =■  NPLUS ( J» I, N) 

TMINUS  -  TMINUS  ♦  TERM1 
7020  TPLUS  -  TPLUS  +  TERM2 

TOTAL! I)  -  TMINUS  ♦  TPLUS 
N TOTAL  -  TOTAL (I) 

T5000  ■  TOTAL! I)  *0.5 
N5000  -  T5000 

IF  (  T5000  .NE.  TMINUS  )  GO  TO  7024 
DO  7022  K  -  1,50 
NPOSIK)  •  NPLUS (K, I,N) 

7022  NNEGIK)  *  NMINUStK, I, N) 

GO  TO  7200 

7024  IF  (  T5000  .LT.  TMINUS  )  GO  TO  7025 
DIFF  »  T5000  -  TMINUS 
DO  7030  J  *  1,50 
PLUS  -  NPLUS ( J, I,N) 

IF  (  PLUS  .LT.  DIFF  )  GO  TO  7035 
IFRCTN  -  1000.0  *  DIFF  /  PLUS 
FRCTN  -  FLOAT(IFRCTN)  /  1000.0 
AVRG  -  AVRG  ♦  FRCTN+20.0 
JEND  =  100  -  KTIMES 
LTIMES  -  50  -  KTIMES 
LEND  -  1  +  KTIMES 
DO  7040  K  -  1,50 
NNEG(K)  -  0 
7040  NPOSIK)  ■  0 
NCMLTN  -  0 
NCMLTP  -  0 
L  -  0 

7050  IF  !  L  .EQ.  100  )  GO  TO  7200 
L  -  L  +  1 
LL  -  L 

IF  (  L  .LE.  50  )  LL  -  51  -  L 
M  -  L 

IF  l  L  .GT.  50  )  M  -  L  -  50 
IF  l  LL  .LE.  JEND  )  GOTO  7055 
NPOSIM)  -  0 
GO  TO  7060 

7055  IF  1  LL  .LT.  JEND  )  GOTO  7065 

NTERM  -  ( IFRCTN*NPLUS ( 50, I, NJ+500)  /  1000 
NPOS (M )  -  NPLUS  !50,I,N)  -  NTERM 
NCMLTP  -  NCMLTP  ♦  NPOS!M) 

GO  TO  7060 

7065  IF  (  LL  .LE.  LTIMES  )  GOTO  7075 
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970 

975 

980 

985 

990 

995 

1000 

1005 

1010 

1015 

1020 

1025 


PROGRAM  LISTING  (CONTINUED) 


IF  {  LL  ,LE.  50  )  GOTO  7070 

NTERM1  -  CIFRCTN*NPLUS(L-LTIMES,I,N)*500)  /  1000 
NTERN2  *  { IFRCTN*NPLUS (L-LTIMES+1* I# N ) +500)  /  1000 
NPOS(M)  *  NPLUS(L-LTIMES,I,N)  -  NTERM1  ♦  NTERM2 
NCMLTP  -  NCMLTP  +  NPOS(M) 

GO  TO  7060 

7070  NTERM1  -  (I FRCTN*N PLUS (LEND-M,I,N)+500)  /  1000 

NTERM2  -  { IFRCTN*NPLUS ( LEND-M+1, I,N)+500)  /  1000 
NNEG(M)  *  NPLUS(LEND-M»I,N1  -  NTERM1  +  NTERM2 
NCMLTN  -  NCMLTN  ♦  NNEG(M) 

GO  TO  7060 

7075  IF  {  LL  .LT.  LTIMES  )  GOTO  7085 

NTERM1  -  ( IFRCTN*NMINUS (l,I,N)+500)  /  1000 
NTERM2  -  (IFRCTN*NPLUS(1, I,N)+500)  /  1000 
NNEG(M)  -  NMINUS(1,I,N>  -  NTERM1  ♦  NTERM2 
NCMLTN  *  NCMLTN  +  NNEG(M) 

GO  TO  7060 

7085  IF  (  LL  .EQ.  1  )  GO  TO  7095 

NTERM1  -  ( IFRCTN+NMINUS (LTIMES-LL+l,I,N)+500)  /  1000 
NTERM2  -  (IFRCTN*NMINUS(LTIMES-LL,I,N)+500)  /  1000 
NNEG(L-KTIME S+l )  -  NMINUS (LTIMES-LL+1, I* N )  -  NTERM1  +  NTERM2 
NCMLTN  *  NCMLTN  +  NNEGIL-KTIMES+l ) 

GO  TO  7060 

7095  NNEG150)  -  ( IFRCTN*NMINUS (49-K TIMES, I,N )+500)  /  1000 
NXTRA  *  KTIMES  +  1 
DO  7100  JXTRA  -  1, NXTRA 

7100  NNEG(50)  >  NNEG ( 50 )  ♦  NMINUS149-KTIMES+JXTRA, I, N) 

NTERM  -  (IFRCTN*NMINUS(50,I,N)+500)  /  1000 
NNEG(50)  *  NNEG ( 50 )  «■  NTERM 
NCMLTN  -  NCMLTN  +  NNEG ( 50 ) 

7060  IF  (  NCMLTP  .EQ.  N TOTAL  )  L  *  100 
IF  (  L  .GT.  50  )  GO  TO  7050 
IF  (  NCMLTN  .LT.  N5000  )  GO  TO  7050 
NCMLTP  «=  NCMLTN 
L  -  50 
GO  TO  7050 

7035  DIFF  «  DIFF  -  PLUS 
KTIMES  *  KTIMES  +  1 
AVRG  -  AVRG  ♦  20.0 

7030  CONTINUE 

7025  DIFF  «  TMINUS  -  T5000 
DO  7110  J  -  1,50 
AMINUS  -  NMINUS (J, I, N) 

IF  (  AMINUS  .LT.  DIFF  )  GO  TO  7115 

IFRCTN  »  1000.0  *  DIFF  /  AMINUS 

FRCTN  «  FLOATII FRCTN)  /  1000.0 

AVRG  «.  AVRG  -  FRCTN* 20.0 

JEND  -  100  -  KTIMES 

LTIMES  -  50  -  KTIMES 

LEND  «  1  +  KTIMES 

DO  7120  K  *  1,50 

NPOSIKJ  «  0 

7120  NNEG(K)  -  0 
NCMLTP  -  0 
NCMLTN  -  0 
L  *  0 

7130  IF  {  L  .EQ.  100  )  GO  TO  7200 
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1030 

1035 

1040 

1045 

1050 

1055 

1060 

1065 

1070 

1075 

1000 


PROGRAM  LISTING  (CONTINUED) 


L  -  L  +  1 
LL  -  L 


IF 

(  L  .LE. 

50  )  LL 

•  51 

-  L 

M  - 

L 

IF 

(  L  .GT. 

50  )  M 

-  L  - 

50 

IF 

(  LL  .LE. 

JEND  ) 

GOTO 

7135 

NNEG(M)  -  0 

GO 

TO  7140 

IF 

(  LL  .LT. 

JEND  ) 

GOTO 

7145 

NTERM  -  ( IFRCTN*NM INUS ( 50* I jN) +500)  /  1000 
NNEG(M)  -  NHINUS(50*I*N)  -  NTERM 
NCMLTN  -  NCMLTN  +  NNEG(M) 

GO  TO  7140 

7145  IF  (  LL  .LE.  LTIMES  )  GOTO  7155 
IF  (  LL  .LE.  50  )  GOTO  7150 

NTERM1  *  (IFRCTN+NMINUS (L-LTIMES* I*N)+500)  /  1000 
NTERM2  -  ( IFRCTN+NM INUS {L-LTIMES+1*I»N)+500)  /  1000 
NNEG(M)  -  NMINUS(L-LTIMES*I*N)  -  NTERM1  +  NTERM 2 
NCMLTN  -  NCMLTN  +  NNEG(M) 

GO  TO  7140 

7150  NTERM1  -  (IFRCTN*NMINUS ( LEND-M* I*  N)+500 )  /  1000 

NTERM2  -  (IFRCTN+NMINUS (LEND-M+1# I, N ) +500 )  /  1000 
NPOS(M)  -  NMINUS (LEND-M* I# N )  -  NTERM1  +  NTERM2 
NCMLTP  -  NCMLTP  +  NPOS(M) 

GO  TO  7140 

7155  IF  (  LL  .LT.  LTIMES  )  GOTO  7165 

NTERM1  -  ( IFRCTN*NPLUS ( 1# I*N) +500)  /  1000 
NTERM2  -  (IFRCTN+NMINUS (l*I*N)+500)  /  1000 
NPOS(M)  -  NPLUS (1* I*N)  -  NTERM1  +  NTERM2 
NCMLTP  -  NCMLTP  ♦  NPOS(M) 

GO  TO  7140 

7165  IF  (  LL  .EQ.  1  )  GO  TO  7175 

NTERM1  *  ( IFRCTN+NPLUS (LTIMES-LL+l*I*N)+500)  /  1000 
NTERM2  -  (IFRCTN*NPLUS (LTIM  ES-LL* I*  N) +500 )  /  1000 
NPOS(L— KTIMES+1)  -  NPLUS (LTIMES-LL+1#I*N)  -  NTERM1  +  NTERM2 
NCMLTP  «  NCMLTP  +  NPOS (L-KTIMES+1 ) 

GO  TO  7140 

7175  NPOS ( 50 )  -  ( IFRCTN+NPLUS ( 49— KTI MES*I#N)  +  500)  /  1000 
NXTRA  -  KTIMES  +  1 
DO  7180  JXTRA  -  1* NXTRA 

7180  NPOS( 50 )  -  NPOS  (50  )  +  NPLUS ( 49-KTIMES +JXTRA* I*N) 

NTERM  -  (IFRCTN*NPLUS(50,I*N)+500)  /  1000 
NPOS (50  )  -  NPOS ( 50 )  +  NTERM 
NCMLTP  -  NCMLTP  +  NPOS (50) 

7140  IF  (  NCMLTN  .EQ.  N TOTAL  )  L  -  100 
IF  (  L  .GT.  50  )  GO  TO  7130 
IF  (  NCMLTP  .LT.  N5000  )  GO  TO  7130 
NCMLTN  -  NCMLTP 
L  -  50 
GO  TO  7130 

7115  DIFF  -  DIFF  -  AMINUS 
KTIMES  -  KTIMES  +  1 
AVRG  -  AVRG  -  20.0 

7110  CONTINUE 

7200  CONTINUE 

DO  7310  K  -  1#  50 

IF  (  K  .GT.  1  )  GO  TO  7315 
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PROGRAM  LISTING  (CONTINUED) 


NSUM1  -  NPOS(K) 

1085  NSUM2  -  NNEG(K) 

GO  TO  7320 

7315  NSUM1  *  NSUM1  +  NPOS(K) 

NSUM2  -  NSUM2  +  NNEG(K) 

7320  NCMPOS(K)  -  NSUM1 

1090  NCMNEG(K)  -  NSUM2 

NRFLCT(K)  -  NSUM1  +  NSUM2 
7310  CONTINUE 

DO  7330  J  -  1,50 
SIGMAX  -  J  *  20 

1095  IF  (  J  .GT.  1  )  GO  TO  7335 

SMFRQl  -  0.0 

SMFRQ2  -  NRFLCT(l)  *  10000  /  NTOTAL 
K5  -  5 
K10  -  10 

1100  GO  TO  7350 

7335  SMFRQl  -  SMFRQ2 

IF  (  J  .EQ.  50  )  GO  TO  7345 
SMFRQ2  -  NRFICT(J)  *  10000  /  NTOTAL 
GO  TO  7350 

1105  7345  SMFRQ2  -  10000.0 

7350  FRQNCY  «  SHFRQ2  -  SMFRQl 

IF  (  K5  .NE.  5  )  GO  TO  7355 
IF  {  SMFRQ2  .LT.  3829.2  )  GOTO  7330 
D05  -  (SMFRQ2-3829.2 )  /  FRQNCY 

1110  SIG05  -  (SIGMAX— 005*20.0)  /  0.5 

K5  -  0 

7355  IF  {  K10  .NE.  10  )  GO  TO  7365 

IF  {  SMFRQ2  .LT.  6826.8  )  GOTO  7330 
DIO  -  (SMFRQ2— 6 826.8)  /  FRQNCY 
1115  SIG10  -  SIGMAX  -  D10*20.0 

K10  *=  0 

7365  IF  (  SMFRQ2  .LT.  8663.8  )  GOTO  7330 
D15  -  ( SMFRQ 2-8663. 6)  /  FRQNCY 
SIG15  -  (SIGMAX— D1 5*20. 0)  /  1.5 
1120  GO  TO  7370 

7330  CONTINUE 
7370  WRITE  (  6,1020  ) 

WRITE  (  6,7380)  N, I, SIG05, S IG10, S IG15 
7380  FORMAT  (  10X, 26H  N, I, SIG05, SIG10, SIG15  ,2I10,3F10.4  ) 
1125  BIASXY ( I,N)  -  AVRG 

IF  (  I  .EQ.  1  )  SSQRRE  -  SIGXL**2  ♦  SIGXR**2 
IF  (  I  .EQ.  2  )  SSQRRE  -  SIGYL**2  +  SIGYR**2 
SIGMXY ( 1, I, N )  -  SORT ( SI G05**2+SSQRRE ) 

SIGMXY(2, I,N )  »  SQRT( SI G10**2+SSQRRE ) 

1130  SIGMXY(3,I,N)  -  SORT (SIG15**2+SSQRRE ) 

7010  CONTINUE 

_  GO  TO  7005 

7406'  CONTINUE 

IF  (  NSMDTR  .EQ. 

1135 _  IF  (  ISPLIT  .EQ. 

NCLLS1  -  99999 

_  NCLLS2  -  99999 

IF  (  NSMDTR  .EQ. 

_  IF  (  ISPLIT  .EQ. 

1140  NCLLS1  -  0 


0  .AND.  NHTKLL  .EQ.  0  )  GOTO  8000 
4  .AND.  NHTKLL  .EQ.  0  )  GO  TO  8000 


0  )  GO  TO  7405 
4  )  GO  TO  7405 
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PROGRAM  LISTING  (CONTINUED) 


- -  DO  7410  N  »  1/5 

BX (N )  *  BIAS XY ( 1/N ) 

BY ( N)  *  BIASXY(2/N) 

DO  7410  K  *  1*  3 

1145  PH(K/N )  -  0.0 

SIGX(K/N)  »  SIGMXYIK/ 1» H ) 
SIGY(K/N)  »  SIGMXY (K/2»  N) 

DO  7410  L  *  1*4 
PKILKL/K/N)  *  0.0 
1150  7410  CONTINUE 

- 75153  IF  <  NHTKLt  .EQ.  0  )  GO  TO  7415 

-  NCLLS2  «  0 

-  DO  7420  N  -  1/5 

PPH(N)  »  0.0 

1155  IF  (  N  .EQ.  5  )  PPH(N)  *  1.0 

DO  7420  L  -  1/4 
PKSHOT ( L/ N )  »  0.0 
PPKHIT ( L/N )  -  0.0 


1160 


7420  CONTINUE 

7415  DO  7430  J  *  1/NY 

-  IF  <  AY ( J )+H  .LE.  YBASE  )  GO  TO  7430 


Y2  -  AY ( J )  ♦  H  -  YC 
Y1  -  AY( J)  -  H  -  YC 

IF  (  Y1  .LT.  YBASE-YC  >  Y1  -  YBASE  -  YC 
TT65  DO  7440  I  -  1/NX 

- - K22  -  IK ( J/ 1 ) 

IF  (  K22  .LE.  0  )  GO  TO  7440 

“  X2  -  AX(I)  +  H  -  XC 

XI  -  AX (I)  -  H  -  XC 

TITO  IF  (  NPRHIT  .EQ.  1  )  GO  TO  7445 

-  N  -  -49 

DO  7450  L  -  1/4 

IF  (  NBETA  .EQ.  NANGLE  )  GO  TO  7452 
1ST  -  SHIFT(MPK(K22/IA)/N)  .AND.  MASK11 
1175  ST  -  1ST 

PK ( L)  -  0.001  *  ST 

GO  TO  7450 

7452  IST1  *  SHIFT(MPK(K22/IB)/N)  .AND.  MASK11 
ST1  -  I ST1 

1180  PK(L*5 )  -  0.001  *  ST1 

7450  N  «  N  ♦  11 

- 7551>  IF  (  NSMDTR  .EQ .  0  )  GO  TO  7475 

IF  (  ISPLIT  .EQ.  4  >  GO  TO  7475 
-  N  -  0 

1185  7455  IF  (  N  .EQ.  5  )  GO  TO  7475 

N  -  N  +  1 

IF  (  KDMSTC(N)  .EQ.  0  )  GO  TO  7455 
- - -  DO  7460  K  -  1/3 

“  Q0FZ1  *  CNORM( (X2-BX(N) )/$IGX(K/N) > 

1190  Q0FZ2  -  CNORM( ( XI— BX(N ) )/SIGX(K/N ) ) 

Q0FZ3  -  CNORM( ( Y2-BY(N) )/SIGY(K/N ) ) 

Q0FZ4  *  CNORMU  Yl-BY(N)  )/SIGY(K»N)  ) 

POFXY  *  (Q0FZ1— Q0FZ2)  *  ( Q0FZ3-Q0FZ4 ) 

PH( K/N)  -  PH(K/N>  +  POFXY 

TT35 -  IF  t  NPRHIT  .EQ.  1  )  GO  TO  7460 

- - -  DO  7470  L  »  1/4 

PKCELL  -  PK< L ) 
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EXPLANATIONS  (CONTINUED) 


1136-1137 

1138-1139 

1140 

1141-1150 

1151 

1152 

1153-1159 

1160- 1227 

1161- 1164 

1165- 1226 

1166- 1167 

1168-1169 

1170 

1171-1181 

1182-1183 

1184-1187 

1188-1201 

118*5-1194 


PROBABILITIES  AND  KILL  PROBABILITIES.  LINES  NEED  TO  BE  SKIPPED  WHEN  NSMDTR 
AND  MHTKLL  CONTROLS  ARE  BOTH  0.  CONDITION  INVOLVING  ISPLIT  EQUALLING  4 
CAUSES  LIMES  TO  BE  SKIPPED  ONLY  FOR  RANGES  REQUIRING  INTERPOLATION  OF 
VULNERABILITY  DATA  IF  NSMDTR  IS  NOT  0. 

INITIAL  SETTINGS. 

PROGRAM  LINES  1140  THROUGH  1150  CAM  BE  INVOLVED  ONLY  IF  NSMDTR  EXCEEDS  0  AND 
RANGE  OF  CONCERN  DOES  NOT  REQUIRE  INTERPOLATION  OF  VULNERABILITY  DATA. 

NEW  INITIAL  SETTING  OVERRIDING  99999. 

INITIAL  SETTINGS. 

PROGRAM  LINES  1152  THROUGH  1159  CAN  BE  INVOLVED  ONLY  IF  NHTKLL  EXCEEDS  0. 

NEW  INITIAL  SETTING  OVERRIDING  99999. 

INITIAL  SETTINGS. 

CONSIDER  IM  TURN  EACH  HORIZONTAL  STRIP  OF  CELLS  IN  SMALLEST  RECTANGLE 
ENCLOSING  TARGET. 

DETERMINE  Y  COORDINATES  Y1  AND  Y2  CORRESPONDING  TO  LOWER  AMD  UPPER  EDGES 
RESPECTIVELY  OF  PARTICULAR  STRIP.  AY(J)  REPRESENTS  Y  COORDINATE  OF 
HORIZONTAL  CENTERLINE  OF  STRIP.  IGNORE  STRIP  IF  ITS  CENTERLINE  IS  BELOW  BASE 
OF  TARGET.  IF  NECESSARY*  RESET  Y1  TO  MAKE  IT  CORRESPOND  TO  TARGET  BASE. 

CONSIDER  IN  TURN  EACH  CELL  IN  PARTICULAR  HORIZONTAL  STRIP. 

K22  IS  USED  FOR  TEMPORARY  STORAGE  OF  IK(J*I)  VALUE  ASSOCIATED  WITH  PARTICULAR 
CELL.  IGNORE  CELL  WHEN  K22  EQUALLING  0  INDICATES  CELL  IS  NOT  PART  OF  TARGET. 

DETERMINE  X  COORDINATES  XI  AND  X2  CORRESPONDING  RESPECTIVELY  TO  LEFT  AND  RIGHT 
EDGES  OF  TARGET  CELL. 

SKIP  PROGRAM  LINES  1171  THROUGH  1181  IF  ONLY  HIT  PROBABILITIES  APE  INVOLVED. 

UNPACK  MPK  ARRAY  QUANTITY  THAT  CORRESPONDS  TO  TARGET  CELL  HIT  AND  ASSOCIATED 
TARGET  ORIENTATION  TO  RETRIEVE  PROBABILITIES  OF  M  KILL*  F  KILL*  M  OR  F  KILL* 
AND  K  KILL  PER  HIT.  STORE  THESE  IM  PK  ARRAY.  THESE  PROGRAM  LINES  ARE 
ESSENTIALLY  LIKE  LINES  753  THROUGH  763.  EXCEPT  FOP  EXPECTED  CASUALTIES  NOT 
BEING  OF  CONCERN  HERE*  EXPLANATIONS  FOR  LINES  753  THROUGH  772  APPLY. 

SKIP  PROGRAM  LINES  1184  THROUGH  1203  IF  NSMDTR  IS  0  OR  IF  RANGE  OF  CONCERN 
INVOLVES  INTERPOLATION  OF  VULNERABILITY  DATA. 

PROCESSING  INDICATED  BY  PROGRAM  LIMES  1188  THROUGH  1203  IS  DONE  OR  SKIPPED  FOR 
N  EQUALLING  1*  2*  3*  4*  OR  5  ACCORDING  TO  WHETHER  KDMSTC(N)  EXCEEDS  OR 
EQUALS  0. 

K  CORRESPONDS  TO  THREE  ALTERNATIVE  SETS  OF  STANDARD  DEVIATIONS  ESTIMATED  FOR 
HORIZONTAL  AMD  VERTICAL  OFFSET  DISTANCE  DATA  FROM  SIMULATED  ENGAGEMENTS. 

Q0FZ1  IS  PROBABILITY  THAT  PARTICULAR  PAIR  OF  BIAS  AND  STANDARD  DEVIATION 
VALUES  RESULTS  IN  X  COORDINATE  OF  PROJECTILE  IMPACT  POINT  NOT  EXCEEDING  X 
COORDINATE  OF  RIGHT  EDGE  OF  CELL  BEING  CONSIDERED.  Q0FZ2  IS  PROBABILITY  X 
COORDINATE  OF  IMPACT  POINT  DOES  NOT  EXCEED  THAT  OF  LEFT  EDGE  OF  CELL. 
DIFFERENCE  Q0FZ1-Q0FZ2  EQUALS  PROBABILITY  THAT  IMPACT  POINT  LIES  ON  VERTICAL 
LINE  THROUGH  SOME  CELL  POINT.  SIMILARLY*  DIFFERENCE  Q0FZ3-Q0FZ4  EQUALS 
PROBABILITY  THAT  IMPACT  POINT  LIES  ON  HORIZONTAL  LINE  THROUGH  SOME  CELL 
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1200 


1205 


1210 


1215 


1220 


1225 


1230 


1235 


1240 


1245 


1250 


7470 

7460 


7475 


7485 


74^5 


7510 

7480 

T5313 


7540 

7530 

7600 


7610 

7620 

7640 

7650 

7660 

7630 


IF  {  NBETA  .EQ.  NANGLE  )  PKCELL  -  PKtL+5) 
PKILL(L,K,N)  -  PKILL(L,K,N)  ♦  POFXY*PKCELL 
CONTINUE 


CONTINUE 

NCLLS1  -  NCLLS1  ♦  1 
GO  TO  7455 

IF  (  NHTKLL  .EQ.  0  )  GO  TO  7440 

DO  7480  N  *  1#  5 

IF  (  N  .EQ.  5  )  GO  TO  7485 

IF  {  SSIGX(N)  .EQ.  0.0  .OR.  SSIGY(N)  .EQ.  0.0  )  GO  TO  7480 
QQ0FZ1  -  CNORMt ( X2-BBX( N) )/ SSIGX (N ) ) 

QQ0FZ2  -  CNORMt (Xl-BBX(N) )/SSIGX(N) ) 

QQ0FZ3  -  CNORM( (Y2-BBY(N))/SSIGY(N) ) 

QQ0FZ4  -  CNORM< (Yl-BBY(N) )/SSIGY(N) ) 

PPOFXY  -  (QQ0FZ1-QQ0FZ2 )  *  (QQ0FZ3-QQ0FZ4) 


PPH(N)  *  PPHtN)  ♦  PPOFXY 
GO  TO  7490 

IF  (  NRDTYP  .EQ.  0  .AND.  NRANGE  .GT.  1  )  GO  TO  7480 
PPOFXY  -  1.0 

PPH(5)  -  PPH(5)  ♦  PPOFXY 

NCLLS2  -  NCLLS2  +  1 

IF  (  NPRHIT  .EQ.  1  )  GO  TO  7480 

DO  7510  L  *  1,  4 

PPCELL  -  PK  t  L ) 

IF  (  NBETA  .EQ.  NANGLE  )  PPCELL  -  PK(L+5) 

PKSHOT ( L,  N)  -  PKSHOT(L»N)  +  PPOFXY+PPCELL 
CONTINUE 


CONTINUE 

CONTINUE 

CONTINUE 

IF  (  NHTKLL  .EQ.  0  )  GO  TO  7600 

IF  (  NPRHIT  .EQ.  1  )  GO  TO  7600 

DO  7530  L  -  1*4 

DO  7540  N  -  1,5 

IF  (  PPHtN)  .GT.  0.000001  )  PPKHIT(L,N)  »  PKSHOT(L,N)  /  PPHtNI 

CONTINUE 

CONTINUE 

CONTINUE 

IF  (  NSMDTR  .EQ.  0  )  GOTO  7700 
IF  (  ISPLIT  .EQ.  4  )  GO  TO  7700 
WRITE  (  6,1020  ) 

WRITE  (  6,7610  )  (  BX(I), 1-1,5  ) 

FORMAT  (  10X, 11H  BX(5)  ,10X,5F10.4  ) 

WRITE  (  6,7620  )  (  BY  (I),  I- 1,5  ) 

FORMAT  (  10X, 11H  BYt 5)  ,10X,5F10.4  ) 

DO  7630  1-1,3 
WRITE  (  6,1020  ) 

WRITE  (  6,7640  )  (  SIGXt I, J ), J-l, 5  ) 

FORMAT  (  10X, 11H  SIGX<3,5)  ,10X,5F10.4  ) 

WRITE  (  6,7650  )  (  SIGY (I,J ),J-1, 5  ) 

FORMAT  (  10X, 11 H  SIGY(3,5)  ,10X,5F10.4  ) 

WRITE  (  6,7660  )  (  PHt I, J ) , J-l, 5  ) 

FORMAT  (  10X, 21H  PH(3,5)  ,5F10.4  ) 

CONTINUE 


IF  (  NPRHIT  .EQ.  1  )  GOTO  7700 
DO  7670  J  -  1,3 
WRITE  (  6,1020  ) 
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POINT.  PRODUCT  POFXY  OF  TWO  DIFFERENCES  IS  PROBABILITY  OF  HITTING  CELL.  EACH 
PH(  K,  N)  VALUE  IS  CUMULATIVE  HIT  PROBABILITY  FOR  ALL  CELLS  ALREADY  PROCESSED 
AND  EVENTUALLY  APPLIES  TO  ENTIRE  TARGET. 

1195 

SKIP  PROGRAM  LINES  1196  THROUGH  1200  IF  ONLY  HIT  PROBABILITIES  ARE  INVOLVED. 

1196-1200 

L  CORRESPONDS  TO  FOUR  KILL  CRITERIA  CONSIDERED  IN  TURN.  PKCELL  IS  USED  FOR 
TEMPORARY  STORAGE  OF  PK(L)  OR  PK(L+5).  THESE  TWO  VALUES  APPLY  TO  FORWARD  OR 
REVERSE  ANGLE  RESPECTIVELY.  PRODUCT  OF  POFXY  AND  PKCELL  IS  PROBABILITY  OF 
HITTING  PARTICULAR  CELL  AND  THEREBY  KILLING  TARGET.  EACH  PKILL ( L, K»  N)  VALUE 

IS  CUMULATIVE  PROBALITY  ACCOUNTING  FOR  ALL  CELLS  ALREADY  PROCESSED  AND 
EVENTUALLY  EQUALS  PROBABILITY  OF  HIT  AND  KILL  FOP.  ENTIRE  TARGET. 

v  1202 

NCLLS1  IS  UPDATED  AND  EVENTUALLY  EQUALS  PRODUCT  OF  1)  TOTAL  NUMBER  OF  TARGET 
CELLS  AND  2)  NUMBER  OF  NONZERO  QUANTITIES  IN  SET  KDMSTC(l),  KDMSTC(2), 

KDMSTC ( 3) ,  AMD  KDNSTC(4). 

1203 

PROCESSING  HAS  BEEN  COMPLETED  FOR  PARTICULAR  VALUE  OF  N  AMD  MAY  NEED  TO  BE 
REPEATED  FOR  SOME  OTHER  VALUE. 

1204 

PROGRAM  LINES  1205  THROUGH  1225  CAN  BE  INVOLVED  ONLY  IF  NHTKLL  EXCEEDS  0. 

1205-1225 

N  EQUALS  1  FOR  FIRST  ROUNDS,  2  FOR  ROUNDS  IMMEDIATELY  FOLLOWING  ROUND  THAT 

HITS  TARGET ,  3  FOR  ROUNDS  IMMEDIATELY  FOLLOWING  SENSED  MISS#  4  FOR  ROUNDS 
IMMEDIATELY  AFTER  MISS  THAT  IS  MOT  SENSED,  AND  5  FOR  HITTING  ROUNDS  WHOSE 
IMPACT  POINTS  ARE  CONSIDERED  SPREAD  UNIFORMLY  AT  RANDOM  OVER  VERTICAL  TARGET 
AREA.  CALCULATE  HIT  PROBABILITIES  AND  PROBABILITIES  OF  HITTING  AND  KILLING 
TARGET  FOR  N  EQUALLING  1  THROUGH  4  AND  KILL  PROBABILITY  OF  RANDOM  HIT  FOR  N 
EQUALLING  5. 

1206-1207 

CALCULATIONS  CORRESPONDING  TO  PROGRAM  LINES  1208  THROUGH  1214  ARE  NOT 
APPLICABLE  FOR  M  EQUALLING  5.  FOP.  N  EQUAL  TO  1  THROUGH  4,  PROCESSING 
INDICATED  BY  LINES  1208  THROUGH  1224  IS  DONE  WHEN  NONZERO  STANDARD  DEVIATION 
INPUTS  HAVE  BEEN  PREVIOUSLY  ESTABLISHED  OP  IS  OTHERWISE  SKIPPED. 

1208-1214 

SIMILAR  TO  PROGRAM  LINES  1189  THROUGH  1104.  EACH  PPH(N)  VALUE  IS  CUMULATIVE 

HIT  PROBABILITY  FOP  ALL  CELLS  ALREADY  PROCESSED  AND  EVENTUALLY  APPLIES  TO 
ENTIRE  TARGET.  CONTINUE  AT  PROGRAM  STATEMENT  2490. 

1215 

FOR  NRDTYP  EQUALLING  0,  SKIP  PROGRAM  LINES  1216  THROUGH  1224  WHEN  ASSOCIATED 
CALCULATIONS  HAVE  ALREADY  BEEN  DONE  FOR  FIRST  RANGE.  PROCESSING  INVOLVED  IS 
REPEATED  FOR  ALL  OTHER  RANGES  WHEN  TARGET  VULNERABILITY  DATA  CAN  VARY  WITH 
RANGE. 

1216-1217 

FOR  N  EQUALLING  5,  PROBABILITY  OF  HITTING  EACH  CELL  IS  NOT  INVOLVED  IN  SAME 

WAY  AS  FOR  OTHER  N  VALUES  AMD  IS  SET  TO  1.  THIS  SPECIAL  SETTING  CAUSES 

PPH (5 )  TO  BECOME  CELL  COUNTER. 

1218 

NCLLS2  IS  UPDATED  AND  EVENTUALLY  EQUALS  PRODUCT  OF  1)  TOTAL  NUMBER  OF  TARGET 
CELLS  AND  2)  NUMBER  OF  SETS  OF  NONZERO  HORIZONTAL  AND  VERTICAL  STANDARD 
DEVIATIONS  USED  FOP  HIT  PROBABILITY  CALCULATIONS  PLUS  1  IF  SETTING  OF  PPOFXY 
TO  1.0  AND  PROGRAM  LINES  1220  THROUGH  1224  APE  NOT  SKIPPED  FOR  N  EQUAL  TO  5. 

1219 

SKIP  PROGRAM  LINES  1220  THROUGH  1224  IF  ONLY  HIT  PROBABILITIES  ARE  INVOLVED. 

1220-1224 

SIMILAR  TO  PROGRAM  LINES  1196  THROUGH  1200.  FOP.  N  EQUALLING  1  THROUGH  4,  EACH 
PKSHOTt  L, N)  VALUE  IS  CUMULATIVE  PROBABILITY  ACCOUNTING  FOR  ALL  CELLS  ALREADY 
PROCESSED  AND  EVENTUALLY  EQUALS  PROBABILITY  OF  HIT  AND  KILL  FOP  ENTIRE 
TARGET.  PKSHOTt L, 5)  VALUES  ARE  SIMILAR  BUT  DO  NOT  REPRESENT  PROBABILITIES. 

1228 

PROGRAM  LINES  1230  THROUGH  1234  CAN  BE  INVOLVED  ONLY  IF  NHTKLL  EXCEEDS  0. 
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EXPLANATIONS  (CONTINUED) 


122° 

SKIP  PROGRAM  LIMES  1230  THROUGH  1234  IF  ONLY  HIT  PROBABILITIES  ARE  INVOLVED. 

1230-1234 

DIVIDE  PKSHOT ( L#  M )  VALUES  BY  PPH(N)  TO  OBTAIN  KILL  PROBABILITIES  PER  HIT  ON 
TARGET.  FOR  N  EQUALLING  5#  DIVISOR  IS  TOTAL  NUMBER  OF  TARGET  CELLS  BECAUSE 
UNITY  CHANCE  OF  HIT  HAS  BEEN  PREVIOUSLY  USED  FOR  EACH  CELL. 

1236-1237 

PROGRAM  LINES  1238  THROUGH  1258  CAN  BE  INVOLVED  ONLY  IF  NSMDTR  EXCEEDS  0  AND 
RANGE  OF  CONCERN  DOES  NOT  REQUIRE  INTERPOLATION  OF  VULNERABILITY  DATA. 

1238-1258 

WRITE  OUTPUT. 

126Q 

PROGRAM  LINES  1261  THROUGH  1292  CAM  BE  INVOLVED  ONLY  IF  MHTKLL  EXCEEDS  0. 

1261-1292 

WRITE  OUTPUT. 

1261-1263 

SPECIAL  TITLING  INVOLVING  IC  AMD  MSET  APPLIES  ONLY  FOR  MOVING  TAPGET  OR  MOVING 
FIRING  WEAPON  RUM. 

1275-1278 

FOR  MOVING  TARGET  OR  MOVING  FIRING  WEAPON  RUN*  WRITE  FIRST  ROUND  HIT 

PROBABILITY  ONLY  AFTER  ADJUSTING  TO  ACCOUNT  FOR  RELIABILITY  FACTOR  RELT. 

1280 

SKIP  PROGRAM  LINES  1281  THROUGH  1292  IF  ONLY  HIT  PROBABILITIES  APE  INVOLVED. 

1285-1288 

FOR  MOVING  TARGET  OR  MOVING  FIRING  WEAPON  P.UN#  WRITE  KILL  PROBABILITIES  PER 

SHOT  ONLY  AFTER  ADJUSTING  TO  ACCOUNT  FOR  RELIABILITY  FACTORS  RELT  AMD  PELF. 

1294-1295 

WRITE  NCLLS1  AND  NCLLS2  FOP.  POSSIBLE  USE  IN  CHECKING  OUTPUT.  EITHER  QUANTITY 
STILL  EQUALLING  99999  INDICATES  ASSOCIATED  CALCULATIONS  WERE  SKIPPED. 

1297-1358 

IF  SIMULATED  ENGAGEMENTS  AP.E  INVOLVED  IN  RUM  AMD  HAVE  BEEN  COMPLETED  FOR  BOTH 
TARGET  ORIENTATION  ANGLES  OF  CONCERN#  PROCESS  DATA  IN  Z  AND  AKIL  ARRAYS  AND 
PRINT  OUTPUT. 

1299-1315 

CALCULATE  AVERAGE  VALUES  PER  SAMPLE  ENGAGEMENT  FOR  EACH  QUANTITY  IN  Z  ARRAY. 
AVERAGE  Z ( I 2# I )  VALUES  PER  ENGAGEMENT  FOR  I  EQUALLING  1  THROUGH  1C  CAN  BE 
INTERPRETED  AS  CORRESPONDING  PROBABILITIES  OF  HITTING  OR  KILLING  TAPGET. 
PERSONNEL  CASUALTIES  ASSOCIATED  WITH  I  EQUALLING  11  AND  12  REPRESENT  AVERAGE 
EXPECTED  CASUALTIES  APPLICABLE  TO  TARGET  VEHICLE  PASSENGERS. 

1317-1344 

IF  RUN  INVOLVES  SUBSEQUENT  ROUNDS  AS  WELL  AS  FIRST  ROUNDS  AND  BASIC  FIFE 

ADJUSTMENT  PROCEDURE  FOR  MISSING  ROUNDS#  CALCULATE  AND  PRINT  AVERAGE  NUMBERS 
OF  ROUNDS  NEEDED  TO  HIT  TARGET  AMD  TO  ACHIEVE  TARGET  KILL  ACCORDING  TO 
VARIOUS  CRITERIA  CONSIDERED. 

1322-1344 

ACHIEVEMENT  OF  FIRST  HIT  AND  KILLING  OF  TARGET  ACCORDING  TO  VARIOUS  CRITERIA 
ARE  CONSIDERED  FOR  RELATED  FORWARD  AMD  REVERSE  ORIENTATION  ANGLES. 

1323-1325 

CONSIDER  ONLY  FIRST  HIT  ON  TARGET  WHEN  NPRHIT  EQUALS  1. 

1327-1328 

INITIAL  SETTINGS. 

1329-1340 

AVERAGE  (OR  EXPECTED)  NUMBER  OF  ROUNDS  NEEDED  FOR  TARGET  HIT  OR  TARGET  KILL 

CAM  BE  OBTAINED  BY  ADDING  ALL  PRODUCTS  OF  1)  SPECIFIC  NUMBER  OF  ROUNCS  AND 

2)  CORRESPONDING  FRACTION  OF  ALL  ENGAGEMENTS  FOR  WHICH  EXACTLY  THAT  HUMBER 

OF  ROUNDS  ARE  REQUIRED  TO  HIT  OR  KILL  TARGET.  ONLY  FIRST  NP.DS  ROUNDS  AP.E 
PROCESSED  EXPLICITLY  IN  THIS  FASHION.  APPROXIMATE  EXTRAPOLATION  FORMULA  IS 
APPLIED  AS  NECESSARY  TO  ACCOUNT  FOR  ALL  OTHER  ROUNDS. 

1329-1332 

FOR  I  EQUAL  TO  2#  SMIK  =  SMIK  +  2 ( Z ( 2# K )-Z ( 1# K) ) .  TERM  SMIK  APPEARING  TO  RIGHT 
OF  EQUAL  SIGN  IM  THIS  EQUATION  IS  EQUIVALENT  TO  1(Z(1#K)).  SINCE  Z(2»K)  IS 
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1255 


1755 


7680 

7670 


7755 


7702 


7705 

1265 

7710 

7720 

1270 

7730 

7740 

7750 

1275 

7753 

7755 

12  80 

_ 7770 

1785 


7773 

1290 

7775 

7780 

7760 

7800 

1295 

7810 

5555 

1300  8110 


8120 

1305  8125 


1310 


DO  7670  1-1,4 

WRITE  (  6,7680  )  (  PKILL ( I, J, K),K-1, 5  ) 

FORMAT  {  10X, 21H  PKILL(4,3,5)  ,5F10.4  ) 

CONTINUE 

CONTINUE 

IF  {  NHTKLL  .EQ.  0  )  GO  TO  7800 
IF  {  ISTMOV  .NE.  1  )  GO  TO  7705 
WRITE  (  6,7702  )  IC,MSET 

FORMAT  {  /  /,  20X, 6H  IC  -  ,I3,10X,8H  MSET  -  ,13  ) 

WRITE  (  6,1020  ) 

WRITE  (  6,7710  )  (BBX  (I ),  1-1,  4  ) 

FORMAT  {  10X, 11H  BBX(4)  ,10X,4F10.4  ) 

WRITE  (  6,7720  )  ( BBY(I ), 1-1,4  ) 

FORMAT  (  10X,11H  BBY{4)  ,10X,4F10.4  ) 

WRITE  {  6,7730  )  { SSIGX  { I ),  1-1, 4  ) 

FORMAT  {  10X, 11H  SSIGX14)  ,10X,4F10.4  ) 

WRITE  {  6,7740  )  { SSIGY { I ), 1-1, 4  ) 

FORMAT  {  10X, 11H  SSIGY( 4 )  ,10X,4F10.4  ) 

WRITE  (  6,7750  )  ( PPH(I), 1-1, 5  ) 

FORMAT  (  10X, 21H  PPH<5)  ,5F10.4  ) 

IF  {  ISTMOV  .NE.  1  )  GO  TO  7755 

PPH(l)  -  PPH(l)  *  RELT 

WRITE  (  6,7753  )  (PPH(I ), 1-1, 5  ) 

FORMAT  (  10X,21H  PPH(5)  ,20X,5F10.4  ) 

CONTINUE 

IF  (  NPRHIT  .EQ.  1  )  GO  TO  7800 
DO  7760  I  »  1,4 
PKSHOT { I, 5)  -  0.0 

WRITE  (  6,7770  )  { PKSHOT { I, J ), J-l, 5  ) 

FORMAT  {  10X, 21H  PKSH0T{4,5)  ,5F10.4  ) 

IF  {  ISTMOV  .NE.  1  )  GO  TO  7775 
PKSHOT (1,1)  -  PKSH0T(I,1)  ♦  RELT  *  RELF 
WRITE  (  6,7773  )  ( PKSHOTt I, J), J-l, 5  ) 

FORMAT  {  10X,21H  PKSH0T(4,5)  ,20X,5F10.4  ) 

CONTINUE 

WRITE  {  6,7780  )  (PPKHIT(I, J ),J-1,5  ) 

FORMAT  (  10X, 21H  PPKHIT<4,5)  ,5F10.4  ) 

CONTINUE 

WRITE  (  6,1020  ) 

WRITE  (  6,7810  )  NCLlSl, NCLLS2 

FORMAT  {  10X, 18H  NCLLS1,  NCLLS2  -  ,2110  ) 

CONTINUE 

IF  {  NHTKLL  .EQ.  9  )  GO  TO  9000 
IF  (  NBETA  .LT.  NaNGLE  )  GO  TO  9000 
WRITE  (  6,8110  ) 

FORMAT  {  /  /,  11X, 28H  NRDS  PROB.  (EXP.  CAS.),  /  ) 

DO  8120  I  -  1,12 

IF  {  I  .GT.  10  .AND.  PASSN  .EQ.  0.0  )  GO  TO  8125 

DO  8120  12  -  1, NRDS 

2(12,1)  -  2(12,1)  /  SAMP 

CONTINUE 

DO  8130  N  -  1,3 

IF  (  N  .EQ.  3  .AND.  PASSN  .EQ.  0.0  )  GO  TO  8200 
WRITE  (  6,1020  ) 

IL  -  5*N  -  4 
IR  -  IL  ♦  4 

IF  (  N  .EQ.  3  )  IR  -  12 
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DO  8130  12  -  1*  NRDS 

WRITE  (  6,8140  )  12, (Z(I2, I),I*IL»IR) 

0140  FORMAT  (  11X, I4,9F13. 7  ) 

1315  8130  CONTINUE 

0200  CONTINUE 

-  IF  (  NRD1  .EQ.  1  )  60  TO  9000 

IF  (  NADJST  .GT.  0  )  GO  TO  8400 
WRITE  (  6,8310  ) 

1320  0310  FORMAT  (  /  /,  31X,  17H  AVG.  NO.  OF  RDS.»  /  ) 

WRITE  (  6,1020  ) 

DO  8320  K  -  1,10 
IF  (  NPRHIT  .EQ.  0  )  GO  TO  8325 
IF  (  K  .EQ.  5  .OR.  K  .EQ.  10  )  GO  TO  0325 
1325  GO  TO  8320 

0325  CONTINUE 

-  SMIK  -  Z { 1, K ) 


AUL  -  2.0 

- -  DO  8330  I  -  2, NRDS 

1330  SMIK  -  SMIK  +  AUL* ( Z( I, K )-Z ( 1-1, K ) ) 

0330  AUL  -  AUL  ♦  1.0 

IF  (  Z (NRDS, K)  .LE.  0.9999  )  GO  TO  0335 

-  AV  -  SMIK 

GO  TO  8340 

1335  8335  AR  -  5.0  -  ( Z ( NRDS-5, K) +Z (NRDS-4, K )  +  Z (NRDS-3, K  )  +Z (NRDS-2,  K 

A  *Z(NRDS-1,K)  ) 

XX  -  Z ( NRDS, K)  -  Z (NRDS-5, K ) 

IF  (  XX  .LE.  0.0  )  GO  TO  0335 
AZ  -  AR  /  XX 

1340  AV  -  SMIK  ♦  (AZ+RDS)*(1.0-Z(NRDS,K) ) 

0T4O  IF  (  AV  .GT.  999.0  )  AV  -  999.0 

-  WRITE  (  6,8350  )  K,AV 

0350  FORMAT  (  20X, 15, 1  OX, F7.  2  ) 

0320  CONTINUE 
1345  8400  CONTINUE 

“  WRITE  (  6,8402  > 

0402  FORMAT  (  /  /,  22X, 24H  PROBABILITY  VERSUS  TIME,  /  ) 

DO  8410  J  »  1,10 
IF  (  NPRHIT  .EQ.  0  )  GO  TO  8415 
1350  IF  (  J  .EQ.  5  .OR.  J  .EQ.  10  )  GO  TO  8415 

GO  TO  8410 

8415  WRITE  (  6,1020  ) 

DO  8420  K  »  1,61 

8420  AK IL (K, J )  -  AKIL(K,J)  /  SAMP 
1355  WRITE  (  6,8430)  (  AKIL ( K, J ) ,K*1, 61  ) 

0430  FORMAT  (  11X,10F8.5  ) 

8410  CONTINUE 

WRITE  (  6,1010  ) 

9UIT5  CONTINUE 

1565 -  IF  (  ISTMOV  .EQ.  1  )  GO  TO  9305 

90T0  IF  (  ISPLIT  .EQ.  4  )  ISPLIT  •  0 

-  IF  (  ISPLIT  .NE.  2  )  GO  TO  9100 

NANGLE  -  NANGLE  -  1 
IF  (  NRANGE  .EQ.  1  )  NEJECT  -  1 
1365  GO  TO  3000 

9100  CONTINUE 

-  IF  (  NSTCRD  .GT.  0  )  GO  TO  9302 

- -  NANGLE  *  NANGLE  +  1 


) 
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EXPLANATIONS  (CONTINUED) 


1333-1340 


1341 

1346-1357 

1360 

1361 

1362-1365 

136? 

1368-1370 

13'*1-1374 

1375-13^9 


FRACTION  OF  ALL  ENGAGEMENTS  I'M  WHICH  TARGET  WAS  HIT  OR  KILLED  BY  FIRST  0P. 
SECOND  ROUND*  DIFFERENCE  Z ( 2* K)-Z (1*K )  REPRESENTS  FRACTION  OF  ALL 
ENGAGEMENTS  IN  WHICH  ACHIEVING  HIT  OR  KILL  REQUIRED  EXACTLY  TWO  ROUNDS.  FOR 
I  EQUAL  TO  3*  SMIK  EQUALS  ITS  PREVIOUS  VALUE  PLUS  TERM  3 ( Z ( 3, K ) -Z ( 2, K ) ) . 
EVENTUALLY  SMIK  ACCOUNTS  FOR  ALL  OF  FIRST  NRDS  ROUNDS.  FURTHER  CALCULATING 
IS  UNNEEDED  IF  Z(NRD$*K)  IS  EQUAL  TO  1  OR  VERY  NEARLY  SO*  INDICATING  THAT 
NRDS  OR  FEWER  ROUNDS  SUFFICED  IN  ALL  OR  PRACTICALLY  ALL  SAMPLE  ENGAGEMENTS. 
IN  THAT  CASE  AVERAGE  NUMBER  OF  ROUNDS  AV  EQUALS  SMIK. 

APPLY  APPROXIMATE  AND  RATHER  INTRICATE  EXTRAPOLATION  FORMULA  TO  ACCOUNT  FOR 
ENGAGEMENTS  WHERE  NRDS  OR  FEWER  ROUNDS  DID  NOT  SUFFICE  TO  HIT  OR  KILL 
TARGET.  FORMULA  IS  BASED  ON  CONDITIONAL  PROBABILITIES  OF  HIT  OR  CONDITIONAL 
PROBABILITIES  OF  KILL.  SUCH  CONDITIONAL  PROBABILITIES  HAVE  NOT  BEEN  INVOLVED 
SO  FAR  IN  PROGRAM  CALCULATIONS  BUT  ARE  INFERRABLE  FROM  PROBABILITIES  IN  Z 
ARRAY.  IF  NRDS  IS  NOT  TOO  SMALL*  CONDITIONAL  PROBABILITIES  FOR  EACH  ROUND 
FOLLOWING  FIRST  NRDS  ROUNDS  CAN  BE  CONSIDERED  APPROXIMATELY  EQUAL  TO  AVERAGE 
OF  FIVE  CONDITIONAL  PROBABILITIES  APPLYING  TO  LAST  FIVE  OF  FIRST  NRDS 
ROUNDS.  LOWER  BOUND  OF  10  FOR  NRDS  HAS  BEEN  CHOSEN  AS  ACCEPTABLE  IN  THIS 
CONTEXT.  AZ  REPRESENTS  AVERAGE  NUMBER  OF  ROUNDS  NEEDED  IN  ADDITION  TO  FIRST 
NRDS  ROUNDS  WHEN  NRDS  ROUNDS  ARE  INSUFFICIENT.  CORRESPONDING  FRACTION  OF 
ENGAGEMENTS  IS  INDICATED  BY  DIFFERENCE  1.0-Z (NRDS, K) .  AVERAGE  NUMBER  OF 
ROUNDS  AV  IS  OBTAINED  BY  ADDING  SMIK  VALUE  ASSOCIATED  WITH  ENGAGEMENTS  WHERE 
NRDS  ROUNDS  SUFFICE  AND  PRODUCT  ACCOUNTING  FOR  ALL  OTHER  ENGAGEMENTS. 

PRIMARILY  TO  AVOID  HAVING  TO  PRINT  VERY  LARGE  NUMBERS*  REPLACE  ANY  NUMBER 
LARGER  THAN  999  BY  THIS  LIMIT.  THERE  IS  OF  COURSE  NO  PRACTICAL  INTEREST  IN 
MANY  AVERAGE  NUMBERS  OF  ROUNDS  THAT  ARE  MUCH  SMALLER  THAN  999  BUT  STILL  TOO 
LARGE  TO  CORRESPOND  TO  REALISTIC  COMBAT  ENGAGEMENT  SITUATIONS. 

CALCULATE  AVERAGE  VALUES  PER  SAMPLE  ENGAGEMENT  FOR  EACH  QUANTITY  IN  AKIL 
ARRAY.  CONSIDER  ONLY  FIRST  HIT  ON  TARGET  WHEN  N PR HIT  EQUALS  1.  RESULTING 
AVERAGE  AKIL ( K*  J )  VALUES  CAN  BE  INTERPRETED  AS  CORRESPONDING  PROBABILITIES 
OF  HITTING  OR  KILLING  TARGET  BEFORE  2(K-1)  SECONDS. 

SKIP  PROGRAM  LINES  1361  THROUGH  136°  FOR  MOVING  TARGET  OR  MOVING  FIRING  WEAPON 
RUN. 

SINCE  CALCULATIONS  INVOLVING  RANGE  INTERPOLATION  OF  TARGET  VULNERABILITY  DATA 
HAVE  BEEN  COMPLETED*  RESET  ISPLIT  TO  INITIAL  SETTING. 

PROGRAM  LINES  1363  THROUGH  1365  APPLY  ONLY  IF  ISPLIT  IS  2.  CALCULATIONS  FOR 
FIRST  PART  OF  SIMULATIONS  HAVE  JUST  BEEN  COMPLETED  FOR  PEVERSE  ANGLE  OF 
CONCERN  AND  MANGLE  NEEDS  TO  BE  RESET  BEFORE  SECOND  PART  OF  SIMULATIONS  FOR 
CORRESPONDING  FORWARD  ANGLE.  NEJECT  IS  ASSOCIATED  WITH  LINE  SKIPPING 
OCCASIONALLY  NEEDED  WHEN  TARGET  VULNERABILITY  DATA  ARE  READ  FROM  TAPE  OR 
DISC.  CONTINUE  AT  PROGRAM  STATEMENT  3000. 

SKIP  PROGRAM  LIMES  1368  THROUGH  1370  IF  RUN  INVOLVES  TARGET  SHAPE  DATA  FROM 
CARDS . 

CONTINUE  EITHER  AT  PROGRAM  STATEMENT  4020  FOP.  NEXT  VALUE  OF  MANGLE  IF  ANY  OR 
AT  STATEMENT  9310. 

IF  PUN  INVOLVES  TARGET  SHAPE  DATA  FROM  CARDS*  CONTINUE  AT  PROGRAM  STATEMENT 
4210  FOR  NEXT  NTRGT C  VALUE  OF  CONCERN  OR  AT  STATEMENT  °310. 

APPLICABLE  FOR  MOVING  TARGET  OR  MOVING  FIRING  WEAPON  RUN  ONLY.  MSET  EQUALS  1 
THROUGH  6  FOR  VARIOUS  TARGET  SPEED  AND  EVASIVE  MANEUVERING  CONDITIONS. 

NANGLE  VALUES  FOR  MOVING  TARGET  OR  MOVING  FIRING  WEAPON  RUN  ARE  RESTRICTED 
TO  1  FOR  0  DEGREES*  3  FOR  30  DEGREES*  AND  5  FOR  60  DEGREES.  CONTINUE  AT 
PROGRAM  STATEMENT  4035  IF  ALL  SIX  MSET  VALUES  HAVE  BEEN  CONSIDERED. 
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PROGRAM  LISTING  (CONTINUED) 


IF  {  NANGLE  .LE.  8  )  GO  TO  4020 

1370 _  GO  TO  9310 

9302  NTRGTC  »  NTRGTC  +  1 

IF  (  NTRGTC  .GT.  4  )  GO  TO  9310 

IF  {  KTRGTC (NTRGTC)  .GT.  0  )  GO  TO  4210 

_  GO  TO  9310 

1375  9305  MSET  -  MSET  +  1 

IF  {  MSET  .LE.  6  )  GO  TO  4035 
MSET  -  1 

NANGLE  -  NANGLE  +  2 
IF  {  NANGLE  .LE.  5  )  GO  TO  4020 
1380  TTF5  NRANGE  -  NRANGE  +  1 

IF  (  NRDTYP  .EQ.  1  . AND. '  NE JECT  .EQ.  0  )  NEJECT  -  1 
IF  (  NRANGE  .LE.  NCASES  )  GO  TO  2010 
9900  CONTINUE 
CALL  EXIT 

1385  END 


1  SUBROUTINE  NRAN{R1*R2) 

A  -  SORT (—2 • O+ALOG ( RANF IX  ) )  ) 
B  -  6.28318530718  *  RANF(X) 
R1  »  A  *  SIN (B ) 

5  R2  *  A  *  COS { B ) 

RETURN 

END 


1  FUNCTION  CNORM(X) 

F  -  0.0 
AX  -  ABS(X) 

IF  (  AX  .GE.  5.0  )  GO  TO  10 

5  F  -  {(U(0.5383E-5*AX+0.488906E-4)*AX+0.380036E-4)*AX  +  0. 0032776263 

A  )* AX+O. 0211410061 )* AX +0.0498 673 4 69 ) *AX  +  1.0 
F  -  .5  /  {(F**8)**2) 

10  IF  {  X  .GE.  0.0  )  F  -  1.0  -  F 
CNORM  -  F 

10  RETURN 

END 
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EXPLANATIONS  (CONTINUED) 


OTHERWISE  RESET  MSET  TO  INITIAL  SETTING  BEFORE  CONTINUING  AT  STATEMENT  4020 
OR  9310  ACCORDING  TO  WHETHER  CALCULATIONS  HAVE  OR  HAVE  NOT  YET  BEEN  DONE  FOR 
ALL  THREE  MANGLE  VALUES  OF  CONCERN. 

1380-1382  CONTINUE  AT  PROGRAM  STATEMENT  2010  UNLESS  CALCULATIONS  FOR  LAST  RANGE  HAVE 

JUST  BEEN  COMPLETED.  NEJECT  MAY  MEED  TO  BE  RESET  WHEN  TARGET  VULNERABILITY 
DATA  FROM  TAPE  OR  DISC  ARE  INVOLVED. 

1383-1385  RUN  STOPS*  POSSIBLY  PREMATURELY  IF  PROBLEM  HAS  OCCURRED. 


SUBROUTINE  MRAN{R1*R2 )  “RODUCES  PAIR  OF  RANDOM  NORMAL  DEVIATES  R1  AND  R2  FROM 
NORMAL  DISTRIBUTION  WITH  MEAN  OF  0  AND  STANDARD  DEVIATION  OF  1. 


WHEN  X  EQUALS  (X2-BX  (N)  ) /SIGX  (K*M  )*  FOP.  EXAMPLE*  FUNCTION  CNOP.M(X)  CALCULATES 
PROBABILITY  OF  RANDOM  VARIABLE  LYING  BETWEEN  MINUS  INFINITY  AND  X2  FOR 
NORMAL  DISTRIBUTION  WITH  BIAS  OF  BX(N)  AMD  STANDARD  DEVIATION  OF  SIGX(K*N). 
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APPENDIX  B 


EXPLANATIONS  OF  PROGRAM  STATEMENTS  919  THROUGH  1131 
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EXPLANATIONS  OF  PROGRAM  STATEMENTS  919  THROUGH  1131 


This  appendix  contains  a  paper  presented,  In  early  1979,  at  the 
Am\y  Numerical  Analysis  and  Computers  Conference.  This  paper,  entitled 
"A  FORTRAN  ROUTINE  FOR  ESTIMATING  NORMAL  DISTRIBUTION  PARAMETERS",  pro¬ 
vides  detailed  explanations  of  program  lines  919  through  1131  of  the 
direct  fire  program. 


The  paper  was  subdivided  Into  three  parts,  namely,  the  Intro¬ 
duction,  the  description  of  the  basic  mathematical  procedure,  and  the 
conclusion.  The  contents  of  these  three  parts  follow. 


1.  INTRODUCTION.  A  Monte  Carlo  computer  program  simulating  the  engagement 
of  a  single  target  by  a  t3nk  main  armament  system  has  been  developed  by  the 
Joint  Munitions  Effectiveness  Manual  (JMEM)  Methodology  and  Evaluations  Working 
Group,  a  tri-service  group  responsible  for  establishing  certain  standardized 
estimates  of  weapon  effectiveness.  The  program  is  often  referred  to  as  the 
JMEM  Direct  Fire  program.  Exercise  of  this  program  yields  large  amounts  of 
data  concerning  the  location  of  round  impacts  in  the  target  plane.  First 
rounds  fired  against  a  particular  target  are  of  chief  concern,  and  these  can 
be  subdivided  into  rounds  that  hit  and  rounds  that  miss  the  target.  Chart  1 
illustrates  a  possible  target  and  a  few  conceivable  first  round  impact  points. 


PaPer  describes  a  computer  routine  for  processing  the  horizontal 
coordinates  or  the  vertical  coordinates  of  hitting  rounds  or  of  missing 
rounds.  The  routine  provides,  for  each  set  of  data  considered,  an  estimate 
of  the  mean  and  three  estimates  of  the  standard  deviation  for  a  normal 
distribution  tentatively  assumed  to  fit  the  data.  After  the  computer  run 
an  analyst  can  judge  by  comparing  the  three  standard  deviation  values  whether 
the  tentative  assumption  of  normality  is  sufficiently  substantiated  and  if  so 
make  a  best  estimate  for  the  parameters  of  interest.  The  word  "routine"  is  used 
for  the  logical  processing  documented  in  this  paper.  However,  the  associated 
program  statements  have  not  actually  been  structured  as  a  separate  routine  but 
are  a  portion  of  the  complete  engagement  simulation  program  previously  mentioned, 
separate  routine  could  readily  be  developed  for  other  applications. 


.  1: _ DESCRIPTION  OF  ROUTINE.  The  program  instructions  of  interest  are  listed 

m  Chart  2.  Note  that  the  entire  chart  consists  ,of  the  loop  DO  7010  I  =  1,2 
The  index  I  equals  1  for  horizontal  coordinates  and  2  for  vertical  coordinates. 
For  any  other  application,  one  could  allow  for  more  values  of  I  or  establish 
a  single  dummy  setting  for  this  index. 
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CHART  1  IMPACT  POINTS  FOR  HITTING  AND  MISSING 
FIRST  ROUNDS 
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CHART  2  PROGRAM  STATEMENTS 


DO  7010  I  «  1,2 

AVRG  -  0.0 

TMINUS  «  0.0 

TPLUS  «  0.0 

KTIMES  «  0 

DO  7020  J  «  1,50 

TERM1  «  NHINUS( J, I, N) 

TER  M2  «  NPLUS ( J,I,N) 

TMINUS  «  TMINUS  *  TERM1 

7020  TPLUS  «  TPLUS  ♦  TERM2 

TOTAL(I)  «  TMINUS  ♦  TPLUS 
NTOTAL  «  TOTAL ( I ) 

T5000  «  TOTAL(I)  *  0.5 
N5000  »  T5000 

IF  (  T5000  ,NE.  TMINUS  )  GO  TO  7024 
DO  7022  K  «  1,50 
NPOS(K)  -  NPLUS  <K, I, N ) 

7022  NNLG(K)  «=  NMI  NUS  (  K,  I ,  N) 

GO  TO  7200 

7024  IF  (  T5000  ,iLT.  TMINUS  )  GO  TO  7025 
DIFF  -  T5000  -  TMINUS 
DO  7030  J  «  1,50 
PLUS  «  NPLUS ( J, I, N) 

IF  (  PLUS  .IT.  DIFF  )  GO  TO  7035 

IFRCTN  -  1000.0  *  DIFF  /  PLUS 

FRCTN  -  FLOATCIFP.CTN)  /  1000.0 

AVRG  «  AVPG  +  FP.CTN*20. 0 

JEND  «  100  -  KTIMES 

LTIMES  *  50  -  KTIMES 

LEND  «  1  +  KTIMES 

DO  7040  K  «  1,50 

NNEG(K)  «  0 

7040  NPOSCK)  *  0 
NCMLTN  *  0 
NCMLTP  «  0 
L  «  0 

7050  IF  (  L  .LQ.  100  )  GO  TO  7200 
L  *  L  +  1 
LL  -  L 

IF  (  L  .LE.  50  )  LL  «  51  -  L 
M  «  L 

IF  (  L  .GT.  50  )  M  «  L  -  50 
IF  (  LL  .LE.  JEND  )  GOTO  7055 
NPOS(M)  «  0 
GO  TO  7060 

7055  IF  (  LL  .LT.  JEND  )  GOTO  7065 

NTERM  -  (IFRCTN* NPLUS(50, I, N)+500)  /  1000 
NPOS(M)  «  NPLUS (50, I, N)  -  NTERM 
NCMLTP  -  NCMLTP  ♦  NPOS(M) 

GO  TO  7060 

7065  IF  (  LL  .LE.  LTIMES  )  GOTO  7075 
IF  (  LL  .LE.  50  )  GOTO  7070 
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CHART  2  PROGRAM  STATEMENTS  (CONTINUED) 

NTERMl  ■  <IFRCTN*NPLUS(l-LTIMES*I*N)+500)  /  1000 
NTERM2  ■  <IFRCTN*NPlU$(L-LTIHES+l>I*N)+500)  /  1000 
NPOS(M)  -  NPLUS ( L-L TIMES* I#  N)  -  NTERMl  ♦  NTERM2 
NCMLTP  •  MCMLTP  +  NPOS(M) 

GO  TO  7060 

7070  NTERMl  «  ( I F RCTN*NPLIJS(  LEND-M*  I *  N  ) ♦  500 )  /  1000 

NTERM2  *  <IFRCTN*NPLUS(LEND-M+lrI*N)+500)  /  1000 
NNEG(M)  -  NPLUS ( LEND-M* I*  N)  -  NTERMl  +  NTERM2 
NCMLTN  *  HCMLTN  ♦  NNEG(M) 

GO  TO  7060 

7075  IF  (  LL  .£T.  LTIMES  )  GOTO  7085 

NTERMl  -  <IFRCTN*NMINUS(l*I*N)4-500)  /  1000 
NTERM2  =  (IFRCTN*NPLUS(l*I*N)+500)  /  1000 
NNEG(M)  =  NMINUS<1*I* N)  -  NTERMl  +  NTERM2 
NCMLTN  *  NCMLTN  +  NNEG(M) 

GO  TO  7060 

7085  IF  (  LL  .EQ.  1  )  GO  TO  7095 

NTERMl  =  (IFRCTN*NMlNUS(LTIME$-LL+i*I*N)+500)  /  1000 
NTEP.M2  -  (IFRCTN*NMINUS(LTIMES-LL*I*N)+500)  /  1000 
NNEG( L-KTIMES+1 )  »  NMINUS ( LTIMES-LL+1* I*N )  -  NTERMl  ♦  NTERM2 
NCMLTN  *  NCMLTN  ♦  NNEG( L-KTIMES+1 ) 

GO  TO  7060 

7095  NNCG( 50)  «  ( I FRCTN*NMINUS ( 49-KTI KES* I * N )♦ 500 )  /  1000 
NXTRA  *  KTIMES  +  1 
DO  7100  JXTRA  =  1*NXTRA 

7100  NNEG( 50)  *  NNEG ( 50)  +  NMINUSt 49-KTI MES+ JXTRA# I» N ) 

NTERM  =  (IFRCTN*NMINUS(50*I*N)+500)  /  1000 
NNEG( 50)  =  NNEG ( 50)  +  NTERM 
NCMLTN  •  NCMLTN  *  NNEG( 50 ) 

7060  IF  (  NCMLTP  «£Q.  NTOTAL  )  L  =  100 
IF  (  L  .GT.  50  )  GO  TO  7050 
IF  (  NCMLTN  .IT.  N5000  )  GO  TO  7050 
NCMLTP  =  NCMLTN 
L  «  50 
GO  TO  7050 

7035  DIFF  *  DIFF  -  PLUS 
KTIMES  -  KTIMES  +  1 
AVRG  *  AVRG  ♦  20.0 

7030  CONTINUE 

7025  DIFF  =  TMINUS  -  T5000 
DO  7110  J  «  1*50 
AMINUS  =  NMINUS ( J  *  I #N ) 

IF  (  AMINUS  .IT.  DIFF  )  GO  TO  7115 

IFRCTN  =  1000.0  *  DIFF  /  AMINUS 

FRCTN  •  FLOAT  ( IFRCTN)  /  1000.(0 

AVRG  *  AVRG  -  FRCTN+20.0 

JEND  =  100  -  KTIMES 

LTIMES  -  50  -  KTIMES 

LEND  *  1  +  KTIMES 

DO  7120  K  *  1*50 

NPOS(K)  »  0 

7120  NNEG(K)  «  0 
NCMLTP  =  0 
NCMLTN  «  0 
L  *  0 

7130  IF  (  L  .EQ.  100  )  GO  TO  7200 
L  *  L  +  1 
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CHART  2  PROGRAM  STATEMENTS  (CONTINUED) 

LL  «  L 

IF  <  l  .LE.  50  )  LL  ■  51  -  L 
M  *  L 

IF  (  L  .GT.  50  )  M  «  L  -  50 
IF  (  LL  ,ILE.  JEND  )  GOTO  7135 
NNEG(M)  -  0 
GO  TO  7140 

7135  IF  (  LL  .LT.  JEND  )  GOTO  7145 

NTERM  «  ( IFRCTN*NMINUS< 50,I,N)+500)  /  1000 
NNEG(M)  -  NMINUS< 50,I,N)  -  NTERM 
NCMLTN  -  NCMLTN  ♦  NNEG(M) 

GO  TO  7140 

7145  IF  (  LL  .IE.  LTIMES  )  GOTO  7155 
IF  (  LL  • LE*  50  )  GOTO  7150 

NTERM1  -  <IFRCTN*NMINUS(L-LTIMES»I»N)+500)  /  1000 
. NTERM2  =  <IFRCTN*NMINUS(L-LTIMES+l,I,N)+500)  /  1000 
NNEG(M)  -  NMINUSl L-LTIMES. I. N )  -  NTERH1  ♦  NTERM2 
NCMLTN  «  NCMLTN  +  NNEG(M) 

GO  TO  7140 

7150  NTERM1  =  ( I FRCTN*NMINUS ( L END-M, I, N)  +5 00 )  /  1000 

NTERM2  =  <IFRCTN*NMINUS(LEND-M+1, I.NJ+500)  /  1000 
NPOS(M)  -  NMINUS(LEND-M,I,N)  -  NTERM1  +  NTERM2 
NCMLTP  *  NCMLTP  ♦  NPOS(M) 

GO  TO  7140 

7155  IF  (  LL  •  LT •  LTIMES  )  GOTO  7165 

NTERM1  ■  (IFRCTN*NPLUS(l>I/N)+500)  /  1000 
NTERM2  =  (IFRCTN*NMINUS(l*I»N)+500)  /  1000 
NPOS(M)  -  NPLUS (1> I»  N)  -  NTEPM1  +  NTERH2 
NCMLTP  =  NCMLTP  ♦  NPOS(M) 

GO  TO  7140 

7165  IF  (  Ll  •  EQ«  1  )  GO  TO  7175 

NTERM1  -  <IFRCTN*NPLUS(LTIMES-LL+1» I»N)+500)  /  1000 
NTERM2  «  <IFRCTN*NPLUS(LTIME$-LL,I,N)*500)  /  1000 
NP0S(L-KTIMES«-1)  *  NPLUS( LTIMES-LL+1. I#N)  -  NTERM1  ♦  NTERM2 
NCMLTP  -  NCMLTP  +  NPOS( L-KTIMES+1 ) 

GO  TO  7140 

7175  N POS( 50 )  ■  ( IFRCTN*NPLUS(49-KTIMES» I#N)  +  500)  /  1000 
NXTRA  -  KTIMES  +  1 
DO  7160  JXTRA  *  1 t NXTRA 

7180  NPOS ( 50)  -  NPOS ( 5 0)  +  NPLUS (49-KTIMES+ JXTRA. I.  M ) 

NTERM  *  (IFRCTN*NPLU$(50» I#N)+500)  /  1000 
NPOS (50)  =  NP  OS ( 50)  +  NTERM 
NCMLTP  «  NCMLTP  ♦  NP0S(50) 

7140  IF  (  NCMLTN  • EQ*  NTOTAL  )  L  -  100 
IF  (  L  *GT«  50  )  GO  TO  7130 
IF  (  NCMLTP  .LT.  N5000  )  GO  TO  7130 
NCMLTN  -  NCMLTP 
L  «  50 
GO  TO  7130 

7115  DIT F  =  DIFF  -  AMINUS 
KTIMES  «  KTIMES  «•  1 
AVRG  «  AVRG  -  20.0 

7110  CONTINUE 

7200  CONTINUE 

DO  7310  K  «  1>50 

IF  (  K  .GT.  I  )  GO  TO  7315 

NSUM1  «  NPDS(K) 
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CHART  2  PROGRAM  STATEMENTS  (CONTINUED) 


NSUM2  -  NNEG(K) 

GO  TO  7320 

7315  NSUM1  »  N3UM1  ♦  NPOS(K) 

NSUM2  •  NSUM2  ♦  NNEG(K) 

7320  NCMPOS(K)  »  NSUM1 
NCMNEG(K)  »  NSUH2 
NRFLCT(K)  »  NSUM1  ♦  NSUM2 

7310  CONTINUE 

DC  7330  J  »  1* 50 
SIGMAX  -  J  *  20 
IF  (  J  .GT.  1  )  GO  TO  7335 
SMFRQl  »  0,0 

SMFRQ2  »  NRFLCT(l)  *  10000  /  NTOTAL 
K5  ■  5 
K 10  »  10 
GO  TO  7350 

7335  SMFRQl  «  SMFRQ2 

IF  <  J  *EQ*  5 C  )  GO  TO  7345 
SMFRQ2  *  NRFLCT  < J )  *  10C00  /  NTOTAL 
GO  TO  7350 

7345  5MFRQ2  =*  10000*0 

7350  FRQNCY  *  5MFRQ2  -  SMFRQl 

IF  (  K5  • NE*  5  )  GO  TC  7355 
IF  <  SMFRQ2  *LT.  3829*2  )  GOTO  7330 
D05  »  (SMFRQ2-3829*2)  /  FRQNCY 
SIG05  «  (SIGMAX- 005*20*0)  /  0*5 
K5  ■  0 

7355  IF  (  K10  *N£*  10  )  GO  TC  7365 

IF  (  SMFRQ2  *LT.  6026*8  )  GCTO  7330 
DIO  »  (SMFRQ2-6826.8)  /  FRQNCY 
SIG10  -  SIGMAX  -  010*20. C 


K10  »  0 

7365  IF  (  SMFRQ2  *LT*  8663*8  )  GOTO  7330 
015  »  ( SMFRQ2-8663.8)  /  FRQNCY 
SIG15  «  (SIGMAX-D15*20. 0)  /  1.5 
GO  TO  7370 
7330  CONTINUE 
7370  WRITE  (  6*102C  ) 

WRITE  (  6*7380)  N*  I*  SIGC5*SIG10*SIG15 


7380  FORMAT  (  10X*26H  N* I*  SIG05*  SIG10#  SIG15  >2I10,3F1C*4  ) 

BIASXY(I*N)  «  AVRG 


7010 


IF  (  I  «EQ*  1 
IF  (  I  *EQ*  2 
SIGMXY ( 1* I*N ) 
S IGMXY ( 2*1*N ) 
S IGMXY ( 3* I*N ) 
CONTINUE 


)  SSQRRE  *  SIGXL**2  ♦  SIGXR**2 
)  SSQRRE  ■  SIGYl**2  +  SIGYR**2 
«  SQRT(SI€05**2*SSQRRE) 

»  SORT (SIG10**2+SSQRRE) 

•  SORT (SIG15**2+SSQRRE) 
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The  index  N  is  set  before  the  DO  7010  loop  and  presently  has  possible 
values  of  1  through  5.  N  equals  1  when  hitting  first  rounds  are  of  concern 
and  3  for  missing  first  rounds.  Meanings  of  the  values  2,  4,  and  5  need  not 
be  explained  here.  The  range  of  values  of  N  can,  for  other  applications,  be 
increased  ordecreased;  in  particular,  a  single  dummy  setting  can  be  used. 

The  horizontal  and  vertical  coordinate  axes  are  each  subdivided  into  100 
intervals  of  20  inches  each;  positive  and  negative  coordinates  are  each  covered  bv 
50  intervals.  Intervals  are  related  to  the  index  J  in  a  way  that  should  soon 
be  clear.  Assume  that  one  has  run  many  simulated  engagements  and  determined, 
for  hitting  rounds  and  missing  rounds  separately,  how  many  times  a  particular 
interval  contained  the  horizontal  coordinate  of  the  first  round  and,  again 
separately,  the  vertical  coordinate.  Input  data  of  this  sort  are  shown  in 
Chart  3.  Numbers  in  the  array  NPLUS(J,I,N)  indicate  how  many  coordinates  are 
within  the  interval  0  to  20  inches  f or  J  =  1 ,  within  the  interval  20  to  40 
inches  for  J  =  2,  and  so  forth.  Similarly,  the  NMINUS(J, I ,N)  array  contains  the 
number  of  coordinates  within  the  interval  -20  to  0  inches  for  J  =  1,  and  so  forth. 

The  arrays  involved  in  Chart  2  and  their  dimensions  in  the  JMEM  Direct  Fire 
program  are  as  follows: 

NMINUS (50 ,2,5) 

NPLUS(50,2,5) 

TOTAL (2) 

NNEG(50) 

NPOS (50) 

NCMPOS(SO) 

NCMNEG(SO) 

NRFLCT (50) 

The  NMINUS  and  NPLUS  arrays  contain  the  data  to  be  processed. 

The  complete  calculations  for  particular  values  of  N  and  I  yield  an  estimate 
of  the  mean,  denoted  by  AVRG,  and  three  estimates  of  the  standard  deviation, 
represented  by  SIG05,  SIG10,  and  SIG15. 

A  total  of  1000  engagements  were  simulated  to  obtain  the  data  in  Chart  3. 

A  first  round  hit  was  obtained  in  529  instances,  and  the  first  round  missed  the 
target  in  the  other  471  instances.  The  impact  points  for  the  529  hits  were 
(for  an  observer  at  the  firing  weapon)  to  the  right  of  the  Y-axis  illustrated 
in  Chart  1  on  324  occasions,  and  to  the  left  on  the  other  205  occasions.  These 
same  impact  points  were  above  the  X-axis  198  times  and  below  331  times.  Impact 
points  of  the  471  misses  had  positive  horizontal  coordinates  in  278  engagements, 
and  positive  vertical  coordinates  in  369  engagements.  .To  illustrate  the  remain¬ 
der  of  these  explanations,  an  arbitrary  choice  has  been  made  of  the  values  N  =  3 
(missing  first  rounds)  and  1=2  (vertical  coordinates) . 

It  is  useful  to  digress  and  consider  how  the  mean  and  standard  deviation  of 
a  normal  distribution  corresponding  to  the  input  data  selected  can  be'  estimated 
graphically.  Chart  4  shows  the  associated  quantitative  basis  and  Chart  5  the 
plotted  points  to  which  one  would  attempt,  using  judgement  rather  than  calculation, 
to  fit  a  straight  line.  The  computer  routine  being  described  performs  computations 
that  basically  parallel  the  graphical  approach. 
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INTERVAL  (INCHES) 


CHART  3  SAMPLE  INPUT  DATA 
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N  =  1 

N  =  1 

N  =  3 

N  =  3 

1_=JL 

1  =  2 

J  =1 

1  =  2 

(205) 

(331) 

(193) 

(102) 

0 

0 

0 

0 

«  0 

0 

0 

0 

.  0 
»— 1 

0 
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3 
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—> 
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0 

14 

2  ii 

17 

12 

29 

§  34 

SI 

50 

21 

19 

z  65 

118 

59 

17 

93 

145 

101 

20 

0-20 

140 

109 

96 

71 

20-40 

100 

71 

81 

104 

ETC 

*  55 

»— HI 

17 

54 

104 

.  21 
—> 

1 

35 

58 

~  8 

0 

8 

19 

si  0 

0 

2 

11 

0 

0 

2 

2 

0 

0 

0 

0 

(324) 

(198) 

(278) 

(369) 

(529) 

(529) 

(471) 

(471) 
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CHART  A  BASIS  FOR  CHART  5 


N=3,  1=2 

CUMULATIVE 

CUMULATIVE  FRACTION 

SUM  (SUM/471) 


6 

3 

3 

.006 

5 

14 

17 

.04 

A 

2: 

■% 

29 

46 

,10 

3 

> 

19 

65 

.14 

2 

“O 

CO 

17 

82 

,17 

1 

13 

2: 

20 

102 

.22 

—> 

1 

§2 

71 

173 

.37 

2 

104 

277 

.59 

3 

104 

381 

.81 

4 

■v 

>— • 4 

58 

439 

,93 

5 

3 

19 

458 

,97 

6 

CO 

3D 

11 

469 

.996 

7 

_ I 

D- 

z: 

2 

471 

1.000 
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PROBABILITY 
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CHART  5  DISTRIBUTION  OF  COORDINATE  DATA 
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The  steps  involved  in  estimating  the  mean  are  illustrated  in  Chart  6.  Note 
that,  for  this  example,  the  nonzero  input  data  of  concern  are  NPLUS(1,I,N)  =  71, 
NPLUS (2 , I ,N)  =  104,  .  .  .,  NPLUS (7 , I ,N)  =  2,  NMINUS (1 , I ,N)  =  20,  .  .  .,  NMINUS 
(6,1 , N)  *  3.  The  totals  102  and  369  are  the  values  of  TMINUS  and  TPLUS  after 
completion  of  the  DO  7020  loop.  Since  TPLUS  exceeds  TMINUS,  the  program  attempts 
to  identify  a  positive  value  on  the  Y-axis  as  the  estimate  of  the  mean.  Conse¬ 
quently,  DIFF  is  calculated  by  the  statement  right  before  the  DO  7030  loop, 
rather  than  by  the  statement  preceding  the  DO  7110  loop.  The  DO  7030  loop 
establishes  first  that  the  point  of  central  tendency  is  at  least  as  great  as  20 
inches,  the  right  end  of  the  first  interval  considered.  This  is  so  because  133 
impact  point  coordinates  need  to  be  dropped  from  the  NPLUS  group  before  the 
coordinates  equal  half  of  NT0TAL.  Since  the  71  coordinates  in  NPLUS (1, I, N)  are 
less  than  133,  KTIMES  is  reset  to  1  and  the  difference  is  reduced  to  62.  Next,  . 
the  DO  7030  loop  determines  that  .6  of  the  104  coordinates  in  the  interval  20  to 
40  inches  need  to  be  dropped.  Since  .6  times  20  equals  12,  the  points  to  be 
dropped  are  simply  assumed  to  be  located  in  the  interval  20  to  32  inches,  while 
all  other  points  in  the  interval  20  to  40  inches  are  considered  to  exceed  32. 

One  can  observe  in  Chart  7  that  the  calculated  value  of  AVRG  corresponds  to  the 
point  where  one  of  the  line  segments  joining  adjacent  points  crosses  the  .5 
probability  level. 

Once  the  mean  is  known,  the  next  step  is  to  relate  the  original  input  data 
to  a  new  set  of  20-inch  intervals  centered  about  the  mean.  Chart  8  illustrates 
how  this  is  done.  Each  NMINUS(J, I ,N)  and  NPLUS(J,I,N)  value  is  first  subdivided 
into  two  subelements  according  to  the  value  of  FRCTN.  For  example,  NMINUS (6 , I ,N) 
is  broken  up  into  2,  associated  with  the  12  inches  to  the  left,  and  1,  considered 
in  the  right  8  inches  of  the  interval  -100  to  -120  inches.  Subelements  from 
adjacent  intervals  are  then  paired  appropriately  and  added  to  get  NNEG(M)  and 
NP0S(M)  values,  where  the  index  M  is  associated  with  the  adjusted  set  of  intervals. 
Note  how  NPLUS (2 , I ,N) ,  which  involves  the  original  interval  containing  the 
estimated  mean  32,  contributes  62  to  NNEG(l)  and  42  to  NP0S(1).  The  cumulative 
totals  MCMLTN  and  NCMLTP,  where  NCMLTP  always  includes  the  maximum  value  of 
NCMLTN,  enable  the  computer  to  determine  when  all  the  original  nonzero  input 
information  has  been  processed.  The  interval  adjustment  just  described  is  done 
by  the  statements  beginning  with  JEND  =  100  -  KTIMES  that  follow  the  final 
determination  of  AVRG  in  the  DO  7030  loop,  or  in  the  DO  7110  loop. 

After  the  interval  adjustment  calculations  have  been  completed,  the  DO  7310 
loop  of  the  program  computes  the  values  in  the  NRFLCT  array  as  shown  in  Chart  9. 

The  NRFLCT  array  represents  an  equal  weight  combination  of  the  NP0S  and  NNEG 
data.  Any  NRFLCT(K)  value  indicates  how  often  the  absolute  values  of  the 
differences  between  the  vertical  coordinates  of  missing  first  rounds  and  the 
estimated  mean  are  equal  to  or  less  than  20  K  inches.  The  fractions  obtained 
when  one  divides  the  NRFLCT  values  by  NTOTAL  are  not  calculated  by  the  program, 
but  are  included  in  Chart  9.  Chart  10  illustrates  how  these  fractions  are 
related  to  the  distribution  of  concern.  ' 

The  DO  7330  loop  determines  the  three  alternative  estimates  for  the 
standard  deviation.  These  estimates  are  based  on  the  probabilities  associated 
in  a  normal  distribution  with  the  mean  plus  or  minus  0.5  standard  deviation, 
plus  or  minus  1.0  standard  deviation,  and  plus  or  minus  1.5  standard  deviation. 
Linear  interpolation  is  applied,  as  necessary,  to  the  distribution  implied  by 
the  NRFLCT  values  to  infer  estimates  of  0.5,  1.0,  and. 1.5  times  the  standard 
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CHART  6 


CALCULATION  OF  MEAN 


104  104 


TMINUS  =  102,  TPLUS  =  369. 

MTOTAL  =  471 
N5000  =  235 
DIFF  =  133. 

71.  <  133.  KTIMES  =  1 

DIFF  =  62. 

104.  >  62.  FRCTN  =  ,6 

AVRG  =32.0 
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PROBABILITY 


CHART  7  ESTIMATED  MEAN  (AVRG) 
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CHART  8  ADJUSTMENT  OF  INTERVALS 


N-3,  1=2 


7 

6 

5 


7 

8 


0 

3  x  .6 
14  x  .6 
29  ETC 

19  ' 

17 

20 
71 

104 
104 
58 
19  . 

11 

2 

0 


8 

7 

6 

5 

4 

3 

2 

1 

1 

2 

3 

4 

5 

6 


235 
233 
224 
pj  201 
=  178 
160 
141 
90 
339 
416 
o.  450 
S  465 

470 

471 


104 


NPOS  tK)  NNEG  (K) 


CHART  9  CALCULATION  OF  NRFLCT  AND 
CORRESPONDING  FRACTIONS 


8 

7 

6 


6 

7 

8 


2  235 

9  233 

23  q  224 

23  2  201 

18  |  178 

19  *  160 

51  141 


90 

90 

NRFLCT (K)/471 

104 

104 

194 

.412 

77 

181 

322 

.684 

34 

~  215 

^  375 

.796 

15 

w  230 

GO 

w  408 

i— 

.866 

5 

|  235 

2  436 

Od 

.926 

1 

z  236 

=  460 

.977 

0 

236 

469 

.996 

0 

236 

471 

1.000 

105 


•  POINTS  AS  IN  CHART  5 
■  POINTS  BASED  ON  NRFLCT  ARRAY 


ADJUSTED 


-140 


0 

COORDINATE  (INCHES) 


CHART  10  FRACTIONS  BASED  ON  NRFLCT  VALUES 
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PROBABILITY 


•  POINTS  AS  IN  CHART  5 
■  POINTS  BASED  ON  NRFLCT  ARRAY 
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CHART  11  ESTIMATE  OF  0,5  STANDARD  DEVIATION 
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MEAN  +  2  STANDARD  DEVIATIONS 


CHART  12  ILLUSTRATION  OF  DISTRIBUTION  FOR  MEAN  OF 

32  INCHES  AND  STANDARD  DEVIATION  OF  AO  INCHES 
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deviation.  For  example,  as  is  shown  in  Chart  11,  a  frequency  of  .383  is  asso¬ 
ciated  with  the  mean  plus  or  minus  0.5  standard  deviation.  Vertical  lines 
through  .5  -  (.383/2),  which  equals  about  .31,  and  through  .5  +  (.383/2), 
approximately  .69,  cross  line  segments  joining  points  based  on  the  NRFLCT  array 
at  intersections  that  correspond  to  adjusted  horizontal  scale  values  of  -18.6 
and  18.6  inches.  The  related  standard  deviation  estimate  is  37.2  inches. 
Similarly,  standard  deviation  estimates  based  on  1.0  and  1.5  times  the  standard 
deviation  turn  out  to  be  39.9  and  53.4  inches  respectively. 


3.  CONCLUSION.  Consider  again  the  context  within  which  arose  the  need  for 
the  procedure  explained  in  this  paper.  Chart  12  illustrates  the  information 
concerning  the  vertical  coordinates  of  missing  first  rounds  that  is  conveyed  by 
the  calculated  mean  and  a  standard  deviation  estimate  of  40  inches,  selected  as 
a  best  estimate  for  the  particular  situation  used  as  an  illustrative  example 
in  this  paper. 
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1400  Wilson  Blvd. 

Arlington,  VA  22209 
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1  Marine  Corps  Operations 

Analysis  Group 
ATTN:  Mr.  Love 
2000  North  Beauregard  St. 
Alexandria,  VA  22311 

1  Marine  Corps  Operations 

Analysis  Group 

ATTN:  Document  Ctr  (Dr.  Simmons 
2000  North  Beauregard  St. 
Alexandria,  VA  22311 

1  Commandant 

US  Coast  Guard 
ATTN:  G-0MR-2B/33 
2100  2nd  St. ,  S.W. 

Washington,  DC  20593 

1  Commander 

US  Army  Intelligence 
Analysis  Center 
ATTN:  IAX-OT-P  (Mr.  Alston) 
Arlington  Hall  Station 
Arlington,  VA  22212 

1  Commander 

Naval  Weapons  Center 
ATTN:  Code  40701 

(Mr.  Melvin  Keith) 

China  Lake,  CA  93555 

1  Reliability  Analysis  Center 

ATTN:  Mr.  I.L.  Krulac 
Griffiss  AFB,  NY  13441 

2  Chief 

Defense  Logistics  Studies 
Information  Exchange 
US  Arn\y  Logistics  Mgnt  Ctr 
ATTN:  DRXMC-D 
FT  Lee,  VA  23801 

1  Commander 

US  Army  Infantry  Center 
ATTN:  ATSH-CD-CS 
FT  Benning,  GA  31905 


No.  of 

Copies  Organization 

1  Commander 

US  Army  Combat  Developments 
Experimentation  Command 
ATTN:  Technical  Library 
FT  Ord,  CA  93941 

3  Director 

US  Army  TRADOC  Sys  Anal  Act 
ATTN:  ATAA-SA 

White  Sands  Msl  Rg,  NM  88002 

1  Director 

US  Army  TRADOC  Sys  Anal  Act 
ATTN:  ATAA-TF  (Mr.  Curl) 
White  Sands  Msl  Rg,  NM  88002 

1  Director 

US  Army  TRADOC  Sys  Anal  Act 
ATTN:  ATAA-TEA  (Mr.  Burnham) 
White  Sands  Msl  Rg,  NM  88002 

1  Director 

US  Army  TRADOC  Sys  Anal  Act 
ATTN:  ATAA-SL 

White  Sands  Msl  Rg,  NM  88002 

1  Director 

US  Army  TRADOC  Sys  Anal  Act 
ATTN:  ATAA-T 

White  Sands  Msl  Rg,  NM  88002 

1  Commander 

US  Army  Combined  Arms  Combat 
Developments  Activity 
ATTN:  Library 
FT  Leavenworth,  KS  66027 

1  Commander 

US  Army  Combined  Arms  Combat 
Developments  Activity 
ATTN:  ATCA-CFD 
FT  Leavenworth,  KS  66027 
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1  Commander 

US  Army  Combined  Arms  Combat 
Developments  Activity 
ATTN:  ATCA-CAM-V 
FT  Leavenworth,  KS  66027 

1  Commander 

US  Army  Combined  Arms  Training 
Developments  Activity 
ATTN:  ATZL-TDA-DSA 
FT  Leavenworth,  KS  66027 

1  Commander 

TRADOC  Combined  Arns  Test  Agcy 
ATTN:  ATCAT-TD-CON 
FT  Hood,  TX  76544 

1  Commander 

TRADOC  Combined  Arms  Test  Agcy 
ATTN:  Library 
FT  Hood,  TX  76544 

1  Commander 

US  Army  Mobility 

Equipment  R&D  Command 
ATTN:  ORDME-O  (B.  Briggs) 

FT  Bel  voir,  V A  22060 

1  Commandant 

US  Army  Field  Artillery 
School 

ATTN:  ATSF-CD-DA  (J.  Horn) 

FT  Sill ,  OK  73503 

1  Commandant 

US  Army  Field  Artillery 
School 

ATTN:  ATSF-CD-MS  (Mr.  Minton) 
FT  Sill ,  OK  73503 

5  Commander 

US  Army  Operation  Test  & 
Evaluation  Agency 
ATTN:  CSTE-ED 
5600  Columbia  Pike 
Falls  Church,  V A  22041 


1  Commander 

US  Army  Logistics  Eval  Agcy 
ATTN:  COL  Walker 
New  Cumberland  Army  Depot 
New  Cumberland,  PA  17070 

1  Commander 

US  Army  Foreign  Science 
&  Technology  Center 
ATTN:  DRXST-CA3  (Mr.  Bollendcrf) 
220  7th  Street,  NE 
Charlottesville,  VA  22901 

1  Commander-In-Chief 

Research  Analysis  Office 

Box  13,  J77 

ATTN:  Mike  Mulhern 

Camp  H.M.  Smith,  Hawaii  96861 

1  Commandant 

Command  &  General  Staff  College 
ATTN:  Technical  Library 
FT  Leavenworth,  KS  66027 

1  Commandant 

US  Army  War  College 
ATTN:  Library  R-1069 
Carlisle  Barracks,  PA  17013 

1  Commandant 

US  Arniy  Air  Defense  School 
Director  for  Combat  Developments 
ATTN:  Mr.  Don  Vance 
FT  Bliss,  TX  79916 

1  US  Army  Training  Systems  Center 

Training  Devices  Directorate 

ATTN:  ATTSC-TD 
FT  Eustis,  VA  23604 

1  Director 

US  Army  Electronics  Warfare  Lab 
ATTN:  DELEW-P 
FT  Monmouth,  NJ  07703 
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1  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Code  G-31  (Mr.  G.  Seidl) 
Dahlgren,  VA  22448 

1  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Code  DG-14  (Mr.  G.  Hornbaker) 
Dahlgren,  VA  22448 

1  Oklahoma  State  University 

Field  Office 
ATTN:  A.  Peebles 
PO  Box  1925 
Eglin  AFB,  FL  32542 

Aberdeen  Proving  Ground,  MD 

1  Director,  USABRL 

Bldg.  328 

1  Director,  USABRL 

ATTN:  DRDAR-TSB-S 
Bldg.  305 

1  Director,  USAHEL 

Bldg.  520 

1  Commander,  USATECOM 

ATTN:  DRSTE-CM-F 
Bldg.  314 


* 


115 


